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Zur Theorie des bei 
mit besonderer Berücksichtigung 
der Temperaturabhängigkeit von C, 


(Mit 7 Figuren) 


I. Einleitung 


Von Gerhard Damköhler 


Im Gebiete hoher Temperaturen treten bei der statistischen 
Berechnung des thermischen Verhaltens von Kristallen hauptsächlich 
zwei Schwierigkeiten auf; die eine ist mit dem Wesen des Kristalls 
verknüpft und besteht darin, daß sich die miteinander verkoppelten 
Schwingungen der einzelnen Teilchen wegen des anharmonischen 
Charakters nicht mehr wie bei tiefen Temperaturen in verschiedene 
Normalschwingungen separieren lassen. Die zweite Schwierigkeit ist 
charakteristisch für jede anharmonische Schwingung, auch für die 
einfachste, und beruht auf der speziellen, nicht mehr harmonischen 
Form der potentiellen Energie, die eine Auswertung der Zustands- 
summe nur noch nach Näherungsverfahren zuläßt, nicht aber in 
geschlossenen Ausdrücken. 


Bei der theoretischen Behandlung sind im wesentlichen zwei 
Wege eingeschlagen worden. Den einen gingen E. Griineisen’) und 
P. Debye?). Sie approximierten die anharmonischen Schwingungen 
der Gitterteilchen durch harmonische mit volumenabhängiger (und 
zwar nur volumenabhängiger) Frequenz. Die beiden genannten 
Schwierigkeiten waren damit formal umgangen. Der Haupterfolg der 
Griineisen-Debyeschen Theorie bestand darin, daß die Tempe- 
raturabhängigkeit des Ausdehnungskoeffizienten in einem weiten 
Temperaturbereich richtig wiedergegeben wurde. Für C, ergibt die 
Theorie allerdings stets den Dulong-Petitschen Wert. Doch ist 
dies Resultat sicher weniger für den Kristall charakteristisch als 
für die spezielle Art der Approximation, wie sich weiter unten noch 
zeigen wird. 

1) E. Griineisen, Ann. d. Phys. 39. S. 257. 1912. 

2) P. Debye, Wolfskehlkongreß, Göttingen 1913, Vorträge über die 
kinetische Theorie der Materie und Elektrizität. Teubner, Leipzig-Berlin. 1914. 
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Ein zweiter Weg ist von M. Born und E. Brody!) sowie 
E. Schrödinger?) beschritten worden. Sie berechneten direkt die 
freie Energie eines Systems miteinander verkoppelter anharmonischer 
Öszillatoren. Die zweite oben genannte Schwierigkeit bezwangen sie 
dadurch, daß derjenige Teil der Zustandssumme, der die anharmo- 
nischen Glieder der potentiellen Energie enthält, durch die ersten 
Glieder einer Potenzreihenentwicklung dargestellt wurde. Man bleibt 
damit allerdings im wesentlichen auf ein Temperaturgebiet beschränkt, 
in dem die anharmonischen Beträge der potentiellen Energie kleiner 
als kT sind. Die von den Autoren angegebenen Ausdrücke ent- 
sprechen jedoch nicht ohne weiteres der Helmholtzschen freien 
Energie des Gitters, da bei deren Berechnung nur diejenigen Schwin- 
gungen der Gitterteilchen berücksichtigt werden dürfen, die der 
Bedingung des konstanten äußeren Volumens entsprechen. Es läßt 
sich zeigen, daß die in den genannten Arbeiten für C, abgeleiteten 
Ausdrücke nicht für C, sondern für C, gelten. Am leichtesten über- 
sieht man dies am einfachen anharmonischen Oszillator, dessen 
potentielle Energie durch 


(1) =. +f xt 


gegeben sei. Die dafür von M. Born?) abgeleitete Beziehung  __ 


6f 
m? 
entspricht nämlich der spezifischen Wärme des unter dem äußeren 
Druck null schwingenden Oszillators, dessen Länge mit steigender 
Temperatur wächst, also sicher einem C,. 

Zusammenfassend ergibt sich demnach, daß die bisher vor- 
liegenden Theorien des festen Körpers bei hohen Temperaturen über 
den Verlauf von C,, entweder äußerst fragliche oder sogar falsche Aus- 
sagen machen. Vor kurzem hat nun A. Eucken‘) auf Grund quali- 
tativer Überlegungen die Vermutung ausgesprochen, daß C,, im klassi- 
schen Gebiet mit steigender Temperatur wieder abnehmen müsse, — 
was sich auch in einzelnen Fällen durch das Experiment wahrschein- 
lich machen ließ. 

Unter diesen Umständen erschien es angebracht, das thermische 
Verhalten fester Körper bei hohen Temperaturen erneut theoretisch 
zu untersuchen. In Abschn. II der nachfolgenden Arbeit soll zu- 


1) M. Born, E. Bordy, Z. f. Phys. 6. S. 152. 1921. 

2) E. Schrödinger, Z. f. Phys. 11. S. 170. 1922. 

3) M. Born, Geiger-Scheel, Hdb. d. Phys. XXIX, 2. S. 675. 1933. 
. Eucken u. W. Dannöhl, Z. f. Elektrochem. 40. S. 789. 1934. 
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nächst die Grüneisen-Debyesche Theorie erweitert werden. Es 
wird sich dabei zeigen, daß durch Hinzunahme einer reinen Tempe- 
raturabhängigkeit der Schwingungsfrequenzen Abweichungen von C, 
gegenüber dem Dulong-Petitschen Wert erklärt werden können. 
Eine exakte statistische Rechnung für den Fall der linearen Kette 
wird in den Abschnitten III und IV durchgeführt, und zwar in III 
unter Benutzung der Helmholtzschen freien Energie, in IV unter 
Benutzung des Gibbsschen thermodynamischen Potentials. Nach 
letzterer Methode werden in Abschnitt V verschiedene Potential- 
kurventypen behandelt, wobei auch Temperaturgebiete erfaßt werden 
können, in welchen die anharmonischen Energiebeiträge größer als 
kT sind. Auf eine exakte Behandlung des Kristalls selbst wurde 
vorerst wegen der großen Komplikationen, die sich bei der Rech- 
nung ergaben, verzichtet. Dennoch ist es wahrscheinlich, daß sich 
die an der linearen Kette gewonnenen Ergebnisse nicht nur quali- 
tativ sondern auch bis zu einem gewissen Grade quantitativ auf das 


Kristallgitter übertragen lassen. ae = 
II. Ausbau der Grüneisen-Debyeschen Theorie be: 


1. Die Grüneisen-Debyesche Theorie läßt sich am einfachsten, 
wenn auch nicht ganz der historischen Entwicklung entsprechend, 
folgendermaßen skizzieren: es wird angenommen, daß sich die gesamte 
innere Energie eines Kristallgitters zerlegen lasse in eine lediglich 
vom Volumen V abhängige Dehnungsenergie ®,(V) des schwingungs- 
freien Gitters und eine harmonische Schwingungsenergie, deren 
Frequenz vom Volumen (nicht dagegen unmittelbar von der Tempe- 
ratur) abhängen soll. Bezeichnet man mit N die Zahl der Teilchen 
im Gitter, mit v, die Grenzfrequenz des Debyeschen Spektrums 
und mit h,k, T wie üblich das Plancksche Wirkungsquantum, die 
Boltzmannsche Konstante und die absolute Temperatur, so erhält 


man für die Helmholtzsche freie Energie den Ausdruck): — re 
(3) | 0 A 
und bei hohen Temperaturen 
F =, —NET 


®,(V) sowie v, (V) wurden von Griineisen und Debye verschieden 
angesetzt: 


1) Vgl. M. Born, Atomtheorie des festen Zustandes, Teubner, Leipzig- 
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, (V) A 


(4a) Griineisen: 


(4b) Debye: 


5. Folge. Band 24. 1935 


B 


V, bezeichnet das Volumen am absoluten Nullpunkt. Auf den inter- 
essanten Versuch Grüneisens, das y der Formel (5a) aus dem 
Verlauf der Potentialkurven im Gitter abzuleiten, soll hier nicht 
eingegangen werden; die Beziehung (5b) entspricht einer formalen 
Taylorentwicklung. Aus Formel (3) und den Beziehungen (4) und (5) 
ergeben sich dann in üblicher Weise die thermische Zustandsglei- 
chung, Ausdrücke für den Ausdehnungskoeffizienten und die Kom- 
pressibilität, sowie für die innere Energie und die spezifische Wärme ( . 

2. Für die Grenzfrequenz v„ wollen wir nun, um bei hohen Tem- 
peraturen eventuelle Abweichungen von C, gegenüber dem Dulong- _ 
Petitschen Grenzwert beschreiben zu können, eine zusätzliche reine 
Temperaturabhängigkeit annehmen: 


(6) 


f(V) soll hier lediglich von V abhängen, g(T) lediglich von T. Natür- — 
lich ist dieser Ansatz nur näherungsweise richtig und gegebenenfalls 
auf der rechten Seite noch durch ein Glied von der Form h (V, T) 
zu ergänzen; doch soll davon vorerst abgesehen werden, weil einer- 


_ 


dln V 
Debye: (1 — + ( 
2 


(5a) 


(5b) 


= 


“poe 


0g (T) 
3°: 
T aT: 


seits h(V, T) vermutlich klein ist gegenüber f(V) und g(T), anderer- 
seits weil bei Berücksichtigung eines solchen Zusatzgliedes die End- __ 
formeln wesentlich komplizierter und unübersichtlicher werden. Aus(3) 
und (6) ergeben sich damit für die innere Energie U sowie für C, die _ 
Ausdrücke: 
aF 
M) U-=F-T(37),- 
69 (7) 
__°_ [ma 
(8) 0,=3Nkli ale 
In welcher Richtung C, bei hohen Temperaturen vom Dulong- | 
Petitschen Wert abweicht, hängt also von g(T) ab. Es läßt sich 
nun qualitativ entsprechend der Überlegung von A. Eucken (a. a. 0.) by 
einsehen, daß g(T) > 0 ist und mit steigender Temperatur monoton 
wächst. Denkt man sich ein Teilchen im Gitter schwingen, so wird | 
es bei einem bestimmten festgehaltenen äußeren Volumen J’ im Zeit- 
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mittel die Potentialkurve 0 (der Fig. 1) durchlaufen, die bei verschie- 
denen Temperaturen verschieden hoch ausgeschwungen wird. Kurve 0 
wird noch bei einer endlichen Amplitude z,, die vom Gittervolumen V 
abhängt, die potentielle Energie + 00 
liefern, da sich die einzelnen Gitter- 


. bd \ 
teilchen nicht durchdringen können. Im | _ 
Gegensatz hierzu wird jede harmonische \ Energee 
Parabel, mit welcher wir die Schwin- ‚_\ 
gungen auf Kurve 0 approximieren, \ 


erst bei z2=+m» die potentielle 
Energie + © besitzen, so daß wir bei 
höheren Temperaturen zur Approxi- 
mation steilere Parabeln verwenden 
müssen als bei tieferen (T, > T, > T,); 
dies bedeutet aber, daß v, bei fest- 
gehaltenem Volumen V mit steigender 
Temperatur wächst. Setzen wir quali- 
tativ an 


r n > Eu 

wo q und » positive Konstanten sein Somingungsanpivde Er 

sollen, so ergibt sich Fig. 1. Approximation der an- | 


(10) C,=3 Nk{l —n(n + 1)q Tr}, harmonischen Schwingung eines 
Gitterteilchens bei konstant ge- 

also ein Abfall von C, mit steigender haltenem Gittervolumen durch 
Temperatur. harmonische Schwingungen von 
3. Für die thermische Zustands- temperaturabhängiger Frequenz 


gleichungerhält man in bekannter Weise: 


(11) — +3NkT 


Olnv„ 
OV 


In diesem Ausdruck wird man stets in einem gewissen 'Tiemperatur- 
bereich, der eventuell experimentell näher abzugrenzen ist, formal 


6D, (V) PN = 
9 - 0 = — . 


(13) Iny,, Inv, ’o (r=) 7,(1 
Hier sollen x,, 7, 79, 7, konstante Größen sein, während y, (T) und 
g (T) zunächst noch unbestimmte Funktionen der Temperatur sind. 
Vermutlich ist 7,(T)<y,. Führt man noch die Abkürzung ein 

Q 3 NK x, an 


(14) 


j 
7 
ER j \ 
J \ “ts 
\ / 
\ / 
\ 
= » 
| \ \ / J 4 


6 Annalen der Physik. 5. Folge. Band 24. 1935 


so erhält man aus (11), (12) und (13) nach kurzer Rechnung für 
die Kompressibilität und ihre Druckabhängigkeit 
li) 


1 - on +- QT + 27%,» 


(15b) 


at 
Yo 
woraus sich im Fall 


e nM <ı der Ausdehnungskoeffizient zu 


Yo 


av 
(17) (3 r),= = )e- T 


ergibt. Man bekommt also auch unter De einer zu- 
sätzlichen reinen Temperaturfunktion g (T) in dem Ansatz (13) für », 
die schon von Grüneisen diskutierte Beziehung (17). Gl. (15a), die 
Grüneisen nicht in dieser Form angegeben hat, und Gl. 16h) ya 
überdies noch unabhängig von dem Glied (= )- 7,(T) in An- 
satz (13), so daß die Kompressibilität durch die abgeleitete Gl. (15a) 
auch dann noch richtig dargestellt werden dürfte, wenn der Aus- 
dehnungskoeffizient nicht mehr exakt durch (17) wiedergegeben 
werden kann. Diese relativ starke Unempfindlichkeit des Aus- 
drucks (15a) für die Kompressibilität gegenüber einer zusätzlichen 
Temperaturabhängigkeit der Frequenz v, ist beachtenswert, be- 
sonders im Hinblick auf die praktische Berechnung von c,— c, 
auf Grund experimenteller Daten. Hier ist man nämlich bei hohen 
Temperaturen, wo Kompressibilitätsmessungen zur Zeit noch gänzlich 
fehlen, auf eine Extrapolation der im mittleren Temperaturgebiet (meist 
30 bzw. 75° C) gewonnenen Kompressibilitätsdaten angewiesen, wozu 
jedoch eine Formel, die vermutlich auch noch im anharmonischen 
Schwingungsgebiet Gültigkeit besitzt, bisher nicht theoretisch ab- 
geleitet worden war. Durch Gl. (15a) werden die Grüneisenschen 
Kompressibilitätsmessungen an Cu, Fe und Pt zwischen 0 und 165° © 
richtig wiedergegeben !) 

Uber die Konstanten x,, 7, 7, 7, haben wir bisher keine 
theoretischen Annahmen gemacht, teilweise im Gegensatz zu 


= 


Griineisen, Ann. d. Phys. 33. S. 1239. 1910. 
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Griineisen und Debye. Vielmehr faßten wir die Entwicklungen (12) 
und (13) rein formal auf, die vier Konstanten sollten aus experimen- 
tellen Daten bestimmt werden. Dies scheint heute in der Tat der 
einzig mögliche Weg zu sein. Es zeigte sich nämlich, daß weder 
der Grüneisensche noch der Debyesche Wert für 1, die sich aus 
den Ansätzen (4a) und (4b) unmittelbar zu ee 


0 m+n+9 3,17 fir m=1, n=9 © 
Debye = VY; Grüneisen = = 3,50 9 


» M=3, n=! 
ergeben, mit den experimentellen Daten in Einklang sind. Wie man 
an Hand der Gl. (15b) sieht, scheidet der Debyesche »-Wert aus, 
da er die experimentell beobachtete Druckabhängigkeit der Kom- 
pressibilität nicht zu beschreiben "vermag. Aber auch der „-Wert 
Grüneisens ist offenbar nicht ausreichend, was schon aus der 
Zusammenstellung der experimentellen Daten bei diesem Autor!) 
hervorgeht. Eine zweite, allerdings nicht ganz so exakte Prüfungs- 
möglichkeit der 7-Werte besteht darin, daß man aus der Temperatur- 
abhängigkeit des Ausdehnungskoeffizienten, der näherungsweise durch 


Gl. (17) dargestellt sei, und (1 -Q ermittelt, aus der Tem- 


peraturabhängigkeit der Kompressibilität nach Formel (15a) x, und 


(20+ =| Q bestimmt, und schlieBlich durch Elimination von Q 
v 


und 7,/7, einen Näherungswert für 7, der mit 7’ bezeichnet sei, findet. 
Die so berechneten 7’-Werte sind mit den aus Gl. (15b) berechneten 
n-Werten in Tab. 1 verglichen. Man kann nicht verlangen, daß »/' 
genau mit 7 übereinstimmt, da bei hohen "Temperaturen die Dar- 
stellung des Ausdehnungskoeffizienten durch Gl. (17) meist unvoll- 
kommen ist und da ferner die experimentellen Daten für die 
Temperaturabhängigkeit der Kompressibilität zur Zeit noch ziemlich 
unsicher sind?. Doch sieht man, daß 7 größenordnungsmäßig den 
nach Gl. (15b) berechneten y-Werten näher kommt als den Grün- 
eisenschen Werten 3,17 bzw. 3,50. 


1) Geiger-Scheel, Handb. d. Phys. 10,1. 8. 36ff. (Artikel vonGrüneisen). 

2) Ein weiterer Grund dafür, daß die nach Gl. (15b) berechneten n-Werte 
mit den »'-Werten nicht genau übereinstimmen, könnte eventuell darin liegen, 
daß für die Temperaturabhängigkeit des Ausdehnungskoeffizienten sowie der 
Kompressibilität hauptsächlich die Anziehungsäste der Potentialkurven zwischen 
den Gitterteilchen maßgebend sind, für die Druckabhängigkeit der Kompressi- 
bilität hingegen vorwiegend die entsprechenden Abstoßungsäste. Es wäre denk- 
bar, daß in beiden Potentialbereichen zwar dieselben formalen Ansätze (12) 
und (13) gelten, daß aber die darin vorkommenden Konstanten etwas ver- 
schieden sind, 
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Die wichtige Frage, die mit den Ausgangspunkt zu der ganzen 2 
obigen Untersuchung bildete, nämlich, ob man aus der Temperatur- 
abhängigkeit des Ausdehnungskoeffizienten allein die der Kompres- 
sibilität berechnen könne, muß damit verneint werden, und zwar “4 
erstens, weil der Ausdehnungskoeffizient meistens nicht durch G1.(17) 
exakt darstellbar ist, und zweitens, weil selbst in dem Falle, 
wo Gl. (17) erfüllt wird, sich daraus nur die Größen Q und tine 


(n+ =2).Q bestimmen lassen, so daß die Berechnung der in 


9 
Gl. (15a) vorkommenden Größe (2 y+ sts) -Q nur möglich ist, wenn 
0 


(oder y,/7,) anderwärts bekannt sind. 
Es bleibt nun zu untersuchen, ob sich ein Grund dafiir ein- Eu 
sehen läßt, daß sowohl der Grüneisensche als auch der Debyesche rive 
Ansatz fiir 7 bzw. ®, (V) versagen. Zu diesem Zwecke vergegen- 
wärtigen wir uns noch einmal die Art der Approximation der an- a 


harmonischen Schwingung durch eine harmonische an Hand der 
inneren Energie des Gitters, wobei wir für die zeitlich periodisch 
veränderlichen Energieanteile die zeitlichen Mittelwerte (durch Über- m 
streichen gekennzeichnet) einführen. Es gilt stets und exakt: a ER 
(18) | U Ori. + Tot. Schwingung + Upot. Dehnung 
| U = Ükin. + (Opot. barm. Schwingung 7 Ükorrektion) + U pot. Dehnung 
U=3NkT + Ü gorrektion + Upot. Dehnung » 
Durch Kombination von (18) mit (7) folgt für „ bzw. ®, (V) die 
stets exakt gültige Gleichung 


(19) ®, (V) = Upot. Dehnung + Ü korrektion + 3NkT? (Fr) 
Vergleicht man hiermit die Formeln (4a) und (4b), so sieht man, 


daß den Ansätzen von Grüneisen und Debye die viel spezielleren 


Annahmen 
(20a) D, (V) Graneisen Upot. Dehnung » AG : 
(20b) D, (V vebye = harm. Anteil von Upot. Debnung 


zugrunde liegen, also Annahmen, die keineswegs als theoretisch 
ausreichend zu bezeichnen sind. Damit wird ein Versagen der 
Formeln (15) und (17) bei einer Rechnung mit den speziellen 
Grüneisen- bzw. Debyeschen 7-Werten verständlich. Gleichzeitig 
sieht man aber auch, daß die Ableitung irgendeines theoretischen 
Ansatzes für ®, (V), der natürlich Gl. (19) erfüllen müßte, mit den 
größten Schwierigkeiten verbunden ist. Die rein formale Behandlungs- 
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weise, wie sie oben durchgeführt wurde, erscheint vorläufig als die 
einzig mögliche. 

4. Eine Neuberechnung von C,—C, wurde für einige Sub- 
stanzen durchgeführt. In Fig. 2 ist jeweils der Verlauf von C,, 
wie er heute am sichersten erscheint, angegeben und der daraus 
errechnete Verlauf von C,. 

Die C,-Werte wurden weitgehend der Zusammenstellung von 
Kelley') entnommen. Die von ihm angegebenen Formeln wurden 


* 


3 


Y 


Fig. 2. Verlauf der spezifischen Wärmen C, und C, für einige Stoffe 
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an Hand der Originalliteratur kritisch gesichtet und konnten im 
Falle des Cu, Ag, NaCl und KBr ohne weiteres verwandt werden: 


Cy + 146-1079 T, 
cr Ag: C, = 5,60 + 1,50-10-* 7, 
Nal: C, = 10,79 + 420. 10T, 
€, = 11,49 + 3,60 - T. 


1) K.K. Kelley, U.S.Dep. of Commerce Bur. of Mines 371, Washington 
1934. Vgl. auch Chem. Zentralbl. 1934, II, 3359. 
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Beim Pb wurden auch direkt die Werte von Klinkhardt') benutzt, 
die insbesondere zwischen 150°C und dem Schmelzpunkt sicherer 
zu sein scheinen als die der anderen Autoren. Beim KCl fand 
die Formel 


C,= 11,18 + 3,76 10-3 T 


Verwendung, die konstant um 0,2 cal/°C höhere Werte liefert, als von 
Kelley offenbar versehentlich angegeben wurde. Die Formel für 
KBr, die sich auf die Magnusschen Messungen ?) bis 270° C gründet, 
wurde versuchsweise bis zum Schmelzpunkt hinauf benutzt, was 
bei der im allgemeinen linearen Temperaturabhängigkeit von C, 
nicht allzu fehlerhaft sein kann. 

Bei der Berechnung von C, —C, wurden die von A. Eucken 
und W. Dannöhl (a. a. O.) bestimmten Ausdehnungskoeffizienten 
benutzt. Den Kompressibilitäten wurden die von P. W. Bridgman 
bzw. J. C. Slater (a. a. O.) bei 30 und 75° C gefundenen Werte 
zugrunde gelegt; aus ihnen wurden die zwei Konstanten x, und 


(24 + en ( der Formel (15a) ermittelt und damit dann die Kom- 
pressibilitäten bei höheren Temperaturen berechnet. 

Wie aus Fig. 2 ersichtlich ist, nimmt C, beim NaCl, KCI, 
KBr und Ag mit steigender Temperatur wieder ab, während beim 
Cu und Pb eine geringe Zunahme bestehen bleibt. Es ist denkbar, 
daß dieser letztere Umstand auf die Metallelektronen zurückzuführen 
ist, doch bedarf es hier zur endgültigen Klärung sicher noch weiterer 
experimenteller wie theoretischer Arbeiten. 


~ 
III. Die Helmholtzsche freie Energie der linearen Kette u 
bei hohen Temperaturen on 


In Abschnitt II hatten wir qualitativ gesehen, ob und in welcher 
Richtung Abweichungen von C, gegenüber dem Dulong-Petit- 
schen Gesetz bei hohen Temperaturen zu erwarten sind. Quanti- 
tativ ließ sich allerdings aus der Theorie allein ohne Zuhilfenahme 
experimenteller Daten nichts aussagen; vielmehr mußten die formal 
eingeführten Konstanten aus jenen bestimmt werden. Um auch 
quantitative Aussagen auf Grund lediglich theoretischer Überlegungen 
zu erhalten, mußte man versuchen (etwa in Anlehnung an die schon 
zitierten Arbeiten von Born und Brody bzw. Schrödinger), ein 
System miteinander verkoppelter anharmonischer Oszillatoren exakt 
zu behandeln, allerdings, wie schon eingangs betont wurde, unter 


1) H. Klinkhardt, Ann. d. Phys. 84. S. 167. 192 
2) A. Magnus, Physikal. Ztschr. 14. 8. 5. 1913. 
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Berücksichtigung der Bedingung des gegebenen äußeren Kristall- 
volumens. Nachdem sich jedoch eine derartige theoretische Be- 
handlung des Kristallgitters selbst als äußerst kompliziert erwiesen 
hatte, sei vorerst die statistische Berechnung der linearen Atomkette 
gegeben. 

1. Definitionen: Die lineare Kette bestehe aus 2 N Teilchen, 
die durch fortlaufende Indizes gekennzeichnet seien. Die Lage- 
koordinate des jten Teilchens sei r , bezogen auf einen allen Teilchen 
gemeinsamen Koordinatenanfangspunkt. r, kann man sich zu- 
sammengesetzt denken aus einem Betrag L,, der die Lage des 
Teilchens im Zeitmittel angibt und nur von der gesamten Ketten- 
länge % abhängt, und einem Betrag z,, der die jeweilige Schwin- 
gungsamplitude des Teilchens, gerechnet von seiner Mittellage aus, 
darstellt. Also 
(21) 


Der Abstand zwischen dem j ten und kten Teilchen betrage 


(22) r— = 7, =(L,-L)+@,— 2) = 
Dabei gilt insbesondere 

2: -* 
on 


Fiir die kinetische Energie der Kette folgt dann: 


Die potentielle Energie sei durch eine Potenzreihenentwicklung 
dargestellt: 
7 (2 ” 
2 


Mit den Abkürzungen: 


a 2 4 
3 
} (‘ 
( 
4 
| 
= 
°. 
| 
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1 | i 1 73 
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1,1% 8) i 
(26d) B,=z1- U9+ L, Ps 
gelangt man zu 
ot. D jk 
- 
| + 0, + + E 
Das erste Glied auf der rechten Seite entspricht der poten- . 54 


tiellen Dehnungsenergie der Kette, die folgenden der potentiellen 
Schwingungsenergie. Der Wertbereich der hierin auftretenden x, 
ist nun durch die Bedingung der konstanten gegebenen Kettenlänge 
einzuschränken. Man berücksichtigt jene am einfachsten durch die 
Annahme der zyklischen Kette), d.h. durch die Annahme, daß die 
Amplitude sich nach jeweils 2N Teilchen periodisch wiederholt. 
Damit ergibt sich zwanglos die lineare Transformation 


+N 
(28) = > wärs wobei p,= 
Ersetzt man auf diese Weise in Gl. (27) die Koordinaten x. durch 
die y,, so verschwinden die in y, linearen Glieder. Die Glieder 


zweiter Ordnung in y, lassen sich durch die lineare Transformation an 


Y, = 12Z,, wobei 2,=2_, und 
in eine tll von Quadraten verwandeln. Schließlich erhält man Er u 
auf diesem Wege für die potentielle Energie den Ausdruck cr ~ 
(30) = Up+ A+B+6 
mit den Abkiirzungen 
N j,k 
(la) A= wobei A, > 84; A., sin? 
p=0 
Pi PaPs - 
B= a worin nur die Kombinationen p, + p, +7, = 10, 
| zu berücksichtigen sind, 
(31b) B > > zZ.) 
123 = (21 2 23 — = 29 23 — % 
Bus” (sin g, + sin 9, (j—k) + sing, (jh), 
1) Vgl. hierzu M. Born, Theorie des festen Zustandes (a. a. O., S.587f. 


Bs. 
| i 
“ 
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C= > Ciog,, worin nur die Kombinationen 


2 N e weg 
-2N 


+ 2, z,+ 2, 2,232, 


) | % % 2, 25 2,)} 
Js d 
+ cos — k) + cos — ky] 


“y + 608 (p,+ —k) + cos (p, + 
+ cos (p,+ + cos + — k) 
+ cos (py + —k) +008 k). 


ses! a 2. Die Zustandssumme läßt sich im Gebiete hoher Temperaturen 
zerlegen in einen kinetischen Anteil, in einen. Anteil, der von der 
potentiellen Dehnungsenergie herrührt, und einen solchen, der der 


potentiellen Schwingungsenergie entspricht. 


‘ 


POL. Schwingg. 


(33) = /m, di... ekl = const TY, 

j | 


(35) > = 141 [ de, 


Die Determinante 4 ist in bekannter Weise aus den Ab- 
leitungen der x, nach den z, baw. 2, gebildet. Entwickelt man in 
Gl. (35) den anharmonischen Anteil des Integranden in einer Potenz- 
reihe und vernachlässigt man Glieder von höherer Ordnung als 22 
bzw. Ö, so folgt: 


| 


u pot. Schwingg. 


Bei der weiteren Rechnung, die nur kurz angedeutet sei, ist zu berück- 
sichtigen, daß bei der Integration des Ausdruckes (36) nur solche Glieder der }} 
nicht verschwinden, die die z, bzw. %, in gerader Potenz enthalten. Da ferner 
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7 | 
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Kette angreift, folgen: 
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unter dem Integranden der Grad der Mehrfachsummen von B*? und C [man vgl. 
die Formeln (31a) bis (31e)] jeweils durch die Anzahl der darin auftretenden 
voneinander verschiedenen |p|-Werte gegeben ist, so brauchen bei hinreichend 
großer Kettenlänge nur diejenigen Glieder berücksichtigt zu werden, die eine 
Höchstzahl voneinander verschiedener |p|-Werte aufweisen, d. h. nur diejenigen 
Glieder, die z, bzw. x, quadratisch enthalten. Damit ergibt sich: 


1 Ag 
—0o TT 
p=0 
+00 _4 1,2 + 00 1,2 rm» 
Y = (x k (k 1 
(38) e 12-1-2 N 
* A, 
TT A, 
oo — oo p=0 


1,2, 3 
Dabei bedeuten >* die Summation über alle Kombinationen, die der 
ON 


2 1,2 
Bedingung p, + p, + ps = 0 genügen, >, die Doppelsumme über p, und 
-2N 
“Py, welche unabhängig voneinander alle Werte zwischen — N und + N an- 
nehmen sollen. 


Durch Einsetzen dieser Ausdrücke findet man schließlich: 


x 
PO. Schwing. | N 
a (39) ‘ Ze p=0 


1.43 L2G | 
P 123 12-1-2 
| 4 - AN, 
3. Zur Berechnung der freien Energie führen wir noch die 
Abkürzung ein 
und erhalten 
N 
(41) F=U,„—2NkThmT+kT > In k T)?-W—k T- const. 
p=v 
In diesem Ausdruck hängen die W, sowie W lediglich von den 
Potentialkurven zwischen den einzelnen Kettenteilchen ab und von 
der Kettenlänge 2. Für die spezifische Wärme Cy, bei konstanter 
Kettenlänge sowie für die Druckkraft 8, die an den Enden der 


oT? 


(42) oll W). 


(43) 
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Die Formeln vereinfachen sich für den speziellen Fall, daß 
jeweils nur die Potentiale ® zwischen benachbarten Teilchen wesent- 
lich sind. Für die Gl. (26) schreiben wir dann: 


1 3) 2 
5 M3) + 13 LL, C= 18 
wobei L = 


(45) 


so folgt 
(46) 
und nach einfacher Rechnung 


(47) Cy = fi + 
hig 


Ferner gilt für & = 0 in dieser Näherung: 


0g 


so daß man für Cg, die spezifische Wärme unter konstantem Druck, : 
unter Benutzung von 


(ar) 

oT 

(50) Cg —Cyg=T-— a2 

1) Gl. (48) und (49) besagen, daß die mittlere Kettenlänge gleich ist der 


mittleren Länge des 2-atomigen Moleküls multipliziert mit der Zahl der Ketten- 
glieder. Dies ist zunächst nicht selbstverständlich, da man sich vorstellen 


könnte, daB während der Schwingung in jedem Augenblick die Ausdehnung | 
einzelner Kettenglieder durch die gleichzeitige Verkürzung anderer teilweise __ 
kompensiert würde, so daß die mittlere Länge eines Kettengliedes kleiner _ 
wäre als die des 2-atomigen Moleküls. Daß dies, solange man Kräfte zwischen 


übernächsten und entfernteren Teilchen vernachlässigt, nicht der Fall ist, ergibt 
sich zwangsläufig aus der obigen Ableitung sowie den Entwieklungen des 
Abschn. IV. An anderer Stelle ‘wird gezeigt werden, wie sich der bei hohen 
Temperaturen anharmonische Schwingungszustand aus den bei niedrigen 
Temperaturen vorhandenen Normalschwingungen entwickelt, und welche Bei 

träge diese einzelnen Normalschwingungen zur thermischen Ausdehnung liefern. 
Auch hier ergibt sich, wie mir Herr Prof. Eucken freundlich mitteilte, für 
die mittlere Länge eines Kettengliedes die des 2-atomigen Moleküls. | 
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erhält: 
K 
12 4 f r 


Man überzeugt sich leicht, daß der Ausdruck (51) fiir Cg identisch 
wird mit der Formel (2) in der Einleitung, wenn man setzt 


0 0 
mw? pw 6a 24f 
und damit 
1 (5 7)_ 1 (5 1 _ _.6f 
\ 12 mo!:\12 mtw* 4 mo!) m’w® m? 


so daß die dortige Behauptung, daß Gl. (2) nicht C, sondern C, 
darstelle, zu Recht besteht. Ferner ergibt sich aus Formel (47), 
3 

‘ wohl für alle Potentialkurven, wie wir noch in Abschn. V sehen 
werden, erfüllt ist. 

Die hier benutzte Methode zur Berechnung der freien Energie 
ist grundsätzlich auch beim Kristall anwendbar, sofern man die 
dort noch komplizierteren Rechnungen nicht scheut. Allerdings 
läßt sich die Methode stets nur in einem beschränkten Temperatur- 
gebiet praktisch durchführen, auch bei der linearen Kette, da zwei 
wesentliche Annahmen gemacht wurden: sowohl die potentielle 
Energie als auch die anharmonischen Glieder der Zustandssumme 


daß C, mit steigender Temperatur sinkt, wenn 


(im obigen Falle ¢ #7 G+2)) sollten sich durch die ersten Glieder 
von Taylorentwicklungen wiedergeben lassen. Die zweite Forderung 
bedeutet jedoch wegen der schlechten Konvergenz von e-*, daß 

U ‚nharm. =B C<k T. 

Über das thermische Verhalten bei noch höheren Temperaturen, 
man möchte sagen in dem „eigentlich anharmonischen Temperatur- 
gebiet“, können wir also vorerst nichts aussagen. Wie man auch 
hier weiterkommen kann, wird im nächsten Abschnitt gezeigt. 


IV. Das Gibbssche thermodynamische Potential der linearen Kette 
bei Vernachlässigung der Kräfte zwischen übernächsten 
und weiter entfernten Teilchen | 
Eine exakte Separierung der miteinander verkoppelten an- 
harmonischen Schwingungen ist, wie schon eingangs bemerkt wurde, 
beim Kristall nicht möglich. Dies gilt auch für die lineare Kette. 
Vernachlässigt man jedoch hier die Kräfte zwischen übernächsten 
und weiter entfernten Teilchen, so wird die Zahl der für die 
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potentielle Energie maßgeblichen Teilchenabstände gleich der Zahl der 
voneinander unabhängigen Lagekoordinaten der Kette, und man 
kann deren potentielle Energie in einzelne anharmonische Potentiale 
separieren, wenn man die Abstände jeweils benachbarter Teilchen 
als neue Lagekoordinaten wählt. 

Wir betrachten eine lineare Kette von 2N Teilchen, die ab- 
wechselnd die Massen m und wi besitzen sollen. Bezeichnen wir 
mit N, und N, die Zahl der Teilchen in der Impulszelle 4, mit n 
die Zahl der gleichen Abstände L, zweier benachbarter Teilchen 
und mit @,, @,, 8, die entsprechenden a priori-Wahrscheinlichkeiten, 
so erhält man für die thermodynamische Wahrscheinlichkeit der 


Führt man die Abkürzungen ein 
N; N, 


und nimmt man wie üblich an 

u =. 

(54) 
=8,=8, é 


so ergibt sich für die ntropie 


| S=k IN [In « Zum] +N 3 
aN 2 g, ing, |} 


wobei &, bzw. &, die kinetische Energie eines Teilchens bezeichnen, 
das sich in der Impulszelle A befindet, und ®, jeweils die potentielle 
Energie zweier benachbarter Teilchen von dem Abstande L,, An 
den ‚Enden der Ketten greifen die Druckkraft % an. Für eine 
virtuelle Verrückung folgt dann aus der thermodynamischen Gleich- 
gewichtsbedingung 


= 
= 
2 
= 
© 
E U=N > Ja + N > +2N > g,®,, 
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und den Gleichungen 
A 


‘ad 
(59) > 4f,=9, Dan-=0, S'dy,=0 


in bekannter Weise 
kT kT 


a 
Setzt man diese Ausdrücke in (55), (56) und (57) ein, so erhält man 
für das Gibbssche thermodynamische Potential der Kette 
(61) G=U—TS+8&-2=— 2NkTiln(# - o) + const.}, 


wenn wir definieren 


e 


const = Inf. 


s und o stellen nichts anderes dar, als den kinetischen bzw. poten- 
tiellen Anteil der Zustandssumme eines einzelnen anharmonischen 
Oszillators, der unter dem äußeren Druck $ schwingt. Sie ist, wenn 
auch nicht in geschlossener Form, so doch durch brauchbare Nähe- 
rungsverfahren, stets auswertbar, und zwar nicht nur durch Potenz- 
reihenentwicklungen, die bei hohen Temperaturen z. T. sehr schlecht 
konvergieren, sondern auch graphisch, so daß auch das „eigentlich 
anharmonische Temperaturgebiet* erfaßt werden kann. Hierin ist 
die Methode des Gibbsschen thermodynamischen Potentials der 
Methode der Helmholtzschen freien Energie, wie sie in Abschn. III 
entwickelt wurde, überlegen. Während die Gibbs sche Methode, die an 
und für sich nur schwache Verkopplung der einzelnen Oszillatoren bei 
der linearen Kette voll ausnützt, wird jene bei der Helmholtz schen 
Methode durch die hier notwendige Annahme der zyklischen Kette 
künstlich verstärkt [man vergleiche die Ausdrücke (31) für B und (). 
Aus (61) erhält man in bekannter Weise die Beziehungen 

ag ag ag ag 

| 
RK fal 7’ 

so daß man Cy nach Formel (50) berechnen kann. Entwickelt 
man ® und o in Potenzreihen und berücksichtigt man nur die 
ersten Glieder derselben, so erhält man wiederum die Ausdrücke (47) 
bis (51), wie man sich leicht überzeugt. 


(63) 


= 
Br; 
| 
(62) af = 
4 
m 
A > 
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V. Diskussion des thermischen Verhaltens der linearen Kette 
(insbesondere von Cy) unter der Annahme spezieller Potentialkurven') 


Sieht man zunächst von dem „eigentlichen anharmonischen 
Temperaturgebiet“ ab, so lassen sich die Verhältnisse am schnellsten 
an Hand der Formeln (47)—(51) überblicken. Sie wurden für die 
verschiedensten Potentialkurven ausgewertet und die Ergebnisse sind 
in Tab. 2 wiedergegeben. Man sieht, daß Cg ausnahmslos mit 
steigender Temperatur sinkt. 

Eine genauere Untersuchung wurde für das Morsesche sowie 


das Coulombsche Potential durchgeführt, die man bis zu einem 
r Kurve 1 y = e~ 
J sols 
-1 

e-wli+e 


PER 
10 (1 204022) 


+e * +032) 


6y=e 

‚-10(1 +e 


= 
-10(1 +. *+ 042) 


Ty=e 


0 @ 02. 0 02 0 @ M % 1% 


r sd Fig. 3. Kurvenblatt zur graphischen Bestimmung der Zustandssumme 
einer linearen Kette im Falle des Morseschen Potentialansatzes 


a ®=Dl|ı+e e -2e e | zwischen den einzelnen Teil- 


chen. Die Temperatur ist bestimmt durch nd = 10, die Druckkriifte 


0, 0,1, 02, 0,3, 0,4, 0,5. 


dure D 
gewissen Grade als Repräsentanten der unpolaren bzw. polaren 
Bindung ansehen kann. o wurde hier bis zu den höchsten T’empe- 
raturen hinauf graphisch bestimmt. 

Ei Man zeichnet für verschiedene Werte von 7 und & die Kurven 


- (®(L) + RL) - 
e sowie € me 
1) Die Kräfte zwischen übernächsten und weiter entfernten Teilchen 
sind auch weiterhin vernachlässigt, sofern nichts Besonderes vermerkt ist. 
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und mißt die von ihnen und der Abszissenachse eingeschlossenen Flächen aus. 
In Fig. 3 ist ein solches Kurvenblatt, das sich auf eine bestimmte Temperatur 
und verschiedene Drucke bezieht, wiedergegeben. Auf diesem Wege erhält 
man o als Funktion von T und & sowie o,,.., (7), den potentiellen Anteil 
der Zustandssumme für die entsprechende harmonische Schwingung. Man 
kann dann schreiben 


(64) o(T,8) = Oparm. (T) 


m. (T) (T, O) 


o(T, &) 
o (T, O) 
gegen ft auf, wie es in den Figg. 4a und 4b geschehen ist, und kann aus dem 
(T, 
o (T, O) 


| fiir 8 = 0 ablesen. Das allmähliche Flacherwerden 


Zur Auswertung der letzten beiden Beziehungen trägt man In 


Abschnitt auf der Ordinatenachse bzw. der Anfangssteigung 3. I 
0° | o (T, &) 


bzw. n 


Kompressibilität mit steigendem Druck. 


Im Falle des Morseschen Potentials 
( 2(L Le) 
68 ®=D\l+e ® = e 
(68) N, + 


erhält man auf graphischem Wege die bis zu den höchsten Tempe- 
raturen gültigen Formeln 


| = (0,53 + 0,03) 


‘harm. 

wobei die jeweils angegebenen Unsicherheiten bedingt sind durch 

die Fehler bei den besprochenen graphischen Integrationen. Die 

Übereinstimmung mit den Werten der Tab. 2 ist verhältnismäßig 

gut. Allerdings steigt Cg in dem eigentlich anharmonischen Tempe- 
1 


raturgebiet stärker an als nach Gl. (51) und = (3 Te ist etwas 


o (T, O) o (T, 


und erhilt auf Grund der Gleichungen ”* für ® = 0 
( 

(65) Cy = {1 

aT | Sharm. (T) 
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Fig. 4 a. der 
Zustandssumme der linearen Atom- 


kette (Morsesches Potential) 


Fig. 4b. Druckabhängigkeit der 
Zustandssumme einer Ionenkette 
(Born-Mayersches Potential, Be- 
rücksichtigung der Kräfte zwischen 


entfernteren Teilchen) 
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größer als nach Gl. (48). Interessant ist, daß der Ausdehnungs- 
koeftizient und die Kompressibilität in dem gesamten untersuchten 


Temperaturgebiet bs T = . (D = Dissoziationsenergie) nicht zu- 


nimmt, jedenfalls nicht außerhalb der angegebenen Fehler der 
Methode. 

Bei der Ionenkette wurde zunächst daran gedacht, ebenfalls die 
Kräfte zwischen übernächsten und weiter entfernten Teilchen zu ver- 
nachlässigen. Als Potential wurde angesetzt 

-L 


+ Bee + const, 


(69) D= 


worin sich B aus | or) =0 zu B= Ar ee bestimmt und o 
OL Jı=1, 


nach M. Born und J. E. Mayer!) für die Alkalihalogenide etwa 
0,345 AE beträgt. Zur Beurteilung der Anharmonizität dieses 
Potentials wurde es mit dem Morseschen Potential in Fig. 5 ver- 
glichen. Es sind gleichzeitig die Schmiegungsparabeln der ent- 
sprechenden harmonischen Schwingungen eingezeichnet und ferner 
diejenigen Energien durch Niveaulinien gekennzeichnet, die bei be- 


Wess 


2 


t 


= - 

.4 


+ stimmten Temperaturen praktisch gerade noch erreicht werden 
I p 5 

— Pmax 

[es wurde = 7,5 gewählt, woe * 5,5-.10”* wird]. Man 
u _ erkennt, daß im Temperaturgebiet bis 1200° abs. die Anharmoni- 
Ri. j’ zität der Kurve II (Coulombpotential), gemessen an dem Verhältnis 


von anharmonischer zu harmonischer Schwingungsamplitude, praktisch 
übereinstimmt mit der Anharmonizität der Kurve J im Tempe- 


raturgebiet bis = 103). 

Für beide Potentialkurven ist also hier ein gleichartiger Verlauf 
der thermischen Ausdehnung vorhanden, was bedeuten würde, daß 
bis 1200° abs. auch der Ausdehnungskoeffizient der Ionenkette 
praktisch konstant wäre. Dies ist nun aber höchst unwahrscheinlich, 
denn der Ionenkristall ist gerade dadurch ausgezeichnet, daß er im 
Vergleich mit anderen festen Stoffen (z. B. Metallen) nicht nur einen 
besonders großen Ausdehnungskoeffizienten, sondern auch eine be- 


sonders große Temperaturabhängigkeit desselben besitzt; es ist nicht 
anzunehmen, daß sich hier Kette und Kristall so wesentlich von- 


einander unterscheiden. 
1) M. Born u. J. E. Mayer, Ztschr. f. Phys. 75. S. 1. 1932. sa pe 
2) Die Temperaturen 600 bzw. 900 bzw. 1200° abs. bei Kurve // würden 
= 0,043 bzw. 0,064 bzw. 0,085 auf Kurve J entsprech 
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| ‘ae der Behandlung der Jonenkette dürfte also die Berück- 
sichtigung der Kräfte zwischen entfernteren Teilchen unerläßlich sein. 
Um die in Abschn. IV abgeleiteten Beziehungen auch weiterhin benutzen 
zu können, müssen wir allerdings die potentielle Energie so an- 
setzen, daß sie in separierbarer Form von den Abständen nur un- 
mittelbar benachbarter Teilchen abhängt. Dies ist zwar nicht exakt, 


\ / 
Vergleich der Anharmonizität \\ Fy —- 


verschiedener Potentialkurven: . \ / 


Kurve I: Morsesches Poten- N Z 


tial, Kurve II: Born-Mayer- 
sches Potential unter alleiniger 
Berücksichtigung der Kräfte 
zwischen aufeinanderfolgen- 


denTeilchen, Kurve//I:Born- — 1200° abs 


Mayersches Potential unter \\ 600 
Berücksichtigung auch über- 
greifender Coulomb-Krifte \ 

\ 


aber doch näherungsweise möglich. Wir denken uns hierzu die 
Kette zunächst in einem Zustand, in welchem alle Abstände 
zwischen aufeinanderfolgenden Teilchen L, betragen, dann dehnen 
wir an einer einzigen Stelle von L, auf L (ohne die anderen Ab- 
stände zu ändern) und fragen, welche Coulombsche Energie dabei 
zu leisten ist. Wie man an Hand der Fig. 6 leicht übersieht, 
ergibt sich 


1) Eigentlich sollte diejenige Coulombsche Energie berechnet werden, 
die aufzuwenden ist, wenn in einer sich im thermischen Gleichgewicht be- 
findenden Kette ein einziger Abstand von L, auf L gedehnt wird, während 
alle anderen Kettenglieder ihre mittlere Länge L behalten. An Stelle von L 
wurde oben der Einfachheit halber L, gesetzt. u 


iw 
4 : 
u. 
> \ pa Ze 
‘on 
4 \ 
a 


Annalen der Physik. 5. Folge. Band 24. 1935 


5 A A A A 
A A 
+ “ther * 


A 
2L,+L + 


Shel 


A 
+ 
A A 


A A 
~ 8h ıL 5 


A 
5L, 
A 
+ 4L 
A 
Ly 


2A A 
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Fiigt man noch ein AbstoBungsglied hinzu, so erhilt man fiir den 
Beitrag, den ein einzelnes Kettenglied effektiv zur potentiellen 
Energie der Kette liefert 


A x (L — L,) (- 1)" 


wobei sich B aus (Sr) 


BE” 0 bestimmt zu B= ; 


Das Potential (71) ist fiir den Fall / = 0,12 (für 9 = 0,345 AE 


ist dann L, = 2,87 AK. = Gitterkonstante A NaCl-Kristalls) in Fig. 5 
als Eures. III gezeichnet. Sie ist viel anharmonischer als Kurve II, 
was man übrigens auch schon auf Grund der Fig. 6 erwarten muß: 
Die Arbeit, die hier bei der Dehnung (Verkürzung des Abstandes 
zwischen 4 und 5 gegen die Anziehung (Abstoßung) dieser beiden 
Ionen zu leisten ist, wird zum Teil durch freiwerdende (auf- 
zuwendende) Abstoßungsenergien zwischen dem linken und dem 
rechten Teil der Kette, insbesondere durch die AbstoBungsenergien 
zwischen 4 und 6 sowie 3 und 5, kompensiert (vergrößert). 
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Das Ergebnis der Auswertung auf Gund der Formeln (47) bis (51) 
ist in Tab. 2 enthalten. Das Resultat der graphischen Auswertung 
zeigt Fig. 7. Cg nimmt bis zu den höchsten Temperaturen prak- 
tisch linear zu. Allerdings ist der Anstieg übereinstimmend mit 
dem Befunde beim nn Potential etwas og als nach 


‘harm 


(97) 


y Cg Ca 
harm 

(al 50) 
— | AT 


a 02 
Fig. 7. Spezifische Wärmen Cy und Cg, a 
_ Ausdehnungskoeffizient sowie Kompressibilität der Ionenkette int 


Formel (51). Cg nimmt mit steigender Temperatur immer stärker ab, 


so daß die Berechnung nach Formel (47) im „eigentlich anharmo- 
nischen Temperaturgebiet* viel zu kleine Abweichungen vom Dulong- 


‚ben würde. 1 (22) 
Petitschen Wert ergeben würde, & \ar)s und (ö8)r 


nehmen mit steigender Temperatur zu, (Fr). ‚deutlich stärker als 
linear und wesentlich stärker als ee (se). beides in qualitativer 
0 


Ubereinstimmung mit den Befunden an Kristallen. a 
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(Juantitativ ist allerdings die Zunahme von 2 (27), bei der 
Xo 
lonenkette erheblich geringer als beim lonenkristall, nämlich 


etwa 20°/, bei der Kette gegenüber 100°/, beim NaCl-Kristall in 
dem Temperaturgebiet von 300—1100° abs. Hier besteht eine Analogie 
zur Atomkette, wo wir unter Zugrundelegung des Morseschen 
Potentialansatzes überhaupt keine Temperaturabhingigkeit des Aus- 
dehnungskoeffizienten gefunden hatten. Möglicherweise liegt hier ein 
für die Kette besonders charakteristisches Verhalten vor. Es könnte 
dadurch bedingt sein, daß die Zahl der nächsten Nachbarn eines 
Teilchens in der Kette geringer ist als im Kristall, so daß, ähnlich 
wie beim Übergang von Ansatz (69) zu Ansatz (71), die effektive 
Anharmonizität beim Kristall größer würde als bei der Kette. 
Doch ist sicher noch ein anderer Umstand für die geringe Tempe- 
raturabhängigkeit der bei der Kette berechneten Ausdehnungs- 
koeffizienten verantwortlich: die Abstoßungsglieder der verwandten 
Potentialansätze steigen nämlich mit sinkendem Teilchenabstand 
langsamer an, als es bei materiellen Teilchen der Fall ist. Man 
erkennt dies insbesondere daran, daß die Zustandssumme o(T, &), 
nachdem sie mit steigendem & zunächst gesunken war, bei größeren 
g-Werten wieder zunimmt (vgl. die Figg. 3 sowie 4a und 4b), was 


in Wirklichkeit jedoch deshalb nicht möglich ist, weil oe 2,0 
sleic 3 = DNk 
gleichbedeutend ist mit ), 2NkT <0. 
Zusammenfassung 


1. Es wird die Griineisen-Debyesche Theorie des festen 
Körpers durch die Annahme einer zusätzlichen reinen 'Temperatur- 
abhängigkeit der Schwingungsfrequenz formal erweitert. Der sich 
daraus ergebende Ausdruck für die spezifische Wärme C, fordert 
einen Wiederabfall derselben mit steigender Temperatur. 

2. Von einer zusätzlichen reinen Temperaturabhingigkeit der 
Schwingungsfrequenz ist in der erweiterten Grüneisen-Debye- 
schen Theorie der Ausdehnungskoeffizient nur wenig, die Kom- 
pressibilität sogar praktisch unabhängig. Man gewinnt damit für 
letztere eine theoretisch wohl begründete Formel, welche die bei 
mittleren Temperaturen gemessenen Kompressibilitäten bis in 
das anharmonische Temperaturgebiet hinein zu extrapolieren ge- 
stattet. Eine direkte Berechnung der Temperaturabhängigkeit der 
Kompressibilität aus der des Ausdehnungskoeffizienten ist nicht 
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3. Es wird gezeigt, daß die in der ursprünglichen Grüneisen- 
Debyeschen Theorie auftretende Dehnungsenergie des schwingungs- 
freien Gitters im anharmonischen Temperaturgebiet noch durch 
Korrektionsausdrücke zu ergänzen ist, die von der Anharmonizität 
der Schwingungen herrühren. 


4. Es wird der Temperaturverlauf von C, beim NaCl, KCI, 


KBr, Ag, Cu und Pb aus den experimentellen Daten für C,, 


I / dv 1 (dv 
(sr), und unter Benutzung der abgeleiteten 


Formel für die Kompressibilitit neu berechnet. Mit Ausnahme 
von Cu und Pb wird mit steigender Temperatur stets ein Wieder- 
abfall von C, gefunden. Das abweichende Verhalten der 
Metalle dürfte wahrscheinlich auf die Metallelektronen zurück- 
zuführen sein. 

5. Es wird für die lineare Kette der Ausdruck der Helm- 
holtzschen freien Energie (als Funktion von Temperatur und 
Kettenlänge) abgeleitet. Die dabei verwandte Methode, die von 
dem Gedanken der zyklischen Kette Gebrauch macht, ist allerdings 
praktisch nur durchführbar in einem Temperaturgebiet, in dem der 
anharmonische Anteil der potentiellen Energie Dgnharm< kT. Die 
Formeln für die spezifischen Wärmen bei konstantem Druck 8 
bzw. konstanter Kettenlänge 2 Cg bzw. Cy, sowie die Ausdrücke für 
2 (#7 ı und z (FE), sind im Falle, daß die Kräfte zwischen über- 
nächsten und entfernteren Teilchen vernachlässigt werden können, 
angegeben. Die mittlere Ausdehnung der Kette wird unter diesen 
Umständen gleich der des entsprechenden 2-atomigen Moleküls 
multipliziert mit der Zahl der Kettenglieder. 


v 


6. Es wird das Gibbssche thermodynamische Potential einer 
linearen Kette (als Funktion von Temperatur und Druck) unter Ver- 
nachlässigung der Kräfte zwischen entfernteren Teilchen abgeleitet. 
Durch graphische Auswertung der Zustandssumme ist dabei eine 
theoretische Behandlung bis zu den höchsten Temperaturen möglich. 


7. Das thermische Verhalten der linearen Kette wird in einem 
Temperaturgebiet Pyunarm. <kT für die verschiedensten Potential- 
kurven untersucht, wobei stets die spezifische Wärme bei konstanter 
 Kettenlänge Cg mit steigender Temperatur wieder abfällt. Die 
Atomkette, wo zwischen benachbarten Teilchen das Morsesche 
Potential gelten soll, sowie die Ionenkette werden bis zu den 
höchsten Temperaturen hinauf graphisch untersucht. Um bei 
letzterer die vermutliche Temperaturabhängigkeit des Ausdehnungs- 
koeffizienten darzustellen, mußten auch die Goulombschen Kräfte 
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zwischen übernächsten und entfernteren Teilchen berücksichtigt 
werden. Es ist dadurch wahrscheinlich gemacht, daß auf diese 
übergreifenden Kräfte die besonders große Tremperaturabhängigkeit 
der Ausdehnungskoeffizienten bei Ionenkristallen zurückzuführen ist. 

8. Es ergeben sich Anhaltspunkte dafür, daß bei hohen, aber 
experimentell noch durchaus zugänglichen Temperaturen die Ab- 
stoBungsglieder des Morseschen, sowie des Born-Mayerschen 
Potentialansatzes den wahren Potentialverlauf zwischen materiellen 
Teilchen nur mehr unvollständig wiedergeben. 


Herrn Professor Dr. A. Eucken danke ich herzlich für die An- 
regung dieser Arbeit, sowie die interessanten Diskussionen während 
der Durchführung. 


Göttingen, Phys.-Chemisches Institut der Universität. 
(Eingegangen 13. Juli 1935) 
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Erweiterung einer früheren Formel 


für die Erdabsorption 
in der drahtlosen Telegraphie Dr ae 


Von K. F. Niessen 


[Natuurkundig Laboratorium der N. V. Philips Gloelampenfabrieken ig 
Eindhoven (Holland)} 


(Mit 8 Figuren) 


1. Einleitung und Ergebnisse 


In einer vorigen Arbeit!) berechneten wir den von einer ebenen 
Erde absorbierten Teil T der Strahlung eines vertikalen Dipolsenders 
in der Höhe h oberhalb der Erdoberfläche. Die allgemeine Formel 
für T ließ sich zu einer ganz einfachen Formel reduzieren, wenn 
noch damals näher präzisierte Annahmen über Bodenart, Wellen- 
länge und Dipolhöhe erfüllt waren. 

Sind «, und o, die Dielektrizitätskonstante und das Leit- 
vermögen der Erde, so daß der für die Erde charakteristische Para- 
meter %,? in der Differentialgleichung 


A Il, + k,? I, = 0 
für die Hertzsche Vektorfunktion JJ, wird: 


& + 
(1) k,? = 


und gilt für die Luft (Index 1) 


(2) 
so ließen die erwähnten Annahmen sich so formulieren: 


6, > & 0°. 


ky 
y/*>1. 


1 
h 
0< 


Nach der Annahme (3) sollte der Leitungsstrom groß in bezug 
auf den Verschiebungssirom und nach (4) der Absolutwert des 


1) K. F. Niessen, Ann. d. Phys. [5] 22. S. 162. 1935; zitiert mit 1. 
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Brechungsindex groß in bezug auf die Einheit sein. Durch (5) — 
wurden Dipolhöhen, die zu groß waren in bezug auf die Wellen- 
länge A, ausgeschlossen. 

Wir stellten damals fest, daß für zwei Größen A und B 2 er 
nur dann 


gelten würde, wenn % 
| 
(6) = 10 
Mittels dieser Beschränkungen ergab sich dann für den in der 5 


Log n + U (e) 1 
(4) f= 


U (ce) = 0,107 — + Log a, 


> 
- 


3 a? 


V (a) = 0,307 „sine _ cos 


(8) 

= 0341-0217 
Y (a) = 0,217 (Cé(a) + em, 


Far 

4ah 

(9) “= . 
A 


Mit Log. wurden immer Briggsche Logarithmen angedeutet 
(im Gegensatz zu lg für die natürlichen). 

Wir wollen nun in der vorliegenden Arbeit die Annahme (3) 
fallen lassen, d.h. wir wollen auch den Fall von Leitungs- und 
Verschiebungsstrom gleicher Größenordnung in Betracht ziehen, also 
voraussetzen: 


Aa 
0<r<z: 


x= bedeutet: Verschiebungsstrom vernachlässigbar, 
= 0 bedeutet: Leitungsstrom vernachlässigbar. 
Im allgemeinen wird: 
0, 
(11) 


1) Die damalige Funktion W («) wurde hier absichtlich zerlegt in 3 


W (a) = X (a) + Ye). 
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Für den in (10) für den allgemeinen Fall definierten Absolut- 
wert n des Brechungsindex gilt nun ein anderer Ausdruck als 


früher. 
Wiederum wollen wir n als groß in bezug auf die Einheit a 
voraussetzen: 


Es erweist sich weiter als bequem, die frühere Beschränkung 
(5) jetzt zu ändern in: 


= 


(3) 

wozu 

(14) 

genügt. 

Für = > gehen die neuen linie in die früheren 
über 


pe den jetzigen finden wir nun als endgültiges Resultat: =» 
= 


V2 cos (Logn + U(a)) Ki) 
(15) T= 

n Va) + 2sin X(c) + V 2cos Y (a) 
mit U(«), V(«), X(a), Y(«) wie in (8), 


(16) K(y) = 0,307 (cos 4 — sin {+ + (1 + cotg 


Diese Formel, welche im Anhang abgeleitet wird, geht offenbar x 


für 7 = > wiederum in die frühere Formel (7) über. 


Sind die Annahmen, mit der Formel (15) abgeleitet wurde, 
nicht erfüllt, so können wir wohl die absorbierten und die durch die 
Luft hinweg ausgestrahlten Teile (T und 1 — T) der Gesamtenergie 
des Senders berechnen, aber müssen dann kompliziertere Formeln 
benutzen, die wir unterwegs bei der Ableitung der Formel (15) an- 
geben werden. Vgl. die Formeln (34) und (36). 


ao 2. Anwendung mit numerischen Größen aus der Praxis 
Führen wir das in CGS-Einheiten gemessene Leitvermögen o 

ein, d. h. setzen wir 

(17) =4 
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und bedenken wir, daß 

o= 


so schreiben sich die früheren Annahmen 8) Pr (4) wegen der 
Voraussetzung (6) als: 


100 (A in ¢ M) 
(20) A> 2co 


Mit den von Zenneck?) angegebenen und auch von Sommer- 
feld?) benutzten Werten für «, und a: 7 


Seewasser . . . 3=80, o =107" (Fall I), 
SiBwasser . . . ,=80, o= 107" (Fall ID, 
Nasser Grund. . ,=10, a=5-10-1% (Fall IID, 
Trockner Grund . ,= 4, o= 10715 (Fall IV) 


bekamen wir dann früher als untere Grenze der Wellenlänge 7* 
. 
(in m): 


(19) nach (20) 
Seewasser 13 > 1,6 
Süßwasser. . . 4 13000 1600 Q*in m). 
Nasser Grund. . 330 | “ 
Trockner Grund . 7000 A* >17000 


Da in der Praxis keine Wellen länger als 20 kM benutzt 
werden, nahmen wir überall (d.h. für jede Bodenart) 


20000, 
so daß die in Betracht kommenden Wellenlängenbereiche waren: 
Seewasser . . . 13 < 4*< 20000 in m) 
Süßwasser . . . 13000 < 20000 
Nasser Grund. . 330 < A* < 20000 T-Formel 
Trockner Grund . 17000 < 4* < 20000 


Da wir uns aber jetzt nicht mehr die Beschränkung (3) oder 


u 7 (19) auferlegt haben, so zerfällt die erste Spalte in (21). 

7. u Die zweite Beschränkung, jetzt in der Form (12) lautet nach 
4 a, (17) und (18) und bei andermaliger Benutzung der Definition (6) 
— für relativ kleine Größen: 

> 10* — 6,7 

2co 


1) J. Zenneck, Ann.d. Phys. 23. S.859. 1907. 
2) A.Sommerfeld, Ann. d. Phys. 28. S. 721. 1909. 
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Diese zweite Beschränkung in ihrer jetzigen Form (22) ist praktisch ct 2 
identisch mit der früheren Form (20) nur im Falle von nassem oder 4 “| 
trocknem Grund (¢, = 10 bzw. 4), aber nicht mehr im Falle von a 


Wasser (&, = 80). a 

Statt der zweiten Spalte in (21) finden wir jetzt in den Fällen I . a = 

> 1, 1000, 330, 17000. 


Die für die neue Formel, d. h. für (15), in Betracht kommenden 


Wellenlängenbereiche sind also: 
Süßwasser: 1000 < 2* < 20000 (A* in m) 

Nasser Grund: 330 < 4* < 20000 | für die neue Formel. 


Trockner Grund: 17000 < 4* < 20000 


Für Wasser als zweites Medium wurde der Wellenlängenbereich 
also nach der Seite von kleineren Wellen hin erweitert, während 
die T-Formel selber nur wenig komplizierter aufzubauen war, wie 
ein Vergleich zwischen (15) und (7) zeigt. 


Wie steht es nun mit der früher und mit der jetzt erlaub- 
ten Dipolhöhe? 


was bedeutete: 
für Seewasser: 


Süßwasser: ice < | (h* und 4* in m) 
GES fiir die 
Nassen Grund: < 400 frühere Formel. 


Trocknen Grund: h* < 

| 
Im erweiterten Falle der vorliegenden Arbeit haben wir aber 

nach (14), (6) und (12) 


2% 


- 
>. 
Früher hatten wir nach (5) und (6): 
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d. h. für: 
Seewasser: h* < /6400 + 3600 


7 * 
Süßwasser: < V6400 + 0,0036 und 2* in m) 
3 | fiir die 
Nassen Grund: h* < 100 7 a neue Formel. 


Trocknen Grund: h* < 7000 + Ie 
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4 10 5 20 
Wellenlänge in km 

mb Fig. 1 

en . a Die höchsterlaubten Dipolhöhen findet man in den Fig. 1, 2, 3 

angegeben'), wo die früher erlaubten Maximalwerte punktiert ein- 


getragen sind. Für nassen und trocknen Grund lassen die alte 
und die neue Formel dieselben Maximalwerte zu. 


1) Die Figg. 2 und 3 geben einige Details aus Fig. 1 vergrößert wieder. 
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Bei Anwendung der erweiterten T-Formel kann also für See- 
oder Süßwasser als zweites Medium nicht nur der Bereich für die 
10 


4 
‘7 


Dipolhöhe in m (h*) 


9 


4 


See 


4 


‘asser 


10 15 20 
Wellenlänge in m 
Fig. 2 


Dipolhöhe in km 
in 


Wellenlänge in km 
Fig. 3 
Wellenlänge, sondern auch derjenige für die Dipolhöhe vergrößert 
werden. 
Die vergrößerte Leistungsfähigkeit der T-Formel erfordert aber, 


wenn wir sie wiederum graphisch darstellen wollen, eine größere 
Mannigfaltigkeit von Kurven. 
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Früher hatten wir „vorläufig“ Kurzwellensender ausgeschlossen 
und gleichfalls Sender in Flugzeugen und fanden dann für den in (9) 
definierten Parameter a: 


0<a<3, 


so daß wir nur T(n,«)-Kurven zeichneten mit den Parameterwerten 
a= _0,25, 0,50, 1 und 2. 

ies des jetzt, speziell für Wasser, vergrößerten Maximal- 
wertes für h und des verkleinerten Minimalwertes von A lohnt es 
sich jetzt, auch den Fall der Kurzwellensender in unsere Betrach- 


Wir schließen aus: 


und den in (23) angegebenen Maximalwerten von h*: 


für Seewasser: max = 6400 + 3600 > 1) 


Süßwasser: max = V6400 + 0,0036 > 1000) 


Trocknen Grund: @max = 37 yar (a* > 17000). 


Betrachten wir Kurzwellensender über Seewasser, d. h. auf einem 
Schiff oder in einem Flugzeug, so brauchen wir, damit die neue 
T-Formel anwendbar sei: 


für A*~ 20:h* < T,a< 4, 
ss für 4* ~ 100:h* < 80, « < 9,6. 


| Die Notwendigkeit, jetzt auch &-Werte bis zu « ~ 10 betrachten 
TER ~ zu müssen, führte uns zu einer Abbildung der U(a)-, V (a)-, X («)- 
und Y («)-Funktionen für 0< <8 (vgl. Fig. 4), welche beim Zeichnen 


der T-Kurven mit «= 1, 2,3....8 zu benutzen sein werden. 
Aber wir müssen jetzt eigentlich sprechen von ; 


T(n, a, z)-Kurven , 


wo xy nicht mehr den früheren Wert 2/2 zu haben braucht, aber 
prinzipiell zwischen 2/2 und 0 liegen kann, was natürlich eine 
zweite Vermehrung der Rechenarbeit zur Darstellung der T-Formel 
bedeutet. 

Die numerischen Verhältnisse in der Praxis liegen so, daß wir 
(wenigstens im Gültigkeitsgebiet der T-Formel) nicht bis zu 7 = 0 
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herunter zu gehen brauchen. Wir haben nämlich nach (11), (17) 
und (18) für den minimalen y-Wert TR 


und also nach (22) 


/ 100 \2 


Xnin = 37° 


so daB fiir Wasser 


der Minimalwert eigentlich entstand). | 
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Fig. 6 


Um auch T-Kurven zeichnen zu können fir 0 < 7 < F haben 


in Fig. 4 auch noch K(z) als Funktion von y angegeben, so daß 
‘ig. 4 die sämtlichen für Formel (15) benötigten Funktionen enthält. 
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Im Anschluß an die vorige Arbeit zeigen wir nun in Fig. 5 
wie z. B. die T-Kurve für «= 2 sich bei einem kleiner werdenden 
Argument x des Brechungsindex verhält, d.h. wenn das Verhältnis 
zwischen Leitungsstrom und Verschiebungsstrom verkleinert wird. 

Wir bekommen aber ein übersichtlicheres Bild der Formel (15), 
wenn wir nicht n als Variable und « und 7 als Parameter wählen, 


10 T 


| 


a | 


Fig. 7 


sondern « als Variable und n und y als Parameter. Für den Para- 


meter 7 haben wir nur den früheren Wert y = = (praktisch nur 


Bie 


Leitungsstrom) und die neuen Werte 7 = 0 (praktisch nur Verschie- 


‘ 


bungsstrom) und = T gewählt, Figg. 6, 7 und 8. Die T-Kurven 


wurden voll ausgezeichnet für diejenigen «-Werte, die in Über- 
einstimmung mit den Annahmen sind, auf die Formel (15) basiert —_ 
wurde. 

Da wir hierzu @ = 1 oder größer voraussetzten, gerade wie in 
der Ableitung der früheren T-Formel, wurden die T-Kurven in dem 


Gebiet <a< punktiert. 
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Nach der Bedingung (13) sollte 
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1 
«<4n- — 
8 V2 cos 


sein, und wenn Größenordnungen durch Faktoren 10 unterschieden 


werden, wie in (6), so bedeutet das er a 
< 2 > 


a7 


& 


“ss. fell NefS |L-- 12 
| 
er 1 2 3 4 E 6 7 8 9 u 
1g. 
ws Hätten wir aber die Größenordnung von « durch Faktoren 5 statt 


Faktoren 10 gekennzeichnet, so würde der zugelassene «-Wert zwei- 
mal größer sein können. Auf dieses Übergangsgebiet beziehen sich 
die rechten punktierten Teile der T-Kurven. 

Zum Gebrauch der T-Kurven erinnern wir noch daran, daß 
nach (12), (17) und 


T gibt die in Pi Erde absorbierte Energie, 1- -T die i in sali Luft — 
ausgestrahlte Energie an, beide ausgedrückt in der vom Sender aus- 


gestrahlten Energie als Einheit. ee, 
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iat. 3..Ableitung der Formel (15) 


fiir den Bruchteil 7 der absorbierten Energie 
_ Diese Ableitung möchten wir nur kurz angeben, da sie an vielen 
Stellen mit der Ableitung der früheren Formel parallel geht. 
Wir denken uns (wie früher) eine Ebene in der Höhe 


4=h+e(e>0, e>0)') 

oberhalb der Erdoberfläche angebracht, welche gerade den Dipol 
oben streift und betrachten die durch diese Ebene nach oben 
| gehende Strahlung 

Diese besteht sowohl aus direkt nach oben emittierter, wie aus 
von der Erde „reflektierter“ Strahleng. Auch denken wir uns eine 
parallele Ebene 3= h—« gerade unterhalb des Dipols angebracht, 
und betrachten die hierdurch nach unten gehende Strahlung: 

6, 
d.h. die Differenz der unmittelbar nach unten emittierten und der 
von der Erde „reflektierten“ Strahlung. 

Die in I, (17a) und (19a), für ©,-„;. und ©,.,_, angegebenen 
Formeln waren dort berechnet unter der alleinigen Annahme «, = 1, 
6, =0 d.h. nicht leitender Luft und können also auch hier noch 
benutzt werden. Da bei diesen Annahmen k, positiv reell wird, 


= 


e wurde damals geschrieben k, = p und wir haben also [vgl. (2)] 
mit: 
f 1 1 k,? ef 2hY 
i=p 
(25) ® i[3( 


ji 
= 
2 
l= 
1) Verwechslung dieser kleinen Länge & mit der Dielektrizitätskonstante s, 


der Erde wird wohl nicht zu befürchten sein. 
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wo {% bedeutet „Imaginärteil von“ und |---| den Absolutwert 
angibt. 

Für 0<l<p ist Yp®— und für oo ist — 
immer positiv gedacht. 

Weiter ist wegen des komplexen Wertes von k,? zu be- 


denken, daß 


immer so zu bestimmen ist, daß sein Realteil positiv ist. Formel (24) 
läßt sich nach der Transformation 


l= p yl — u? = 
schreiben in der Form: | = 
ec? p® = u-A-jB 
(27) 1— wu? O<uct) 
wo A und B reell sind. 
Formel (25) wird bei der Transformation 
py! + u: 
jah-.ı = lu + A +7 B)?| 
| not, | 
(1 + u*)udu , 


wenn gesetzt wird: 


(29) O<uco). 
wo C und D reell sind. 

Die Größen A, B, C und D sind Funktionen von % und die 
Vorzeichen in ihren Definitionsgleichungen (27) und (29) wurden so 
gewählt, daß sie positive reelle Größen werden. Letzteres folgt 


sofort aus dem Umstand, daß TER 


(30) V?—k?= py Fu? — 


einen positiven Realteil haben sollte. 

Zur näheren Berechnung von A, B, C und D müssen wir nun 

den Ausdruck (10) für k,?/k,? benutzen: 
2 
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n? sin 


wobei M, N, w und wurden bei der Substitution: 


] 1 = Me-iv (0<u<l) 


((<u< oO), 
so daß 


(32) 0<u<{l) 
Da k,?/p? im ersten Quadranten liegt und ye —k,? immer einen 


positiven Realteil haben sollte, muß mit Rücksicht ‘oul (30) gewählt 


| 


Bis. jetzt haben noch keine PERS stattgefunden (außer 
der Annahme einer nichtleitenden Luft) und es berechnet sich also 


noch strenge als Quotient von (28) und (26); a aes 
(u — u?) du + _ Du+wv) d 
(34) (u + A) + B? fers “ws OF+ 


u*) || cosau+ — - BY + du 
( (u + A)? + 


Diese komplizierte Formel ist für jede Bodenart bei beliebiger 
Wellenlänge und Dipolhöhe anwendbar. Wie im Anfang erwähnt, 
haben wir für Wasser und nassen Grund: we 


Deshalb kann fiir diese Bodenarten in M [siehe ( (31)], wo N 
u? neben &, fortgelassen werden, und auch in (32) für yw, wo auch 
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Für N und % dagegen handelt es sich um einen Bereich von u 
von 0 bis oo, so daß dort eine Vernachlässigung von u? neben :, 
an und für sich nicht gestattet sein würde. N und @ aber waren 
für C und D benötigt und diese treten nur in einem Integral auf, 
dessen Integrad den Faktor e- «“ enthält. 

Wenn wir also, trotzdem u bis oo gehen kann, in N und gu? 
neben «, vernachlässigen, werden die dadurch entstandenen Fehler 
(speziell beträchtlich bei u>1) doch durch die nihilierende Wirkung 
des Faktors 


au 


a = 1 (oder größer) !) 


im Falle 


zum Verschwinden gebracht. 
Nach obiger Vernachlässigung ergibt sich dann: 


tgy=tp= 
Wir haben aber nach (11) 


0, 
& @ 


tgy = 


und da w und g nach (33) im zweiten Quadranten liegen mußten, 
haben wir also zu nehmen: 


In dieser Näherung sind die Größen A, B,C und D, die — 
Funktionen von wu waren, in Renstenten entartet: 


cos 4 
n 
(35) 1 
B=C =—-—sin<. 


Nach Substitution von (35) in (34) ergibt sich fiir den Teil T 
der Gesamtenergie, der absorbiert wird, falls Wasser oder nasser Grund 
als Erde gewählt wird, bei beliebiger Wellenlänge und nicht allzu 
kleiner Dipolhöhe h*). 


(36) T= WEDER in... = 


1) Da a = 4ah/2 ist, liegt meist der Fall a» 1 (oder größer) auch wirk- 
lich vor, denn A ist die Höhe des’ Dipols und kann bei wirklichen Antennen 
am besten als die Höhe des „Schwerpunktes‘ der vertikalen Antenne oberhalb 
der Erdoberfläche interpretiert werden. 


: oe 2) Wir setzten ja a = 42h/A ~1 oder größer an. 
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wobei unter Benutzung der Abkürzungen: 


1 
p ah 
ı 7 25) Wey utd’ 
0 


2 


2(u? + 7, u +0) 


0 
1 
P, = — sin? Stau, 
i 26 
P, = — u)du. 


Der Term 3 im Nenner von (36) entsteht auf ähnliche Weise 
wie in I, S. 141. 

Nur P, und P, lassen sich in einfacher Weise berechnen. 
Deshalb ziehen wir nun noch zur weiteren Vereinfachung die Be- 
dingung (12) in Betracht, d.h. 


Nach der Formel 
(37) n= Var+ Ver+ 4c* o72?, 


(wo 4 in Zentimeter) ist für Wasser (e, = 80) die Bedingung n> 1 
praktisch immer erfüllt, unabhängig davon, ob wir mit See- oder Süß- 
wasser, mit langen oder mit kurzen Wellen zu tun haben. Bei nassem 
Grund aber muß man sich jedesmal noch davon überzeugen, daß die 
benutzten Wellen lang genug sind, umn>1, z.B. n >10 zu machen, 
damit die nachfolgende T-Formel anwendbar sei, sonst muß man sich 
an Formel (34) halten, die ja für jede Bodenart gültig war. 

Die Vereinfachungen von P, und P, sind bei Hinzunahme der 
Bedingung n> 1 basiert auf die Kleinheit der Größen y,, 7, und 0. 
Wir bedürfen für die benötigten Transformationen und Vernachlässi- 
gungen wohl nach ähnlichen in I verweisen, wo y, und y, einander 
gleich waren. 

Während früher zur weiteren Vereinfachung neben n> 1 auch noch 
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anzunehmen war, wollen wir jetzt neben n>1 "auch noch “65 
=x V2co <1 4 
voraussetzen, woraus die Bedingung (13) entstand. - 
Der Bedingung ie ve jedenfalls genügt, wenn wir wählen 
h 
was die früher Bedingung (14) war. 
Wir finden nun bei nicht zu kleiner Wellenlänge [in (37) sollte vie 
jan>1 und bei nicht zu kleiner und nicht zu großer Dipolhihe 
(wir wählten ja « = er ~ 1 oder größer und in 1497<75): ER: 
P, = „eos (— 5 +lgn—Zcotg 5) 
1 
P, = — cos Flign —Ig e+ (5 
(wo 7 = 0,577..., Eulers Konstante) en 
1 COS « sin « 
P,=z+ a2 a3 > 


P,= — cos lien da 


Werden noch Si(@) und eingeführt: 


und die auf S.146 (I) für die Eulersche Konstante abgeleitete 
Formel benutzt, so kommt nach Substitution der gefundenen Aus- 
drücke für P, bis P, in (36) und bei Benutzung von Briggschen 
statt natürlichen Logarithmen sofort die in (15) erwähnte T-Formel, 


deren Gültigkeitsbereich schon im Anfang besprochen wurde. 7 
Eindhoven, den 17. Juni 1935. 


(Eingegangen 3. August 1935 
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Zur Verwendung von Operatoren 
in der Kontinuumsdynamik BE 


Von Hellmut Baumann 


Inhalt: Einleitung. — Vorbemerkungen. — A. Der Widerstands- 

operator: Herleitung; Die mechanische Widerstandslehre; Die unendliche ; 
Kette; Der Widerstand eines Kontinuums. — B. Der Wellenoperator: Her- 2 
leitung; Diskussion. — C. Der Reflexions- und der Deformationsoperator: a Ahn 
Einfluß eines angekoppelten Systems auf das Wellenbild; Einfluß mehrerer 
angekoppelter Systeme auf das Wellenbild. — D. Beispiele: Berechnung 
spezieller Balkenwellen; Bestimmung der Eigenfrequenzen einer Klaviersaite te « 
während des Anschlags. — Zusammenfassung. BR: 
Nachdem im Anschluß an Heaviside’ in der Elektrizitits- ar : 
lehre schon häufig die Operatorrechnung Verwendung fand, hat man am, 
sich in letzter Zeit auch auf dem Gebiet der Mechanik gelegentlich = 
dieser Methode bedient, und zwar besonders bei der Berechnung er: 
der kleinen Schwingungen materialkundlich komplizierterer Körper. 1 
Die meisten Autoren nach Heaviside haben sich der funktinen- 


theoretischen Seite dieses Kalküls gewidmet. In dieser Hinsicht va, 
erscheint er als eine „symbolische Methode“, d.h. eine Schemati- 
sierung der durch die Laplacetransformation und das Fourierintegral 30 
gegebenen Transformation, und somit schließlich als eine Integral- 
sammlung % 5), 

Demgegenüber soll der Hauptinhalt der vorliegenden Arbeit = 


darauf gerichtet sein, die allgemeinen Eigenschaften dieser Trans- N Se 
formation, die sich weitgehend mit denen der linearen Differential- ia 
gleichungen mit konstanten Koeffizienten decken, zu einem Kalkül en oN 
zu vereinigen, der hier speziell auf die Mechanik zugeschnitten wird. 
1) O. Heaviside, Electromagnetic Theory, London 1893—1912. shes « 
3 2) H. Jeffreys, Proceedings Royal Society London (A) 138. S. 283.192. == 
“4a 3) Vgl. auch W. Prager, Ing.-Archiv 4. S. 16ff. 1932, u 


4) J. R. Carson, Electric Circuit Theory and the Operational Calculus, 

_ New York 1926, erweiterte deutsche Bearbeitung von F. Ollendorf u. K. Pohl- ay 

hausen, Elektrische Ausgleichsvorgänge und Operatorenrechnung, Berlin 1929. 

5) B. van der Pol, Phil. Mag. [7] 8. S. 861ff. 1929. De vs 
Annalen der Physik. 5. Folge. 24. 4 on er 
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Die Operatormethode erscheint nämlich zur Behandlung linearer 
Differentialgleichungen insofern besonders geeignet, als sie direkt 
auf den Prozessen der Differentiation und Integration aufgebaut 
werden kann. 

Das Fourierintegral hat ja vermöge des Kerns + .ert die 
Eigenschaft, daß die Multiplikation des Integranden mit der nten 
Potenz der Integrationsvariablen p der nmaligen Differentiation des 
Integrals nach der reellen Veränderlichen ¢ gleichkommt. Es gilt 


t.f/( t, „ap. 
Diese Differentialbeziehung, die zunächst nur für ganze positive 


n-Werte besteht, läßt sich bei geeigneter Definition des Operators (+) f 


dt 
für alle rationalen n-Werte aufrechterhalten. Somit ist es auch 
möglich, den Operator durch die Beziehung 


t 

se 


zu definieren. Entwickelt man yp (ae) in geeigneter Weise in eine 


Potenzreihe 
d d \n 
(ae) = (ae) > 
so ist die Anwendung dieses Operators auf obige Definition zuriick- 
geführt. Nach Erlangung einiger Sätze über die Konvergenz von 
Operatorreihen ist man dann imstande, lineare Systeme in ganz 
allgemeiner Weise zu berechnen. So kann man mit Hilfe der 
Operatorrechnung jedes Wellenbild der linearen Kontinuumsdynamik 
aus begrifflichen Einheiten aufbauen, denen gewisse Operatoren 
äquivalent sind, ganz so, wie sich das betreffende System aus 
einzelnen Elementen zusammensetzt. 


Vorbemerkungen 


Es folgt nun ohne Beweis eine Zusammenstellung der für diese Arbeit 
wichtigen mathematischen Tatsachen. Man hat sich gewöhnt, in der Operator- 
rechnung das Fourierintegral mit reellem ? 


a+ioo 
1 t 
(D fe f (p) Fit), 
a-ico 


in dem das positiv reelle a eine obere Schranke fiir die Realteile der singulären 
Stellen p = o, des Integranden ist, durch die „symbolische Form“ 


(La) f(p) > FO 
(—> lies: bedeutet) abzukürzen, und die Auswertung des Integrals als ,,Inter- 
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1 O. Heaviside, a. a. O. 
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pretation oder „Deutung“ zu bezeichnen. Die so definierte Funktion F (t) 
verschwindet, solange 
lim |f(p)|s lim const » |p|", 
rl>o 


identisch für negatives ¢. Daher ist diese Darstellung der Funktion F be- 
sonders geeignet, um die Bewegung mechanischer Systeme aus der Ruhe heraus 
zu berechnen. 
Das Integral (I) hat folgende Eigenschaften " % 3): 
1. Es gilt, sofern die benutzten Limiten existieren, 


(IIa) lim F® = lim f(p) 
» t>+0 
(IIb) lim Fo. 
2 +0 


Haben alle singulären Stellen des Integranden negativen Realteil, so gilt ferner 
(III) lim Fi)= lim f(p). 


t—> co 
2. Wie schon in der Einleitung erwähnt, kann nach Angabe der Ent- 
wicklungsvorschrift 


und Erklärung des Operators (+) fiir rationale n die ns 2s 5 
d 
‘> Gi) = |— 
(V) (p)-f(p) > GO= (ar) Fo 
aufgestellt werden. Wir wollen diesen Rechenschritt „partielle Interpretation“ 
nennen. Für ganze positive n ist der Operator bereits erklärt, für ganze 
negative n bedeutet er, wie sich aus der Betrachtung des Integrals (I) un- 
mittelbar ergibt, die nmalige Integration mit der unteren Grenze — e, wo 8 
eine beliebige positive Zahl bedeutet. Mitn=—k(k=1,2...) gilt also 


d k 4 te 


—e 


(N.B. für negatives ¢ verschwindet F identisch.) 
Die TER von (I) für f(p) = p- ergibt bei rationalem r 


(IT ist die Gausssche J7-Funktion.) Hieraus folgt für die ebenfalls rationalen 
Zahlen n und k die partielle Interpretation 


(va + (4 Ik-m’ 


womit der Operator | 2 y bei rationalem n und Anwendung auf eine Potenz 


erklärt ist. 


| 2 
| aa 
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Bs Die allgemeine Erklärung dieses Operators bei Anwendung auf eine 
durch (I) darstellbare Funktion Ft) wird mittels eines Satzes von Borel’) ge- 


ae 2 leistet, der in symbolischer Schreibweise lautet: 

| > 


OW ud f(p) > Fo. 
Da in folgender Arbeit neben den ganzen nur halbzahlige n benutzt werden, 


geniigt es offenbar, den Operator (ae) ” zu erklären, womit dann gleich- 
zeitig alle andern Operatoren mit halbzahligem n erklärt sind. Setzt man also 
speziell » (p) = p"/:, so ergibt sich unter Benutzung von (VII) 


= 
ax) 2 + f(p) —> f F(it—1-dt. 


Betreffs der Entwicklung (IV) lassen sich folgende Konvergenzsätze un- 
schwer herleiten ?): 

A. Man erhält nur dann eine konvergente Entwicklung, wenn die zur 
gliedweisen partiellen Interpretation benutzte Reihe 


IVa) = > An p” 


für beliebig große |p| konvergiert, denn dann konvergiert die Entwicklung des 
Integrals (I) auf dem ganzen Integrationsweg, sofern nur a groß genug ge- 
wählt wird. 

B. Ist ([Va) eine Entwicklung nach fallenden p-Potenzen, so konvergiert 
die partielle Interpretation, für alle reellen ¢-Werte im Definitionsintervall des 
„Operanden“ F. 

Ist dagegen (IVa) eine Entwicklung nach steigenden p-Potenzen, so 
konvergiert die partielle Interpretation für hinreichend große ¢ > T>0. Die 
endliche Zahl T ist im Einzelnen besonders zu berechnen. Ist der Operand 
N + An | 


für positives ¢ beliebig oft differenzierbar, so ist T = lim | - . 
n—> oo | An_ı | 


C. Konvergiert die benutzte Reihe ([Va) nur in einem endlichen Gebiet, 
so führt die partielle Interpretation auf eine asymptotische Lösung, häufig in 
Form einer semikonvergenten Reihe. 

3. Operatoren von der Form 


d i 
rl>o P 


ist, liefern bei Anwendung auf eine für {> (0) nicht identisch verschwindende 
Funktion Ft) ein Ergebnis, das für {< r identisch verschwindet, wie man 


1) E. Borel, Legons sur les Series Divergentes. S. 104. 1901. 
u 2) K. W. Wagner, Arch. f. Elektrotechnik 4, S. 159ff. 1916. 
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bei Integration von (I) über den unendlich großen Halbkreis, auf dem Re(p) > 0 
ist, unmittelbar erkennt. Der spezielle Operator 
d 


?\ at) ~ 
bedeutet eine Verschiebung der t-Achse um r. Es gilt 


e { 
Der wid a 
. Der erstandsoperator 


0, für t< 1, 


a) Herleitung 
Es möge zwischen der Kraft P(t) und der durch sie bewirkten 
Auslenkung Y (t) ihres Angriffspunktes die lineare Differentialgleichung 
mit positiv reellen Koeffizienten 


d\n d\n 
_ bei endlicher oder unendlicher Summation über n bestehen, und es 


, . mögen die Funktionen P(t) und Y(t) für negatives Argument 
identisch verschwinden, dann geht diese Differentialgleichung durch 


die Laplace-Transformation 


(<— lies: dargestellt durch) in die „symbolische Gleichung‘ : 


oder, was dasselbe ist 


über. In ausführlicher Schweibweise bedeutet dies — vgl. Vor- 


 bemerkungen (la) — daß das Integral 


a+ioo 


identisch verschwinden muß. Dies ist der Fall, wenn zwischen den 
symbolischen Funktionen y(p) und f (p) die Besishung oa 


(3b) y(p) = -f(p) baw. f(ip . 


besteht. 
Führt man nun die Abkürzung 


| 
€: c_ 
| 
— 
| | z(t) = — g(p) > 
ay 


_® 


“4 


= 


er 
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at 
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(6) bzw. f(p) = R(p)- y(p) 


an. Durch Anwendung der „partiellen Interpretation“ erhält man 
hieraus die Operatorbeziehungen 
(7) Y@= Pi) baw. PH =R (4) Y (@). a 
(a) 


Den linearen Differentialoperator R (4) wollen wir mit Heaviside 


den Widerstandsoperator des betreffenden Systems nennen, und in 
ihm ein Symbol für die lineare Differentialgleichung sehen, der die 
Bewegung des Kraftangriffspunktes genügt. Aus der Herleitung er- 
gibt sich die überaus wichtige Möglichkeit, mit linearen Operatoren 
nach den Rechengesetzen der Algebra zu verfahren. Daher 


wollen wir auch, soweit sich dazu Gelegenheit bietet, unendliche 
Summen formal durch geschlossene Ausdrücke wiedergeben, worunter 
man dann immer eine auf dem Integrationsweg des Fourierintegrals — 
vgl. Vorbemerkungen (la) — konvergente Reihe zu verstehen hat. 
Z. B. ist mit dem Ausdruck 

die Reihe 

gemeint. 


b) Die mechanische Widerstandslehre 


Man darf, wie soeben gezeigt wurde, die linearen Operatoren 
wie Zahlen kombinieren. Das bedeutet aber, daß man die Betrachtung 
linearer mechanischer Systeme dadurch erleichtern kann, daß man 
bei der Rechnung zunächst jedes Teilsystem als rein elastische Feder 
behandelt, — denn hier vermittelt ja eine Zahl (die Federkonstante C) 
den Zusammenhang zwischen Kraft und Auslenkung: 


C-Y@=P, — 
und erst am Ende der Rechnung die bisherigen Federkonstanten 
durch die richtigen Widerstandsoperatoren ersetzt. Es ist, mit 
anderen Worten, möglich, eine mechanische Widerstandstheorie zu 
formulieren, die der von Heaviside für die Elektrik aufgestellten 
analog ist, und die die verschiedenen „Schaltungen“, d.h. Kopplungen 


zum Gegenstand hat. 


ag 
2 

= 

= 

= 
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Um nicht unnötig neue Bezeichnungen einzuführen, wollen wir 
unsere Ausdrucksweise weitgehend der in der Elektrik üblichen 
anlehnen. Wir wollen auch hier von den beiden „Polen“ eines 
Widerstands sprechen, und darunter 
— wie dies bei einer Feder ohne P— Rk 
weiteres klar ist — die Angriffs- 
punkte der Kraft P(t) und derGegen- Fig.1. Kraft und Gegenkraft an 
kraft Q(t) = P(t) verstehen (Fig. 1). 

schen Widerstands 
Um eine kurze Bezeichnung zu haben, 
wollen wir den einen den Hauptpol, den anderen den Gegenpol 
nennen. 

Die Angabe eines Widerstands soll sich immer auf den Fall 
beziehen, daB der Gegenpol im Bezugsraum fest (raumfest) ist. 
Es sollen nun die beiden wichtigsten Schaltungen der zweipoligen 
Widerstiinde betrachtet werden. 


1. Die mechanische Parallelschaltung 


Werden n mechanische Widerstände 


deren Gegenpole raumfest sind, parallel geschaltet, d.h. in ihren Haupt- 
polen fest verbunden, so daß ihre Auslenkungen zwangsläufig einander 
gleich sind (vgl. Fig. 2), so benötigt 

man, um die Auslenkung Y(f) zu 

bewirken, für jeden einzelnen Wider- 

stand R, die Kraft 


P,$=R, (7) 


und somit für das ganze System die _ 
Kraft Fig. 2. Parallelschaltung von 


zwei mechanischen Widerständen. 
a d (R bedeutet: raumfest, P—> mar- 
®) kiert den Kraftangriff). Resul- 
nr. tierender Widerstand: R= R, + R, 
Durch Vergleich mit der Definitions- 


Parallelschaltung 


_ dieselbe Formel, wie für die elektrische Serienschaltung. Ersetzt 
man die Widerstände durch Federn, so gilt: 


(9a) R> Max(R). 


4, 
x 
| 
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Ps ; ae Denkt man sich eine Folge von n mechanischen Widerständen 

so miteinander daß der Gegenpol des sten Widerstands 
; a r er mit dem Hauptpol des (i + 1)ten fest gekoppelt, und der Gegenpol 
a des nten raumfest ist (Fig. 3), und wirkt auf den Hauptpol des 
B: ersten Widerstands die Kraft Pt, so unterliegen sämtliche Wider- 
> stände dieser Serie derselben 


seitige Verschiebung seiner Pole, 


R > Kraft P(t. Demnach erleidet 
x jeder Widerstand R, eine gegen- 


Fig. 3. Serienschaltung 


von zwei mechanischen Widerständen. die durch Br 

Resultierender Widerstand: 10 1 Pt 

R, + RB, di 


_ gegeben wird 

Die Verschiebung des ersten Hauptpols (Kraftangriffspunktes) gegen 
| den Bezugsraum ist aber die Summe aller Verschiebungen innerhalb 
der einzelnen Widerstände, nämlich 


n 


(11) 


Um den resultierenden Widerstand dieser Anordnung zu finden, 
0 braucht man — gemäß Gl. (7) — nur 


re zu setzen. Man erhält also für die mechanische Serienschaltung 
dieselbe Formel wie für die elek- 


trische Parallelschaltung. Hier be- 
a steht fiir Federn die Ungleichung 


Fig. 4. Widerstandsskizze eines 
Massenpunktes. (Der Gegenpol ist R<Min(R) 
raumfest!) Man beachte, daß bei Trigheits- 


widerständen (Massen, Trägheitsmomenten) 
in der Regel nur über den Hauptpol verfügt werden kann. Der Gegenpol 
ist als raumfest anzusehen. Ein Massenpunkt muß also in der Schaltskizze 
stets wie in Fig. 4 gezeichnet werden. Eine denkbare praktische Ausführung 
eines Trägheitswiderstands mit zwei verfügbaren Polen ist in Fig. 5 angedeutet. 


“= 
1 
E 
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ee Yen nun an sollen alle Widerstände kurz mit R an Stelle von 


R(=;) bezeichnet werden, denn man darf ja mit linearen Operatoren wie 


mit algebraischen Zahlen verfahren. 


Schwungrad 
Blegsame 
Well 
Pp 


<——~ 
Differentialgetriebe 


Fig. 5. u =: = Fig. 6. Eine kombinierte Parallel- 
Modell einer zweipoligen Trägheit Serienschaltung. Resultierender 
Widerstand: 


e) Anwendung: Die unendliche „Kette“ 

Als Beispiel für eine kombinierte Parallel-Serienschaltung soll 
jetzt der Fall betrachtet werden, daß (Fig. 6) zwei Wiederstände R, - 
und R, in ihren Hauptpolen fest gekoppelt sind, während der Gegen- 
pol von R, raumfest und der von R, an einen weiteren Widerstand R, 
angeschlossen ist. Der Gegenpol von R, sei wiederum raumfest. 
Diese Anordnung läßt sich Br als eine Parallelschaltung des 
Widerstandes R, mit der Serie (R,, R,) ansehen, und man findet so 
unter Anwendung von (9) und (12) als resultierenden ng des 


ganzen Systems den Ausdruck ~~” 

Aus dieser Formel läßt sich leicht der Widerstand einer end- 
losen Kette von mechanischen Systemen berechnen, die (vgl. Fig. 7) 
aus einer unendlichen Serie von gleichen Widerständen R, besteht, 
an deren Polen (in der Figur mit 0,1,2,3,... bezeichnet) je 
ein Widerstand R, abgezweigt ist, dessen Gegenpol raumfest sei, 
Eine solche Kette hat nämlich die Eigenschaft, daß ihr Widerstand 
— falls er endlich ist — durch Hinzufügen eines weiteren „Gliedes“ 
nicht geändert wird. Angenommen nun, dieser Widerstand sei R,, 
so kann man noch ein weiteres Kettenglied hinzufügen, indem man 
den Widerstand R, davorschaltet (in Serie) und an den neuen Haupt- 
pol den Widerstand R, ankoppelt, dessen Gegenpol raumfest ist. 
_ R, liegt demnach der ganzen Kette parallel, und man kann obige 
Formel (13) anwenden, um den Widerstand der um ein Glied ver- 


| £ | 
—— 
\ 


58 Annalen der Physik. 5. Folge. Band 24. 1935 


längerten Kette durch R,, R,, R, ausdrücken. Berücksichtigt man 
nun, daß dieser gleich dem Widerstand R, der ursprünglichen Kette 
sein muß, 50 a. man nach Einsetzen von 


Fig. 7. Schaltskizze der ersten drei Glieder einer „Kette“. 
Resultierender Widerstand der „Kette“: 


in (13) die quadratische Gleichung fiir R ia 


(15) R= + 
deren Lösung unter der für Parallelschaltung gültigen Voraussetzung 
— vgl. (9a) — 

R>R, 


den Widerstand der Kette 
(16) 


liefert. 
Die Bewegung des n-ten Kettengliedes läßt sich ebenfalls aus 
der Widerstandslehre schnell ermitteln. Einer im Punkte 0 der Fig.7 
angreifenden Kraft P,, die dort die Auslenkung Y, hervorruft, wird 
ja das Gleichgewicht gehalten durch die Widerstandskräfte von Ry 
und der mittels R, angekoppelten Kette. Diese sind u 


P=RY, ud P,= Y, 


Ihre Summe ist gleich P,. Es wird also durch den Win B, 


auf den Punkt 1 der Kette die Kraft P, 
P,=P,— RY, 
übertragen, die dort an dem Widerstand R die Auslenkung 


Yo P, P,- Rk, Y% 
a R R 


Setzt man nun für P, den Wert RY, ein, so ergibt sich 


erzeugt. 
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zwischen den Auslenkungen am Ende und am Anfang des ersten 
Kettengliedes, Y, und Y,, die Operatorgleichung 


(17) Y, = (1- 

Dieselbe Beziehung besteht zwischen den Auslenkungen am Ende 

baw. Anfang jedes anderen Kettengliedes. Somit berechnet sich die 
_ Auslenkung Y, des Punktes n der Fig. 7, der durch n-Glieder vom 


Punkte 0 getrennt ist, aus der Auslenkung Y, durch n-malige An- 1 
wendung des — 


Diese Beziehung besteht zwischen den Auslenkungen aller tint n- 


Kettenglieder getrennten Pole. Es gilt mit positivem r 


d) Der Widerstand eines Kontinuums 
Im Anschluß an die im vorigen Abschnitt berechnete Kette 
von mechanischen Widerständen kann man leicht zu der Betrachtung 
eines eindimensionalen Kontinuums übergehen. Zu diesem Zweck 
muß man jeden der Widerstände R, gleichmäßig in eine Serie von 


m-Widerständen R® mit dem resultierenden Widerstand R, unter- 


teilen, deren Größe sich zufolge von Gl. (12) auf 


_ Fig. 8. Die ersten zwei Glieder der „Kette“ (Fig. 7) in m-facher Verfeinerung 


berechnet (vgl. hierzu Fig. 8). Gleichzeitig muß man die Wider- 
stände R, auf die alten und die hierdurch neuentstehenden Pole [in 
der Figur mit (0,0), (0,1), ... (0, m)=(1,0), (1,1)...(1,m) usw. bezeichnet] 
gleichmäßig verteilen, so daß jetzt — vgl. Gl. pr — an jedem Pol 


m 1 ; 
(20) Ry = 


Dar 
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angekoppelt ist. Den Widerstand R™ dieser m-fach verfeinerten 
endlosen Kette berechnet man nach (16), wenn man dort wd 


bzw. RS” für R, baw. R, einsetzt, zu 


(m)_ 1 m 1 ptm m) lm 
oder, nach Einsetzen der Werte aus (19) und (20) zu . De Fr 


(22) 


1 


Bei unendlicher Verfeinerung der Teilung, d.h, wenn m gegen un- 
endlich geht, ergibt sich hieraus der Widerstand W eines nach einer 
Seite hin unendlich erstreckten, eindimensionalen homogenen Konti- 
nuums, bei dem auf irgendeiner Teillänge — ein Maßstab wurde 
ja nicht benutzt — die Gesamtheit der Widerstände R/”, die von 
der Feinheit der Teilung unabhängig ist, den Wert R, ausmacht, 
während gleichzeitig die Serie der Widerstände Ry” (ebenfalls von 
der Verfeinerung unabhiingig) den resultierenden Widerstand R, hat. 
Man findet 
(23) W =lim R™ = yR, - R, . 

Wir wollen fortan ein solches Kontinuum einen ,,Stab“ nennen und 
dadurch ein Längenmaß einführen, daß wir alle Widerstände auf 
eine Teillänge von 1 cm beziehen. Den resultierenden Widerstand 
einer 1 cm langen Serie der Widerstände _ nennen wir dann seine 
„Materialeigenschaft‘‘ und bezeichnen ihn mit a,. Ferner wollen wir 
die Gesamtheit der auf einem 1 cm langen Stück des Stabes liegen- 
den Widerstände R)" seine „Raumbindung‘‘ nennen und sie mit a, 
bezeichnen. Mit den neuen Bezeichnungen stellt sich der Wider- 
stand eines unendlichlangen Stabes als geometrisches Mittel von 
Raumbindung und Materialeigenschaft dar, 
en W = 

Als Raumbindung muß man stets den Trägheitswiderstand pro Längen- 
einheit in Rechnung stellen. Bedeutet y die spezifische Masse, so ist in diesem 
Fall bei einem Stabquerschnitt F 


d d \? 
Ay (+7) = YF- (47) 


Darüber hinaus käme etwa noch eine Abstrahlung X "7 (Raum- 


dämpfung proportional der Geschwindigkeit) und eine elastische Bettung e in 
Betracht. Die Konstanten K und e sind so zu verstehen, daß die Ausdrücke 


zu setzen. 


[EZ 
= 


G, (a, t) = e- 2) 


und 


| 
AB 
. 
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die Kräfte pro Längeneinheit darstellen, welche durch diese Raumbindung 
und die Bewegung Y««,t) an der Stelle x entstehen, so daß im allgemeinen 
Fall als Raumbindung der Operator 


Ir F d \* K d 
zu setzen ist. 


Die einfachste Materialeigenschaft ist die Elastizität. Hier ist, wenn man 
mit E den Elastizitätsmodul und wieder mit F den Stabquerschnitt bezeichnet, 


=EF. 


Allgemeinere Materialeigenschaften werden z. B. von Jeffreys’) betrachtet. 
Aus der Defininition des Operators a, folgt, daB stets zwischen der Deh- 
nung &(£) und der Spannung o (f) die lineare Differentialgleichung 


(26) 


d 


B. Der Wellenoperator 


a) Herleitung des Operators aus der endlosen Kette 

Es soll nun im Anschluß an die Operatorbeziehung (18) die 
Auslenkung Y(z,t) an der Stelle x eines Stabes durch den „Wellen- 
operator* auf die Auslenkung Y(0,t)= Y,(t) des Anfangspunktes 
bezogen werden. Zu diesem Zweck muß man von einer solchen 
Verfeinerung der oben berechneten Kette ausgehen, daß sowohl der 
Punkt x (z sei rational) als auch die Punkte mit den Koordinaten 
1,2,3... auf Pole der Kette fallen. Angenommen auf die Längen- 
einheit entfallen r Kettenglieder, so ist der Punkt x durch z-r 
Glieder vom Anfang der Kette getrennt, und seine Auslenkung ergibt 
sich aus Gl. (18), wenn man dort x-r an Stelle von n und — vgl. (20) — 


= +a, an Stelle von R, einsetzt. Für R hat man, entsprechend 
Gl. (22) den Wert 
1 i 
einzusetzen. Geht man dann zu unendlicher Verfeinerung der Kette 


iiber, indem man r gegen unendlich streben läßt, so erhält man die 


Oh) 


1 1 2 


1) Jeffreys, Proc. Royal Soc. London, A 138. S. 283. 1932; auch 


 K.Hohenemser u. W. Prager, ZAMM 12. S. 216 ff. 1932. 
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ergibt sich 


(28) Y (a, t) = lim (1 + (t). 


Der liefert bekanntlich das Ergebnis’) 


d 

Den Operator e =: (2) wollen wir den „Wellenoperator‘‘ nennen. 
Er beherrscht die Wellenausbreitung im Stab für 2>0- Für seine 
Gültigkeit braucht nur vorausgesetzt zu werden, daß der Stab für 
xz > 0 homogen ist, daß er also insbesondere keine Einzelwiderstände 
trägt. Es ist aber nicht notwendig, daß die Stelle x = 0 gerade 
am Anfang des Stabes liegt, sie darf auch irgendwo im Innern 
liegen, denn gleichzeitig mit der benutzten Beziehung (18) gilt ja 
auch (18a), so daß man ein beliebiges Kettenglied als Ausgang der 
Betrachtung wählen kann. 

Man kann den Wellenoperator natürlich auch auf direktem 
Wege aus der allgemeinen homogenen Differentialgleichung für Stab- 
längswellen 


(30) — 4, Y” (a, t) + a, (4) Y(z,t) = 0 


(” bedeutet die zweimalige Ableitung nach x) gewinnen?). Hier ist 
jedoch die formale Einführung der Randbedingungen für x xr =0O 


nicht sehr befriedigend. _ 
7 


b) Diskussion des Wellenoperators . 
Fiir die Diskussion des Wellenoperators ist es von Vorteil, 
die Zeitfunktionen durch ihre Symbole 
(31) y(t,p) > Y(a,f) und y,(p) +> Y,( 
zu ersetzen. Dann hat man gleichzeitig q(p) an Stelle von a(42) 
zu schreiben. 
> 1. Die Front des Wellenzuges 


Setzt man voraus, daß lim 4” ap) existiere, so kann man den 
|p 


Wellenoperator mit dem in den Vorbemerkungen (X) erwähnten 


1) Z. B. R. Courant, Vorlesungen über Differential- und Integralrechnung, 
Bd. 1. Berlin 1927. S. 32. 
z te 2) W. Prager, Ing. Arch. 4. S. 16. 1933. 
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Verspätungsoperator identifizieren, und den Exponenten q in die 
drei Teile zerlegen 


(32) 
wo Fuck 

8) 7 im 2) 
und 


ist. Man erhält so für 2>0 die Auslenkung Y (x, t) in der symboli- 


schen Form 


und abe, daß c die Frontgeschwindigkeit der Welle bedeutet, 
denn die Bewegung beginnt — vgl. Vorbem. (Xa) — an der Stelle x 


um die Zeit = = später als an der Stelle 0. Ferner nimmt die 


Frontintensität nach dem Gesetz e~** ab. Es ist nämlich, wenn 


2.2 
man den Operator e ¢ bereits Tn 


für t>—, und die Welle setzt an der Stelle x zur Zeit t = — 


mit der Auslenkung 
(35) Y (2, =) (az) } ( 


ein. Der Wert für ¢ kann erst nach Ausübung des Verzerrungs- 
operators e= **° substituiert werden. Man erhält so 


wo die symbolische Darstellung von F(é) durch 


(86a) F(t) 


‚gegeben ist. Läßt man zum Zweck ui in (86) Ben Sub- 
stitution den positiven Ausdruck 


7 
-=) 
= 
| 
2 


> 


wah 2. Das (stationäre) Innere des Wellenzuges 


Um das Innere des Wellenzuges bei einem Stab zu berechnen, 
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gegen Null gehen, und berücksichtigt die in den Vorbemerkungen 
zitierte Beziehung (Ila), so erhält man zunächst 


(37) lim F(e)= lim 

und mit (33b) sowie der Umkehrung von (IIa) wi 
(38) lim F()= lim y,(p) = Y, (0). 

|p| 
So ergibt sich schließlich die Fronthöhe aus (36) und (38) zu 
(39) h(x) = (2, = = erkı, Y, (0) = e- kz, h (0) 


Ist der benutzte lim y,(p) gleich Null, so hat man diese Be- 
|p|—>oo 


rechnung nicht fiir die Auslenkung Y, sondern etwa fiir die Spannung 
des Materials a, Y’ anzustellen, um auf den Faktor e- ** zukommen. 


4 


muß man im Wellenoperator d/dt durch iN ersetzen, wo N die 
Kreisfrequenz der Welle ist. Man kommt darauf, wenn man 


(40) Y,() = A(N).eixt 


setzt und den Wellenoperator für große Zeiten betrachtet, indem 
man ihn nach steigenden Potenzen von d/dt entwickelt. Dies 
Verfahren führt nach Satz C der Vorbemerkungen zu einer asympto- 
tischen Lösung, da die Operatorreihe in der Regel nicht auf dem 
ganzen Integrationsweg des Fourierintegrals konvergiert. 

Die so erhaltene Lösung gilt im allgemeinen nur in hinreichender 
Entfernung von der Wellenfront. In der Nähe der Front hat man 
den Verzerrungsoperator e~*-* in eine auf dem Integrationsweg 
überall konvergente Reihe zu entwickeln. 

Bezeichnet nun g(N) den Imaginärteil von q(iN), so erhält 
man für den Wellenanteil mit der Frequenz N den Ausdruck 

(y) 
(41) Y (a, ) 
wo A(N) seine Amplitude bezeichnet. Hier erkennt man, daß die 
Wellengeschwindigkeit durch 


(42) o(N) — 


gegeben ist, woraus sich für die Gruppengeschwindigkeit 


: 
a 2 a 
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C. Der Reflexions- und der Deformationsoperator 
a) Einfluß eines angekoppelten Systems auf das Wellenbild 

Es soll nun der Einfluß einer beliebigen Randbedingung auf 

das Wellenbild betrachtet werden. Es möge in einem Stab von 

links her eine Welle anlaufen, die an einer Stelle s=0 irgendwie 

_ erzeugt wurde. Wir wollen sie die „Primärwelle“ nennen und ihre 

 Laufrichtung durch einen nach rechts weisenden Pfeil andeuten. 

Diese Welle wird also nach (29) für 2>0 und t>0 durch die 
Operatorgleichung 


—> 
dargestellt. Ist nun an den Stab im Punkte x =a>O0 ein mecha- 


nisches System mit dem Widerstand RS) angekoppelt, so kann 


man, da alle Differentialgleichungen linear sind, dessen Einfluß auf 

das Wellenbild dadurch berücksichtigen, daß man die durch die 

Reaktionskraft des Systems entstehenden sekundären Wellen der 
Primärwelle überlagert. 

Bezeichnet man mit Y(a,t) die an der Stelle 2 =a wirklich 

auftretende Bewegung, so wird diese Reaktionskraft Q(f) durch den 

Ausdruck — vgl. (7) — 


(45) 


x: | gegeben. Sie bewirkt, da der Widerstand des nach beiden Seiten 
hin unendlichen Stabes nach Gl. (24) 2W beträgt, an der Stelle x =a 
die sekundäre Auslenkung 
1 
(ae) 


(46) Yı(a, 

so daß sich die nach rechts >a) fortschreitende 
in der Form 

(47a) 


QO, 


t) 


und die nach links (« < a) fortschreitende durch den Ausdruck * 


darstellen läßt. Demnach beträgt die nn — Y (az, t) im 
Intervall OS 2Sa 
(48a) = Yı + Yır 
und fir 
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Speziell gilt für die Auslenkung Y (a, t) der Kopplungsstelle x =a 
49) ¥ (a, t) = Yr(a, ) + Yu, 8), 

a. woraus sich die bis jetzt noch unbekannte Funktion Y (a, t) be- 
rechnen läßt. Durch Einsetzen von (45) in (46) findet man 


-R, 
50 Yn (a, t) = Yla,t 
und unter Verwendung von (49) folgt a nee 


-R 
Yırla, 


(Yı (a, t)+ Yır(a, t)) 


Die Auflösung dieser Operatorgleichung nach der unbekannten 
Funktion Yy (a, t) ergibt 


r -R 

(51) Yun(@, = (a, t). 

Hiermit sind alle Funktionen bekannt, und wir wenden uns der 
Betrachtung des Wellenbildes zu. 


1. Die reflektierte Welle 

Wir wollen fortan die der Primärwelle entgegenlaufende Sekundär - 
welle als reflektierte Welle bezeichnen und ihr den Index R geben. 

Unter Verwendung von (47b), (51) und (44) läßt sie sich in die Form 


(52) 
bringen. Den Operator 


(53) "= 


nennen wir den Reflexionsoperator, denn er charakterisiert die bei 
der Reflexion entstandene wesentliche Veränderung der Welle. Für 


seinen Einfluß auf die Wellenfront ist das Verhalten im Fre r (p) 
|->0o 
maßgebend, wie man im Anschluß an die in Gl. (35 ff.) vorgenommene 


Berechnung der Wellenfront unmittelbar einsieht. Und er geht für 
die stationäre Welle der Frequenz N in den komplexen Faktor (Vertor) 
r (N) über, der die Intensitäts- und Phasenverhältnisse bestimmt. 
Man sieht sofort, daß beim ungestörten homogenen Stab, wo R = 0 ist, 
keine Reflexion stattfindet, während ein festes Lager (R = oo) die 
Welle mit entgegengesetzter Phase total reflektiert, denn hier 
ist r=—1. 

Bricht der Stab an der Kopplungsstelle ab, so ist fiir R der 
Wert = R—W einzusetzen, da der fehlende halbunendliche 
Stabteil den Widerstand W hat. Man erhält in diesem all- 
gemeineren Fall den Reflexionsoperator 


„__W-R 

4 = | 


A 
= 
= 
be 
BA, 
Al 
JE 
- 
ae 
> 
> 
4 
am 
| 
> 
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der im Falle eines freien Endes (R = 0) in die Zahl 1 übergeht, 
entsprechend der hier erfolgenden Totalreflexion mit gleicher Phase. 


2. Der Widerstand des gekoppelten Systems 


Im Intervall 0 = «=a wird die wahre Auslenkung Y (z, t) auf = 


die primäre der Stelle z=0 durch die Summe DE 


4 (Yo Yı@,y+ Ya!) 


bezogen, wie sich durch Einsetzen von (44) und (52) in (48a) ergibt. 
Insbesondere besteht an der Stelle x = 0 zwischen der wahren Aus- 
lenkung Y (0, t) und der primären Y; (0, t) die Beziehung 


d 
— womit übrigens G1. (55) in 


eine Beziehung zwischen den wahren Auslenkungen, übergeht. — 
Bringt man also im Punkte 0 eine Einzelkraft P(t) an, so daß 
sich die — Auslenkung auf 


(0,8) = 


57) (a, 1) = 


l+e 


w( az) 


beläuft, so erfolgt nach (56) an dieser Stelle die wahre iailiiiniaita, 


2a-d,» 


l1+e” 
Y (0, t) = 


P(t). 


dt 
9 
Er enthält den Verspätungsoperator e—?«?, der für die Zeit r< a 


_ gleich Null zu setzen ist, und man sieht, daß sich der Einfluß der 
_ Randbedingung bei x =a erst nach Verlauf der Zeit r im Kraft- 
 angriffspunkt « = 0 geltend macht. Wird a = 0 gesetzt, so ergibt 
sich für das System in Übereinstimmung mit der für Parallel- 
schaltung gültigen Formel (9) der Widerstand 2W + R. TR, 


. 
» 
v . 
s 
hieraus den Widerstand des 
| 
: 
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3. Die durchgehende Welle in 


Nach Gl. (48b), (47a), (51) und (44) findet man für den durch- 
gehenden Wellenzug, dem wir den Index D geben wollen, den 
Ausdruck 


Y,+Yn für >a 
(60) 


= [e-=7 für 


| =e--:.DY,(0,) für z>a. 
Der Operator 


ew 


ist für die beim Passieren der Stelle «=a erfolgte Deformation 
der Welle maßgebend. Im Gegensatz zu dem bei einem Kontinuum 
auftretenden Verzerrungsoperator e-*'° hat der Deformations- 
operator vor allem zwei Eigenschaften: 

1. Er hängt in unstetiger Weise von der Ortskoordinate x ab 
und kann in bestimmten Grenzfällen (wie bei der Herleitung des 
Wellenoperators aus der „Kette“) gegen einen Wellenoperator kon- 
vergieren. 

2. Nennen wir n den Grad der Unstetigkeit der Wellenfront, 
wenn sich die symbolische Funktion y(0, p) bei unbegrenzt wachsen- 
dem p wie c, - p” verhält, und m den Grad der Deformation, wenn 
sich mit wachsendem p die Funktion D(p) der Grenzfunktion c, -p=" 
nähert, so sehen wir, daß der Deformationsoperator den Grad der 
Unstetigkeit an der Wellenfront herabsetzt oder ihn unverändert 
läßt. Denn y(0,p) geht für >a in die symbolische Funktion 


(62) 7(0, p) = D(p) - y(0, p) 
über, die sich im lim 7 (0,p) wie c,+c,-p"-” verhält, und da 
|p| eo 


— vgl.(1) — die Widerstände W und R nur positiv reelle Koeffi- 
zienten haben, ist m nicht negativ. 


b) Einfluß mehrerer gekoppelter Systeme auf das Wellenbild 
1. Der beiderseits begrenzte Stab 


Es ist nun leicht, das Wellenbild bei einem angekoppelten 
System auf das bei mehreren Ankopplungen entstehende zu er- 
weitern. Wir betrachten zunächst einen Stab von der Länge a +b, 
der an seinem linken Ende von einem System mit dem Wider- 


| 
> al) 
| 
f - 
S 
| 
1.2: 
a 
rs 
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stand R, an seinem rechten von einem solchen mit dem Wider- 
stand R, begrenzt ist, und im Abstand b vom linken Ende, a vom 
rechten der Kraft P(t) unterworfen wird. Wir wollen die Be- 
wegung der Stabpunkte auf zwei voneinander unabhängigen Wegen 
tinden. 
Durch wiederholte Reflexion der vom Kraftangriffspunkt z = 0 
nach beiden Seiten ausgehenden Primärwelle an den Enden erhalten 
wir das Wellenbild in Form einer unendlichen Summe der reflek- 
_ tierten Wellenzüge aller Ordnungen. Dieser Weg empfiehlt sich 
besonders, wenn man die Lösung für relativ kleine Zeiten sucht. 
Es genügt, die Rechnung für den rechten Stababschnitt 0=2=a 
durchzuführen. 


Bezeichnet man mit: 
Yra die bei a einmal reflektierte Welle, 


Pau die bei a und b je einmal reflektierte W elle, h 


| 

al EUER die bei a zweimal und bei b einmal reflektierte 

Welle usw., 
7. so hat man in dem Intervall 0=2=a folgende Wellen zu ad- 
_ dieren, um die wahre Auslenkung zu erhalten: . 
< > < lig 
Y = Y,+ Yrat Yraot Yravat 
+ Yrot+ Yroat 

Nach Umordnung der Summanden und Abspalten der gemeinsamen 
_ Wellenoperatoren ergibt sich unter wiederholter Anwendung von (52): 

+ r+ e-2(a +b)g. roe t+ -+-) Y1(0,8). 


fetta. + F*+ ---)¥1(0,t) 
— 
wobei 
(66) F = e-?la+b)g, ret, 
Er wurde, md unter r, bzw. r, der mit dem Endwiderstand R, 
zw. R, nach (54) gebildete Reflenionespereter zu verstehen ist. 
Vorsmagenntak, daß sich, wie in speziellen Fällen leicht möglich, 
für das Fourierintegral [Vorbem. (I)] ein Integrationsweg finden läßt, 
auf dem überall der Betrag der Funktion F(p) kleiner als 1 ist, 


er 
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kann man die geometrische Reihe (65) summieren und erhält 
schließlich: 


7 e @+y,) 
| = (0,8) 
(67) | — e+22Y7(0,1) | 


Ari 


Hier wird die Bewegung für beliebige Randbedingungen durch 
die Primärwelle und die Reflexionsoperatoren ausgedrückt. 

Um den Widerstand dieses Systems zu finden, hat man für 
7 7 
Yı(0,t) den Wert ee 


1 
(0, = OW 


einzusetzen und die wahre Auslenkung des Kraftangriffspunktes 


(1 + (1 + e~ ,,) 1 
mit der Definitionsgleichung (7) zu vergleichen. Es ist also dieser 
Widerstand: 


-2a+b)g, 


(1 + 948. (1 + ony) 
Wir werden dasselbe Ergebnis jetzt unabhängig von der Kon- 
vergenz der Reihe (65) herleiten. 
Durch eine Widerstandsrechnung wird man unmittelbar auf 
Gl. (69) geführt, wenn man den beiderseits begrenzten Stab als 
Parallelschaltung zweier Stabstücke ansieht, an deren freiem Ende 
die äußere Kraft P(é) angreift, und die am anderen Ende an die 
Systeme mit den Widerständen R, bzw. R, gekoppelt sind. Der 
Widerstand eines Stabes, der nach einer Seite hin unendlich lang 
ist, nach der anderen die Länge a besitzt und an diesem Ende mit 
dem Widerstand R, gekoppelt ist, beträgt nach G1.(59): oe 
7 


Eh Hieraus findet man sofort den Widerstand eines Stabes von 
Bir: der Länge a mit dem Endwiderstand R,, indem man den Wider- 
stand W des halbunendlichen Stabes subtrahiert. Er ist: 


_ = 


: | 
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Sei dieses Stabstück der rechte Teil des soeben betrachteten 
Stabes und sei der linke Teil ein anderes Stabstück von der Länge b, 
mit dem Reflexionsoperator r,, das also. den Widerstand 


hat, so erhält man beim Zusammenschalten der beiden Stabteile nach 


Formel (9) den Widerstand des beiderseits begrenzten Be: ion 
—2(a+b)q, 

R = R'+ R’= 


übereinstimmend mit Gl. (69. Kennt man aber die Auslenkung 


so kann man nach (57) die Auslenkung Y (z,t) jeder Stelle x im 
Intervall 0=2=a angeben. Hiermit ist auf jeden Fall die Bich- 
tigkeit der Formeln (67) und (69) bewiesen. 


2. Der in regelmäßigen Abständen gekoppelte Stab 
Wir wollen nun einen Stab betrachten, der in regelmäßigen 
Abständen an unendlich viele einander gleiche Systeme gekoppelt 
ist. Da die Summation der reflektierten und durchgehenden Wellen 
aller Ordnungen recht unübersichtlich ist, soll diese Berechnung nur 
auf der Widerstandslehre aufgebaut werden. 
5 Den Widerstand R eines einseitig unendlichen Stabes, der im 
Kraftangriffspunkt z=06 und an den Stellen 
an je ein System mit dem Widerstand R, gekoppelt ist, finden wir 
nach derselben Überlegung, die bei der Berechnung der unendlichen 
Kette auf die quadratische Gl. (15) führte. Dieser Widerstand ist 
die Summe des Widerstands BR und des Widerstandes eines Stab- 
eek ie 'stückes von der Länge L, an dessen freiem Ende die Kraft an- 
= greift, und dessen anderes Ende mit dem noch unbekannten Wider- 
stand R gekoppelt ist. Nach (9) und (70) erhält man also für R 
Gleichung: 


: (73) 


zu ist. 


= 4 
| 
| 
: 
iy? 
wi. 
rat 
r der Ausdruck =) 
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Diese Gleichung ist in R quadratisch und ergibt unter der fiir 
Parallelschaltung gültigen Bedingung R > R, die Lösung NE 
Die Bewegung Y (x,t) an einer beliebigen Stelle des Stabes kann 
man vermöge Gl. (57) auf die des unmittelbar vorhergehenden n-ten 
Kopplungspunktes 2 = z,= n-L beziehen und erhält für das Inter- 
vall <r=xz,_, den Ausdruck: | 


(75) Y(z,t) = "Y(q,,t), 
wo nach (73) und (74) De 
R, +/+ .B+W + We Gig Le 
Ward, 


zu setzen ist. Speziell besteht für zwei sukzessive Kopplungs- 
punkte z, und 2,,, die Operatorbeziehung: 


Y = Y t) 


Durch wiederholte Anwendung dieser Formel kann man die 
Bewegung Y (z,,t) des n-ten Kopplungspunktes zu der des nullten 


(des Kraftangriffspunktes) in die Beziehung setzen: in 7 
so daß man schließlich für die Bewegung einer beliebigen Stelle x | 


zwischen x" und 2,4, die Formel 


| ( ? ) 


+e r } (0, t) 


(78) 


gewinnt. Man iiberzeugt sich leicht, daB bei kontinuierlicher Ver- 
teilung der Widerstände R, auf den ganzen Stab der Wellenoperator 


| 
— 
4 
(76) | 
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Hiermit soll der theoretische Teil dieser Arbeit beschlossen 
werden. Die Reihe der Betrachtungen ließe sich noch nach ver- 
schiedenen Richtungen hin erweitern. Abgesehen davon, daß natür- 
lich Torsions- und Saitenwellen demselben Formalismus genügen, 
da sie ebenfalls durch eine Differentialgleichung zweiter Ordnung 

beschrieben werden, kann man auch die Transversalwellen eines 
 Balkens (Differentialgleichung vierter Ordnung) mit Hilfe der hier 
entwickelten Begriffe zusammensetzen. Hier ergibt sich nur da- 
durch eine größere Mannigfaltigkeit, daß man zwei Wellen (Druck- 
und Schubwelle) erhält, die miteinander gekoppelt sind. Auch 
kommen ja hier als Randbedingungen gleichzeitig Kräfte und 
Momente in Betracht. 

Während man bei all diesen qualitativen Betrachtungen ohne 
direkte Interpretation der symbolischen Funktionen auskommt, ist 
eine quantitative Berechnung natürlich nicht ohne Interpretation 
möglich. 

D. Beispiele 
a) Berechnung spezieller Balkenwellen 


Es sollen die in einem beiderseits unendlich ausgedehnten, 


-, anfangs ruhenden Balken von der Störungsstelle ausgehenden Wellen 


gefunden werden, wenn an dieser Stelle (x = 0) 
. eine plötzliche Verrückung 

. ein Kraftstoß u 
P,® =1-s,(, 

. ein Kraftsprung 

Pj=P-s() 

vorgegeben wird. aa, 

Die Funktion s, bedeutet den Einheitssprung 


0, wenn t< 0 


und die Funktion s, hat die Eigenschaft 
| CO, ” t = 0 
u so, daß für beliebige positive Zahlen ¢, und ¢, das Integral 
's 
fs ,@)-dt=1 
-h 


ist dann der durch den Stoß übertragene Impuls. 
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Die Wellenausbreitung mége von der symbolischen Differential- 
gleichung 


IV 

(79) 
beherrscht werden. 

Lösung. Die Wellenlösung der Differentialgleichung (79) lautet 


> -z.4 
81) 4nt= 


und den Integrationskonstanten 
—> Y (0, t), y, —>Y' (0,0). 


Da kein äußeres Moment vorhanden ist, muß die Lösung milde 
sein, d. h. 


So erhält schließlich nach Einsetzen von (81) die Wellenlösung 
die Form 


Den Widerstand des Balkens erhält man aus der üblichen 
Randbedingung 


die im vorliegenden Fall =. 


2 EJ (9, + %)* 


lautet. Setzt man fiir q die Werte aus ae ein, so findet man 
den Widerstand 


(83) Re = 8EJn*p Yp. 


v 


Nun sollen die entstehenden Wellen für die drei vorher genannten 
typischen Störungsursachen berechnet werden. 


1. Die plötzliche Verrückung — 
Die symbolische Darstellung der Funktion 
Y (0,t) = A-s, (t) 


lautet, wie sich aus dem are Fourierintegral (I) ohne 


weiteres ergibt, 4 


>= 
4 
| 
3 
: 
fi 
. 
k 
ae 


i. Y (a,f) = A-(1—28(9), 


\ . konvergiert. Fürt > 0 liefern die Glieder mit ganzem pos. Exponenten 
keinen Beitrag, was man dem Integral (J) leicht ansehen kann’), 
und man erhält bei gliedweiser Deutung die Ausdrücke 


H. Baumann. Operatoren in der Kontinuumsdynamik 75 


Damit liegt sofort die symbolische Lösung der — vor: 


(85) yo -(cosan Yp + sinznyp). 
Hieraus folgt nach (IIb) für t = 0 
0, wenn «+0, 


1. lim ¥(2,t) = + lim y=) 
2 +0 2 5.4, » 
Die Lösung für {> 0 erhält man durch Reihenentwicklung von (85 
nach steigenden Man findet ') 


(86) Y(x,0) = 


e-nzVr. cosnzyp=1-7; “nz Vp + Vp — Zn + - 


und 
e-"zV» .sinnz Vp = _ P+ m py p 
+- | 


Die Reihen konvergieren für alle p-Werte und zwar so, daß nach 
Satz B die gliedweise Interpretation für alle positiven t-Werte 


2! na \*is 22 1.3.n° na \*ls 
= 
Vat 3! 2 \yat +5 2° Vat 


(55)"- 
Vat 2 5! 2? 
Durch findet man so die Interpretation t>0 


- 


Vat 


z 


(87a) S(2)= a*®-da 


1) Hayashi, Fiinfstellige Tafeln, Berlin und Leipzig 1928. S. 181. 
2) L. Casper, Archiv für Elektrotechnik 16. S. 368ff. 1926. 
3) E. Jahnke u. F. Emde, Funktionentafeln, Leipzig und Berlin 1933, 
S. 109. 
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ee. Es läßt sich auch sofort die Kraft angeben, welche erforderlich 
Be ist, um einen Balkenpunkt plötzlich um den Betrag A zu verrücken 
und dann festzuhalten. Ihr Symbol lautet 


Für {= 0 findet man bei Auswertung des Fourierintegrals 
| | = 
P (0) = Jp Vp 
100 


eine positiv gerichtete unendlich groBe Kraft. Der hierdurch iiber- 
tragene Impuls, symbolisch dargestellt durch 8 EJn*A Yp, ist eben- 
falls unendlich groß, jedoch so, daß das Integral der Impulsfunktion 
stets null ist. 

Die Interpretation für positives ¢ lautet — vgl. Vorbem. (VII) — 


4EJn°’A 1 
t) = — —_____ - 


2. Der Kraftstoß 
Die symbolische Darstellung des KraftstoBes ist in dem Ausdruck 


gegeben, denn es gilt allgemein tien der Kraft P (t) und dem 
Impuls J(t) die Beziehung 


zu setzen. Die symbolische Lösung der durch den Kraftstoß « er- 
zeugten Balkenwelle lautet somit 


= —nz I 
(90) y=e Vp. - 


Die Interpretation für > 0 erfolgt, indem man dieses Symbol in 
die Gestalt bringt 


I 


91 — —_ —anV d 2 
(91) y SEIN Vp fe P.sinanVp-da 


Vv 


und den EEE mit der zweiten der oben vorgerechneten 
Reihen identifiziert. So erhält man die Lösung 


6 
a 
= 
a 
E 
— 
k 
d 
Pi)=—_J(, 
und hier ist speziell 
= 
4 
3 
> 
iti 


<a H. Baumann. Operatoren in der Kontinuumsdynamik 77 


Die Funktion S ist durch (87a) erklärt, entsprechend ist 
z 


C(2) = [cos tat. de. 


ein, so erhält Y die Gestalt 


> t 1 I Ynt 
(93) y (x, t) = . C, (2) + 8, ) = F, (2). : 


Die Funktion 


(94) = (cos sin + 2z-S(2) — (2) & 


ist in Fig. 9 aufgezeichnet und ihre Werte sind in folgender Tabelle 
zusammengestellt. Man erhält den zeitlichen Verlauf des Wellen- 


. 
Führt man nun zur Abkürzung Py 
R ¢ ~ 
= z-S(z)+ —-(cos>2?—1 = 
7 
E 
var 
q's) 
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bildes, wenn man die Figur sowohl in der z-Achse als auch in der 
F,-Achse nach dem Gesetz Yt auseinanderzieht, so daß das Bild 
sich selbst ähnlich bleibt. 


z F, (2) | z F, () £ F, (z) 
0,0 + 1,0000 2,0 0,0896 4,0 0,021 
01 0,9844 2,1 0,0338 4,1 — 0,011 
02 0,9379 2,2 — 0,0309 4,2 — 0,025 
0,8 0,8620 2,3 — 0,0744 4,3 — 0,001 
0,4 0,7597 2,4 — 0,0738 4,4 + 0,023 
05 0,6349 2,5 — 0,0358 4,5 0,012 
0,6 0,4931 2,6 + 0,0180 4,6 — 0,017 
0,7 0,3418 2,7 0,0534 4,7 — 0,016 
08 | 0,1896 2,8 0,0514 4,8 + 0,015 
0,9 | 0,0482 2,9 0,0135 4,9 0,014 
10 — 0,0732 3,0 — 0,0311 5,0 — 0,012 
came | | Zope 
13 - 0.2210 33 + 0.0207 
14 — 0,1858 3,4 0,0391 
15 — 0,1181 3,5 0,0168 
1,6 — 0,0337 3,6 - 0,0215 
‘= + 0,0468 3,7 — 0,0311 
18 0,1011 3,8 — 0,0052 
19 0,1172 3,9 + 0,0260 


3. Die plötzliche Belastung 


Die symbolische Lösung der durch einen Kraftsprung hervor- 
gerufenen Welle lautet 


(95) Yy = sinnzyp)- 


Sm®EJ.pVp be SnEJ Vp = 


-cosan Vp - (de). 


Dann läßt sich die Deutung in derselben Weise wie oben vollziehen, 


> 
y= 


(96) 


es . 
und man erhält zunächst mit z = —\— den Ausdruck: 


Vat 


Mer tf: 1 
2nEJ \6 2a + — 8, 


= 
En 
% 
Fr 
= 
= 
5 
ee Zur Interpretation stellt man sie als dreifaches Integral dar: -_ 
47) 
- 
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Hier bedeutet 


ebenso 


3 2 z 1 
S,(2) = 2 — 1) #7; (2) — -sinz2?. 


Nach Einsetzen dieser Funktionen ergibt sich mit 


| = (1 — C(z) — S(@)) + (sin 2 — cos— 2? 


2 2 
(98) 
| + S22. —C(e) + cost + sin 
die Lösung 


6n® EI Vn ‘ 

Die Hilfsfunktion F, wurde eben- 
falls tabuliert und in Fig. 10 auf- 
gezeichnet. Durch Dilatation der 


¥ (2, t) = F,(2)- 
bee) 


Figur nach dem Gesetz Yt in der 


z-Achse und tt in der F-Achse 
erhält man die Kinematik. Fig. 10 onal as 


eu 
z F, (2) z (z) 2 
s 0,0 + 1,0000 1,0 — 0,1157 
0,1 0,9577 1,1 — 0,1264 
0,2 0,8490 1,2 
0,4 0,5298 1,4 — 0,0197 Bu 
0,5 0,3598 1,5 + 0,011 
0,6 0,2009 1,6 0,025 
0,7 0,0732 1,7 0,036 
0,8 - 0,0255 1,8 0,026 Leer 
an 0,9 — 0,0870 1,9 0,010 nz : 
2 


b) Bestimmung der Eigenfrequenzen einer Klaviersaite 
während des Anschlags 


‘Oo 

8 


Exe 


Es ist eine Saite gegeben, die an den Stellen z, =a und 
2, =— 6a befestigt und an der Stelle x, = 0 mit einem Hammer 


gekoppelt ist, dessen Befilzung als rein elastisch angesehen werden 


z z . a 2 z rt 2 
C(z) — — - S(z)+ -|1— 008— 2 
( ) 6a 2 2a ( ) 3 n? 2 
=) 
7 
29 
|. 
= 
be 
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soll. Der Anschlag erfolgt, indem man der Hammermasse plötzlich 
den Impuls I erteilt. 

Das Schaltschema der mechanischen Widerstände ergibt hier 
folgendes Bild (Fig. 11): 

Die beiderseits befestigte Saite hat im Punkte 2=0 den 
Widerstand R,. Hierauf wirkt die Stoßkraft erst durch die in Serie 
geschaltete Befilzung R,. Dem ganzen ist die Hammermasse R, 
parallel geschaltet. Die Stoßkraft pI wirkt demnach auf den resul- 
tierenden Senn 

#3 R,R, + + BR, 
und erzeugt die Auslenkung 

y, (R, + R)-pI Am 

des Klavieranschlags Dabei ist 

(99) R, =m p?, m = Masse des Hammers, 


A, 


(100) R, = C = const. = Filzelastizität bei Druck gegen die Saite 
zu setzen. 

Den Widerstand R, kann man aus Gl. (69) entnehmen, wenn 
man einsetzt 
a, = Saitenspannung S, a,=yFp? (y = Dichte, F = Querschnitt) 
und daher 


W=S#, C= 
ferner 
b = 6a, r,=—l=r, 
Man findet so 
- 14a? 
R, = 28 
(101) 


{Gtg a P+ 


Bezeichnet man die Verformung der Widerstände R, und R, mit y, 
bzw. y,, so findet man aus den Beziehungen der Serienschaltung 


R, = By ys 
die der Saite im Punkte 
| 
: 
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oder, nach Einsetzen von y, 


R,-pI 
(103) Ys R, BR + R,R, 
Führt man die Abkürzung Be 
_ R,R, +R, R, +R, RB, 
= mp! + (1+ 5p°) 52. (Giga? + Cig6a?) 


ein, so läßt sich die Bewegung der Saite im Anschlagpunkt für die 
Dauer der Berührung in die Gestalt 


ts 


bringen. 

Durch Partialbruchzerlegung erhält man, wenn p, die einfachen 
Nullstellen von R(p) sind, 
(106) R(p,) = 0, ¥ 

(107) = 2 (p — RB’ (p,) > 
(n) (n) 


Da die Rechnung dämpfungsfrei angesetzt wurde, sind die p, rein 
imaginär, und man erhält die Eigenfrequenzen «, (Kreisfrequenzen), 
indem man 

?, = 
setzt. Wegen des transzendenten Charakters der Gleichung 

Ri a.) = 0 
empfiehlt es sich, die Lösung graphisch zu ermitteln. Dabei soll als 
Einheit der Frequenz die Grund-(kreis)-frequenz der frei schwingen- 
den Klaviersaite 


c nc 

{108) 
eingeführt werden, indem | 

(109) 


gesetzt wird. “Man erhält dann nach elementarer Umrechnun 
Bestimmungsgleichung 


(110) ctg s, + ctg 


Dabei ist 
Tal 
nas 
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In Fig. 12 wurden die Hammerfunktionen 


os 


9 
6 


62 
ctg - 


mit ihren Asymptoten über s aufgetragen und die Schnittpunkte der 
Kurven angezeichnet. Dabei wurden nachfolgende Daten zugrunde gelegt, 
die sich beim normalen Hammer auf Messungen stützen, und es er- 
gaben sich die zu unterst aufgeführten Eigenfrequenzen (mit dem 
Grundton der freischwingenden Saite als Einheit). 

Die durch den Bogen ) verbundenen Werte können als Auf- 
spaltung einer Eigenfrequenz der freischwingenden Saite aufgefaßt 
werden. Man erkennt aus der Figur, daß im Falle des Zusammen- 
treffens der Asymptote der Hammerfunktion mit einer solchen der 
Saitenfunktion — d.h. wenn sy ein ganzes Vielfaches von 7/6 ist — 


keine Aufspaltung erfolgt. mate 


N 
MW 7 AQ Z 
| 
| 
| 
X 
Er 
4 5 -¢ | | 
|| 
A \ Ä 
W 2 A 
Fig. 12 
a für vier verschiedene Wertpaare o, sy sowie die Saitenfunktion 
Zi 
| 
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Hammer weich (w) | normal (n) | hart (h) hen in 
Bezeichnung 
der Schnitte 4 
m=ing 9,32 9,32 9,32 0,932 
4,488 11,22 22,44 11,22 
2 5 10 5 
1,2 3 6 30 
a 0,825 0,850 0,995 
1,194) * 1,280, 1,355 1,920 
2,150 2,255 Hel 2,815 
3,175 3,335 3,440 3,748 
4,155 4,365 4,540 4,763 
5,095 5,278 5,535 5'768) 
6,025 6,076 6,181 | 6,337 
Berechnet für S = 50 kg-Gew., yF = 0,25 g/em, a = 10 cm. 


Die beiden mit * bezeichneten Eigenfrequenzen wurden berechnet. 


Zusammenfassung 


Ausgehend davon, daß die linearen Differentialoperatoren alge- 
braischen Gesetzen genügen, wurde im Anschluß an den Heaviside- 
Kalkül eine synthetische Methode zur Berechnung linearer zeitlich 
unveränderlicher mechanischer Systeme entwickelt, die vermöge der 
allgemeinen Eigenschaften der zugrunde liegenden Integraltransfor- 
mation auf das wesentliche Verhalten der betrachteten Systeme 
schließen läßt, ohne (im Gegensatz zu der sonst üblichen Operator- 
rechnung) von der direkten Interpretation Gebrauch machen zu müssen. 

An dem einfachen Beispiel der Bewegungen eines Stabes (Typ 
der symbolischen Differentialgleichung: 


4, ¥"(%,P) + (P) + (2, p) = 0, 
wo a, und a, von x unabhängig sind), mit Randbedingungen, welche 
die Zeit (d. i. die zu p konjugierte Variable t) nicht explizit ent- 
halten, wird das Verfahren systematisch durchgeführt. Zum Schluß : 
werden einige numerische Beispiele berechnet. a 


Herrn Prof. Dr.-Ing. M. Schuler bin ich fiir Férderung dieser 
Arbeit zu groBem Dank verpflichtet. ER 


Göttingen, Institut für angewandte Mechanik. 
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Von G. Schaitberger 


(Mit 9 Figuren) 
ae ae Die niedrigste Spannung, bei der noch Röntgenstrahlen einer 
vorgegebenen Frequenz entstehen oder die Maximalfrequenz, die bei 


einer konstanten Röhrenspannung erzeugt wird, _ dem Quanten- 


Diese Beziehung ist neben der Bestimmung aus der Bohrschen 
Formel der Rydbergkonstante die zuverlässigste Methode zur Be- 
stimmung des Planckschen Wirkungsquantums. Die genauesten 
Bestimmungen von Ah durch Messung der kurzwelligen Grenze 
stammen von W. Duane!) und seinen Mitarbeitern und von 
H. Feder’). Leider stimmen die Resultate der beiden Autoren 
nicht überein. Die Abweichung liegt außerhalb der angegebenen 
Fehlergrenze. Deshalb konnte eine Neubestimmung von Vorteil 
sein. Dabei kam die sogenannte Isochromatenmethode in An- 
wendung. D.h., es wurden Kristall und Ionisationskammer am 
Braggspektrographen auf eine bestimmte Wellenlänge eingestellt und 
die Intensität der reflektierten Strahlung bei verschiedenen Span- 
nungen in der Nähe der kurzwelligen Grenze gemessen. Aus den 
so erhaltenen Kurven wurde die kurzwellige Grenze durch Inter- 
polation erhalten. Die experimentellen Schwierigkeiten bestanden 
in sehr genauer Hochspannungsmessung und in einer Aa genauen 


Bestimmung der Wellenlänge. En. 


i Hochspannungsquelle und § Spannungsmessung 


Als Spannungsquelle diente nicht wie bei Feder oder Duane 
eine Hochspannungsakkumulatorenbatterie sondern eine Stabilivolt- 
anlage von Siemens, Reiniger und Veifa, die speziell fiir hohe 
Spannungskonstanz gebaut war. Diese Anlage wurde mit Wechsel- 
strom von 400—500 Perioden gespeist. Der Wechselstrom wurde 


1) W. Duane, H.H. Palmer u. Chi-sun- Yeh, Journ. Opt. Soc. Amer. 5. 
S. 213. 1921. 
2) H. Feder, Ann. d. Phys. [5] 1. S. 497. 1929. 
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einem alten Schuckertumformeraggregat von sehr solider Bauart 
entnommen, das mit der 70 Volt-Hausbatterie betrieben wurde. 
Diese alte Maschine zeichnete sich durch große Schwungmasse aus, 
die einen gleichmäßigen Gang garantierte. Das Stadtnetz konnte 
wegen seiner Inkonstanz, die sich natürlich auch auf die Hoch- 
spannungsseite der Stabilivoltanlage übertrug, nicht verwendet werden. 
Das Umformeraggregat hatte außerdem wegen der ungefähr 10 mal 
größeren Periodenzahl des erzeugten Wechselstroms den Vorteil zu 
einem ausgeglicheneren Hochspannungsgleichstrom beizutragen. 

Das Stabilivoltprinzip beruht darauf, daß in jeder Phase 
Kondensatoren aufgeladen werden, die sich über das Röntgenrohr 
teilweise entladen. Die Spannung am Rohr sinkt deshalb nach jeder 
Periode herunter, und zwar proportional der entnommenen Strom- 
stärke. Eine Stabilivoltanlage liefert also keine konstante Gleich- 
spannung, sondern eine Gleichspannung, der eine Wechselspannung 
überlagert ist. Um den eventuellen Einfluß auf das Endresultat 
durch Verschiebung der kurzwelligen Grenze festzustellen, mußte 
die Größenordnung der Schwankungen bestimmt werden. Bei 
50 Perioden Wechselstrombetrieb und 20 Milliampere Belastung, der 
größten bei den Messungen verwendeten Stromstärke, zeigte sich am 
Elektrometerfaden der im folgenden beschriebenen Spannungsmeß- 
anordnung eine Unschärfe. An die Hochspannungsseite der Spannungs- 
meßanlage wurde dann eine kleine Wechselspannung gelegt, die die- 
selbe Unschärfe hervorrief. Hieraus konnte die Größe der über- 
lagerten Wechselspannung geschätzt werden. Die Schwankungen 
waren kleiner als 15 Volt und wahrscheinlich größer als 5 Volt. 
Dies deckt sich mit den Messungen an der Stabilivoltanlage der 
Münchener Technischen Hochschule’), die ähnlich gebaut war. Das 
Siemens Wernerwerk, das diese Messungen ausführte, fand pro Milli- 
ampere Belastung ein halbes Volt Schwankung. Nach einer Näherungs- 
formel von Webster?) sind die Schwankungen umgekehrt pro- 
portional der 3. Potenz der Frequenz des verwendeten Wechsel- 
stroms. Bei 400—500 Perioden waren sie also einige hundertstel 
Volt groß. Gegenüber der verwendeten Gleichspannung von min- 
destens 9kV kann man die überlagerte Wechselspannung unbedenklich 
vernachlässigen. 

Die Hochspannung wurde mittels einer Influenzmethode ge- 
messen. Zwei Kondensatoren, ein großer und ein kleiner, waren 
elektrostatisch geschützt mit einem Fadenelektrometer verbunden 


. Kulenkampff, Ann. d. Phys. 87. S. 597. 1928. 
. L. Webster, Proc. Amer. Nat. Acad. 6. S. 270. 1920. 
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(Fig. 1). Intluenziert man gleichzeitig mit einer Wippe über den 
kleinen Kondensator die Hochspannung, über den großen eine 
Niederspannung mit entgegengesetzten Vorzeichen, so verhalten sich 
die Spannungen umgekehrt wie die Kondensatoren, wenn sich am 

Klektrometer kein Ausschlag zeigt. 


a A 8 Als groBer Kondensator fiir die Nie- 
| derspannung wurde ein Harmskonden- 
f\ sator verwendet; der kleine Konden- 


sator (Fig. 2) wurde mehrmals abge- 


C 
: 7] | ändert, bis es sich zeigte, daß das 


-_ Kondensatorenverhältnis temperatur- 
. D unabhiingig war. Die Niederspannung 
für den großen Kondensator wurde 
r einem Potentiometer entnommen, an 
dem eine Akkumulatorenbatterie lag. 
Hochspannung Meder spa nung Das Potentiometer, ein Wolff- und 


ins) > rheosts armoctlic rleich- 
ein Siemensrheostat ermöglichte gleich 


8 Harmshondensator zeitig die Niederspannungsmessung 
C fsoenelektrometer ER 
D Winpe durch Vergleich mit 5 Normalelemen- 


ten. Diese Normalelemente wären 
mit einem geprüften Original Weston- 
element verglichen worden. Die Ab- 
weichung + 0,12°/,, ist im Endresultat berücksichtigt. Die Spannungs- 
meBanlage wurde mit einem genau bekannten Spannungsverhältnis ge- 
eicht, das einem Potentiometer entnommen wurde. Die Kondensatoren 


Fig. 1. Schaltung der Spannungs- 


meBanordnung 


E Hochsoan- starre 

S/H/1SCh 
— dung mitdern Eek - 

- frometer 


Zsen 

Anzeilan 

Messing 
Harlgumm 


Fig. 2. Hochspannungskondensator zur Spannungsmessung 


verhielten sich dann bei Elektrometerausschlag Null umgekehrt wie die 
Widerstände. Am Potentiometer lagen 2—5.6kV. Als Widerstände 
wurden zwei auf Glas gewickelte Manganindrahtwiderstände von 
zusammen über 100000 Ohm und die oben erwähnten Rheostaten 
verwendet. Das Verhältnis der Widerstände zueinander wurde 


durch Vergleich mit einem 1000 Ohm-Normal genau festgelegt. 


| 
va 
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Für das Kondensatorenverhältnis wurde zu verschiedenen Zeiten 


erhalten 


1:91,94 

1:91,98, 
1:91,92, 
1:91,97, 


| Mittel: 1:91,95, 


| Die Unsicherheit, die durch die Spannungsmessung entsteht, soll mit 


2 0/ 1 ‚rer 
eingesetzt werden. 


Röntgenrohr 
7 Als Röntgenröhre wurde ein Elektronenrohr, das dauernd an 
: der Diffusionspumpe!) lag, verwendet. Die Glühkathode war die- 
| selbe wie am Ottschen Rohr?). Die Heizung des Glühfadens wurde 
einer isoliert aufgestellten Akkumulatorenbatterie von 8 Volt ent- 


nommen. Die Antikathode und das Gehäuse lagen an Krde. Alle 
auswechselbaren Teile: Anode, Kathode und Fensterhaltevorrichtung 
waren mit Bleidichtungen versehen. Zum Strahlenaustrittsfenster 
wurden Cellophan und Aluminiumfolie verwendet. Zwei "°/,,, mm 
starke Cellophanfolien und eine '/,,, mm starke Aluminiumfolie 
wurden mit dünner Schellacklösung zusammengeklebt. Dabei diente 
das Aluminium nur zum Schutz gegen vagabundierende Elektronen. 
Wegen der niedrigen Ordnungszahl der Elemente, aus denen Cellophan 
besteht, kam für die Absorptionsverluste praktisch nur das !/,,, mm 
| starke Aluminium in Frage. 


Spektrograph 
Der Spektrograph war derselbe, den Feder bei seiner h-Be- 
stimmung verwendet hatte. Die Spaltbreiten waren sehr gering, um 
ein möglichst paralleles Röntgenstrahlenbündel zu erhalten. Der 
Spalt bei der Ionisationskammer S, war 0,12 mm, S, war 0,07 mm 
breit (Fig. 3). Die Braggsche Fokussierungsbedingung, d. h. der 
‚gleiche Abstand der Spalte von der Spektrometerachse, war sehr 
genau eingehalten. Als Kristall wurde ein ausgesuchtes Steinsalz 
von Seemann in Freiburg verwendet. Jeder Quadratmillimeter war 
von Seemann nach seiner Schneidenmethode geprüft worden. Um 
die reflektierende Schicht nicht durch die Luftfeuchtigkeit zu 
‘schädigen, wurden folgende Vorsichtsmaßregeln angewendet. Der 


1) Für die zum Ausfrieren der Hg-Dämpfe verwendete flüssige Luft sage 
ich den Vereinigten Sauerstoffwerken Nürnberg meinen Dank. 

2) H. Ott, Phys. Ztschr. 27. S.598. 1926. Herrn Prof. H. Ott danke ich 
an dieser Stelle für wertvolle Anregungen beim Bau des Röntgenrohres. 
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Kristallkopf war als doppelter Exsikkator ausgebildet. Die äußere 
Glocke wurde während der Messung abgehoben. Der innere Exsikkator 
war mit Cellophanfenster für den Durchtritt der Strahlen versehen. 


: “4 | 
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A Cellophanfenster 


Ein EinfluB des Bielibeten durch D Quarzschliff 
Streustrahlung usw. konnte nicht E Gummidichtung 
festgestellt werden. F Innenelektrode 


G Potentiometer 

H 100-Volt-Akkumulatorenbatterie 

Fig. 4. Kompensationsionisations- 
kammer (schematisch) 


Auf demselben Spektrogra- 
phenarm, auf dem sich S, befand, 
waren die Ionisationskammer und 
das Fadenelektrometer montiert 
(Fig. 4). Die Ionisationskammer bestand aus zwei parallel ge- 
schalteten Kammern, von denen nur die eine die Röntgenstrahlen 
aufnahm; die andere mit entgegengesetzter Gehäusespannung diente 
zur Kompensation des geringen von der radioaktiven Strahlung der 
Umgebung herrührenden Restionisationsstroms. Die Innenelektrode 
war mit einem doppelten Quarzschliff eingeführt. Als Füllung 
wurde Methyljodiddampf verwendet. Die Ladungsempfindlichkeit des 
Elektrometers konnte mit einem Influenzierungsring!) bei jeder 
Messung leicht bestimmt werden (etwa 25 Skt. pro Volt - cm). 

Die Röntgenstrahlen und die reflektierende Kristallfläche waren 
nach den bekannten Methoden in die Drehachse justiert. Der Teil- 
kreis des Spektrographen war in ganze Grad geteilt; mit Hilfe eines 
20 teiligen Nonius konnte man 3 Min. ablesen, 1 Min. schätzen. Die 
Teilung wies Fehler auf, sie wurde deshalb nur zur Vorjustierung 


i > 
Ya 
= 
| 
i 
AR 
/ 
. 
A 
_ 
. 
‘© 
Lt 
¢ 4 
a 
mee 


G. Schaitberger. Ein Beitrag zur h-Bestimmung 


benützt. Da die Wellenlängen der Röntgenlinien, hauptsächlich von 
der Siegbahnschen Schule, mit einer relativen Genauigkeit von 
1:100000 bestimmt sind, wurden sie zur Frequenzbestimmung 
herangezogen. Kristallkopf und Spektrographenarm mit Kammer 


und S, konnten mit Mikrometer- 
schrauben gedreht werden. Mit 
diesen Schrauben konnte auf die 
Mitte der Linien einjustiert werden 

(Fig. 5). Als Linien wurden die 

L-Linien der Antikathode ver- 

wendet. Wolfram rührte von der 

Zerstäubung des Glühfadens her. 

Platin wurde galvanisch aufgetra- 

gen. Die Bleilinien wurden da- 

durch erhalten, daß die Anti- 
kathode mit einem stark bleihalti- 

gen Lötzinn verzinnt wurde. Der 

Fehler bei der Einjustierung in 
die Linie wird auf einige Sekunden 
geschätzt. Selbst bei den kleinsten 
doppelten Glanzwinkel kann man 

ihn unbedenklich vernachlässigen. 
Das schwere Elektrometer 

mit der Kompensationsionisations- 
kammer auf dem Spektrographen- 
arm ließ befürchten, es könnten 
_ Deformationen eintreten, die die 
Justierung gefährden würden. 

Ebenso könnten vielleicht durch 


vayıayurz UI 


/rommelteile 


l l 1 l 


9D OV HO 


Fig. 5. Justierung 
in die Linie WLy, = 1096,30 XE. 
Letzte Durchdrehung 
des Ionisationsspektrographenarmes 
vor und nach der MeBreihe. 
100 Trommelteile entsprechen 9 Min. 


das Abheben und Daraufsetzen 


der Exsikkatorglocke Verschiebungen eintreten. Die Justierung in 
die Linie wurde deshalb nach jeder Meßreihe geprüft. Abweichungen 
außerhalb der Meßfehler (einige Sekunden) konnten nicht festgestellt 


werden. 


Der Steinsalzkristall hatte einen Temperaturkoeffizienten, der 

die Gitterkonstante und damit in diesem Falle die Wellenlänge, die 

am Spektrographen eingestellt war, veränderte. Um diesem Fehler 
Rechnung zu tragen, soll vorsichtshalber als Unsicherheit + 0,04°/,, 

entsprechend einer Unsicherheit von 1°C des Mittelwertes der Kri- 


stalltemperatur eingesetzt werden. 


Daß sich bei Wiederholung der 


_ Justierung in die Linie keine wesentlichen Unterschiede zeigten, be- 
_ rechtigt zu der Annahme, daß die Temperatur während der Ver- 


suche genügend konstant war. 
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Gang der Messung 

Für die Isochromaten mußte die Intensität der vom Kristall 
reflektierten Strahlung bei einer bestimmten Wellenlänge in Ab- 
hängigkeit von der Spannung am Röntgenrohr gemessen werden. 

A) Zur Messung der reflektierten Strahlung mit Ionisations- 
kammer und Elektrometer waren erforderlich: 

1. Bestimmung der Nullpunktsverschiebung (Gang des Elektro- 
meters) vor der eigentlichen Messung. 

2. Bestimmung der Ladungsempfindlichkeit des Elektrometers 
mit des Influenzierungsringes. 

3. Ablesung des durch die reflektierten Röntgenstrahlen ent- 
sta denen Ionisationsstromes am Elektrometer. 
4. Ablesung des Entladungsstromes durch das Röntgenrohr. 

5. Wiederholung der Gangbestimmung nach der Messung. 
(In vielen Fällen als 1 des nächsten Kurvenpunktes verwendet.) 

B) Zur Messung der Spannung mit der Spannungsmeßanordnung 
waren erforderlich: 

6. Einstellung einer Niederspannung für den Harmskondensator. 

7. Messung dieser Niederspannung, die mit dem Kondensatoren- 
verhältnis multipliziert den Wert der Hochspannung ergab, durch 
Vergleich mit den 5 Normalelementen. 

8. Influenzierung der Hoch- und Niederspannung über die 
Kondensatoren und Einregulieren auf Elektrometerausschlag 0 

9. Kontrolle der Messung der Niederspannung. 

Diese Messung für einen Punkt der Isochromate wurde zeitlich 
folgendermaßen ausgeführt): 

Während der Gangbestimmung (1.) wurde eine bestimmte Nieder- 
spannung für die Spannungsmeßanordnung eingestellt (6) und ge- 
messen (6.). Dann wurde dieWippe der Hochspannungsmeßeinrichtung 
herumgeworfen und die Hochspannung durch einen in Serie mit 
Stabilivolt und Hochfrequenzmaschine liegenden Widerstand auf 
Elektrometerausschlag 0 einreguliert (8). Nach Ablauf der 2 Min. 
dauernden Gangbestimmung konnte durch Umlegen eines Schalters 
in wenigen Sekunden die Ladungsempfindlichkeit des Ionisations- 
kammerelektrometers bestimmt werden (2). Sofort anschließend 
wurde die Bleiblende am Röntgenrohr geöffnet und 2 Min. lang 
Röntgenstrahlen austreten lassen. Während dieser Zeit wurde die 
Spannung am Rohr dauernd beobachtet. Das Elektrometer der 


1) Wegen der Gleichzeitigkeit der Teilmessungen konnten die Ablesungen 
nicht allein ausgeführt werden, deshalb wurde das Ionisationskammerelektro- 
meter (1, 2, 3, 5) von den Herren K. Dietz oder Dipl.-Ing. W. Franz ab- 
gelesen, denen an dieser Stelle gedankt sei. fs 
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SpannungsmeBanordung reagiert ja fast trägheitslos auf Spannungs- 
schwankungen. Mit dem oben erwähnten Widerstand wurde die 
Spannung dauernd nachreguliert und auf ein halbes Promille genau 
konstant gehalten. Nach diesen 2 Min. wurde der Réntgenstrahl 
wieder abgeblendet, die Milliampere des Röhrenstromes abgelesen (4.), 
die Wippe zurückgeworfen. Das Elektrometer mußte jetzt wiederum 
auf Null stehen, sonst wurde die Messung verworfen. Gleichzeitig mit 
Schließung der Bleiblende wurde am Jonisationskammerelektrometer 
der Ausschlag abgelesen (3.) und anschließend wieder 2 Min. lang 
der Gang bestimmt (5... Da der Gang nach der Messung in vielen 
Fällen bei rascher Aufeinanderfolge der Messungen gleich als Gang 
vor der Messung für den nächsten Kurvenpunkt diente, nahm die 
Zu Tab. 1 
A Ausschlag 
G Mittelwert des Ganges vor und SH -— 95000 Ohm” 
nach der Messung 
mA Röhrenstrom 
E Empfindlichkeit des Elektrometers 
KV ausgerechnete Spannung 
Ak= 25 - 10°, korr. Ausschlag 
') Angabe der Ohmwerte noch ohne 
Korrektion der Fehler der Wider- 


standskästen 


Tabelle 1 


| mA | E. SK. 


26,8 | 3633,1 
26,6 | 3633,0 
26,3 | 3615,6 
26.4 | 3615.6 
263 | 3598.0 
26.4 
26.4 
26,3 
26,5 
26,6 
26,9 562,4 | 12,964 
272 | 3562. 

27,5 
28,0 13,026 
26,0 | 

26,7 
27,9 
26,6 
25,7 


G0 Go te Go GO GO GD GO 


or 


or 


WO Wwe 


2222 


BANNAN ANN 


wot 
| 
bo 


13,090 


- 


13,154 


Oa 


| 
| 
W i 
| 
" 
| if 
73 -0,5 7 
31-01 | -0,1 
4) —0,3 | +0,1 
14 
7| -05 | —03 
9 +02 | -0,35 
1-18) -0,65 
121 -21 | -08 
23) | -0,05 
215) -39 | —0,05 
1216| -30 | -0 Wa 
17) -4,9 | —0,2 
1318| —5,7 | -0,35 
= 19 —6,0 | —0,25 a 
20, -50 | +01 | 
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Messung eines Kurvenpunktes nur 4!/, Min., im Höchstfall bei noch- 
maliger Gangbestimmung 6—7 Min. in Anspruch. Die Kurven sind 
teils mit steigender, teils mit fallender Spannung aufgenommen 
worden. Ein Beispiel einer Messung bringt Tab. 1. Figg. 6 und 7 
stellen Ausschnitte aus Isochromaten dar. Die 3. Kurve auf Fig. 6 
ist die zu Tab. 1 gehörende. Der größte Fehler bei der Bestimmung 
der Grenzspannung aus den Isochromaten dürfte von der Unsicher- 
heit der Gangbestimmung des Elektrometers herrühren. Bei der 
äußersten Elektrometerempfindlichkeit, die durch die Enge der Spalte 
gefordert war, waren die Messungen überhaupt erst nach Einbau der 
S. 88 erwähnten Kompensationsionisationskammer mit der erforder- 
lichen Genauigkeit möglich. Auch bei geringerer Elektrometer- 
empfindlichkeit waren die Schwankungen des Ganges prozentuell 
die gleichen. 

An die aus den Kurven erhaltenen Grenzspannungen (Tab. 2) 
sind einige Korrektionen anzubringen. 

1. Die eine Korrektion war durch die Spaltweite bedingt. In- 
folge der Spaltweite trat keine bestimmte Wellenlänge in die 
lonisationskammer ein, sondern ein ganzes Wellenband. Um diesen 


Wolfram 


(V—) Spannung in kV gemessen mit Voltmeßanlage am negativen Heiz- 
fadenende. 
(V+) Spannung in kV gemessen mit VoltmeBanlage am positiven Heiz- 
fadenende. 
D Differenz in Volt zwischen den Mittelwerten (V—) und (V+). 
(Die Abweichungen, auch des nach der umgekehrten Seite liegenden 
j Wertes --2,3, entsprechen den Meßfehlern.) 
_ MMittelwert aus den Mittelwerten (V—) und (V+)). 


9,043 9,631 | 
9,6491 9,649 963% | 9,644, 
1279,17 9,655 
Wolfram Ly 11,254 
11249 || 11,253, 11,256 -2,3 11,254, 
1096,30 11,258 | | 
Platin Ly, | 12,920 12,891 a 
12,911, { 2802 | 12,891, 19,8 12,901, 
955,99 12,903 12,892 an 
Blei Ly, | 14715 14,721 
14720 | 14,724, | 14,718, 6,4 14,721, 
838,01 14,729 14,701 


kin XE (V-) (V+) D M 
Mittel Mittel 
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Einfluß zu eliminieren, wurden Isochromaten mit verschiedener Spalt- 
weite aufgenommen und auf Spaltweite 0 extrapoliert (Fig. 8). Als 
Wert für die Korrektion ergab sich 0,57°/,, mit einem möglichen 
Fetter von + 0,03°/,,. 
2. Eine andere Korrektion bedingte der Glühfaden. Die Elek- 
tronen traten hauptsächlich von der heißesten Stelle des Glüh- 
_ fadens, der Mitte aus (Fig. 9). Je nachdem das Voltmeter bzw. die 


= 
oe 


l 
S45 6 


Spaltweite 

Fig. 8. Extrapolation a Fig. 9. Voltmeter i 
auf Spaltweite 0 zeigt in diesem Falle zu niedrig 


Spannungsmeßanlage, am negativen oder positiven Ende des 
_ fadens lag, wurde die Voltgeschwindigkeit der austretenden Elek- fae 
- tronen zu hoch oder zu niedrig gemessen. Die Differenz der Volt- wa 
__ geschwindigkeiten muß gleich der Heizfadenspannung (8 Volt) sein. 


A — im Mittel 8,4 Volt erhalten. Es wurde an- 


hebt. Außer der das angelegte elektrische Feld 
p Geschwindigkeit haben die Elektronen auch noch eine thermische 
Geschwindigkeit (ungefähr 0,1—0,3 Volt). Die thermische Ge- 
schwindigkeitkeit kann gegenüber der anderen vernachlässigt werden. 
Die h-Werte für die verschiedenen Wellenlängen sind in Tab. 3 


enthalten. Für die Berechnung wurde verwendet 


2d = 5628 XE in den Siegbahnschen Wel- 
 lenlängen der Linien enthalten. 

Als Mittelwert ergibt sich h = 6,553, + 0,002 - 10727 Erg/sec. 
n In der Fehlergrenze sind enthalten + 0,2°/,, Spannungsmeßfehler, 
+ 0,03 °/,, Spaltweitenkorrektionsfehler; 0,07°/,, entstanden durch 
den Temperaturkoeffizienten des Kristalls usw. Nicht berücksichtigt aa 
sind die Fehler von d, e und ce. 
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2 Tabelle 3 
Mittelwert der Grenzspannung | Planeksches 
Wellenlänge in AE in kV | Wirkungsquantum - 10? 
Wolfram L | 
1279,17 9,644, 6,553, 
Wolfram L 
12901 6,552, 
838.01 14,721, 6,553, 


Zum Vergleich seien die anderen Werte des Planck schen 
Wirkungsquantums nach der Methode der kurzwelligen Grenz- 
bestimmung der Röntgenstrahlen angeführt. Beide Werte sind hier 
ebenfalls mit 2 d = 5,628 XE, e = 4,774-10—" und c = 2,9986.10" 
berechnet: 
Duane und Mitarbeiter’) 6,561 + 0,003. 10=?" Erg/sec 
Feder. . ... . . 6,547 + 0,003 - 1072" Erg/sec. 
In den so berechneten Werten steckt natürlich die Unsicherheit 
der Gitterkonstante, der Voltangabe (internationale Einheit), des 
Elementarquantums und der Lichtgeschwindigkeit. Die Licht- 
geschwindigkeit kann mit c = 2,99774 . 10192) für den hier in Betracht 
kommenden Fall als genau genug bestimmt angenommen werden; 
eine Neubestimmung würde ohne wesentlichen Einfluß auf das End- 
resultat sein. Der Fehler der Voltangabe (internationale Einheit), 
der von einem Fehler der internationalen Widerstandseinheit her- 


Br: a rührt, kann durch Multiplikation mit 1,0005 beseitigt werden. Die 
zu 2: Angabe erfolgt nunmehr in absoluten Volt. Die verbleibende Un- 
ae ow sicherheit rührt im wesentlichen von der Gitterkonstante und vom 

Elementarquantum her. 

> Als Gitterkonstante hat Siegbahn den Mosleyschen Wert 
Br 7 ou d = 2,814 XE fiir Steinsalz genommen. Da die Genauigkeit wegen 
der hohen relativen Genauigkeit der Réntgenlinien (heute 1:100000) 

nicht ausreichte, wurden dem Mosleyschen Wert zwei Nullen an- 


le gehängt und dieses d als Gitterkonstantennormal definiert. Alle 
 Wellenlängenangaben der Röntgenlinien enthalten diesen Wert. 


Br, ®: 1) Journ. Opt. Soc. Amer. 5. S. 213. 1921. Der von Duane angegebene 
Wert 6,556-10-%7 ist von ihm mit 2d,, = 6,056 (Kalkspat), e = 4,774-1071°, 
Bere: c = 2,9986-10'° gerechnet. Zur Umrechnung von elektromagnetischen in elek- 
trostatische Volt wurde offenbar c = 3,00-10'° verwendet. 

2) A. A. Michelson, Pease und Pearson, Science Service in Sci- 
ence 78. Dec. 29. 1933. _ 
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Weil es die gebräuchlichste Gitterkonstante ist, sei es auch hier 2 
verwendet. Es soll aber dadurch kein Entscheid über den 
wirklichen Wert von d getroffen werden zwischen dem Sieg- 
bahnschen Wert, dem Wert von Baerden (aus Strichgitter- 
reflexionen berechnet), dem von Kirkpatrick und Roß!) (nach 
dem Gesetz von Duane und Hunt unter Verwendung eines wahr- _ ; 
scheinlichen h/e-Wertes berechnet) und den anderen Gitterkonstanten- _ = oe 
angaben. Das Siegbahnsche d wird in dieser Arbeit zwecks — ; 
leichterer späterer Korrektur mit dem unbekannten Faktor & ver- 
sehen, der angeben soll, um wieviel das konventionelle Siegbahn- 
sche d von dem „wirklichen Wert“ der Gitterkonstante entfernt ist. _ 

Bei Vergleich von Konstanten findet man öfter diese Kon- 
stanten an irgendeiner versteckten Stelle mit verschiedenen Werten 
anderer Konstanten gerechnet. Aus Gründen der Übersicht über — 
die verwendeten Konstanten und zwecks späterer leichterer Um- | 
rechnungsmöglichkeit erfolgt in Tab. 4 die Angabe des Endresultates | 
he vi Al „de = 


Auf der rechten Seite dee Gleichung steht das soins a el 
bestimmte V in absoluten Volt, das dabei verwendete 2 (mit dem : 
konventionellen Siegbahnschen d gerechnet) und die Lichtgeschwin- . 
digkeit c. Auf die linke Seite der Gleichung wurden zu dem ge- = 
suchten 4 die Größen gebracht, die wegen ihrer Unsicherheit h - DB 
verfälschen können und infolgedessen einer späteren Korrektur aus- 
gesetzt sind. Es handelt sich um das verhältnismäßig 7 
Elementarquantum und um «, den Faktor, der das konventionelle — 


- 
abelle 
h _ VA 
em e 


1,3726 + 0,0002. 1077 


Peder ....<.... 1,3727 + 0,0008-10-"7 
Schaitberger.... + 0,0006 - 1077 


1,3759 + 0,0008-10-"" 


wobei 2 mit 2,81400 XE berechnet, V in absoluten statischen Volt und c mit = 
2,99774-10" angegeben ist. « ist ein unbekannter Faktor, der das konventionelle 
d von Siegbahn zum „wirklichen d“ ergänzt. u 


1) P. Kirkpatrick u. P. A. Roß, Phys. Rev. 45. S. 454. 1934. 
2) P. Kirkpatrick u. P. A. RoB, Phys. Rev. 45. S.454. 1934. 1,3726 wird 
erhalten nach Abzug einer fraglichen Korrektur von 4,5 Volt für die An- 


ziehungskraft der Antikathode auf die Elektronen, sonst 1,3729. ts un 
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Wegen der Unsicherheit des Elementarquantums 
‚= 4,770 + 0,005 . 10-103) 
soll die Berechnung von h aus VA/c nur der Vollständigkeit halber 
in Tab. 5 folgen. 
Tabelle 5 


6,5473 + 00006-1077 
6,548 +0,003 +1077 

Behaitberger.... | 6,555 £0,002 


mit e = 4,770-10~ 1° aus VA/c von Tab. 4. Fehlerangabe ohne Fehler von e. 


Vorliegende Arbeit wurde auf Anregung und unter Leitung von 
Herrn Prof. Dr. F. Harms im Physikalischen Institut der Universität 
Würzburg ausgeführt. Herrn Prof, Dr. F. Harms bin ich für sein 
stets förderndes Interesse und seine zahlreichen Ratschläge zu 
tiefstem Dank verpflichtet. 

Für die verwendeten Akkumulatorenbatterien, die einer Hoch- 
spannungsbatterie, welche Herrn Prof. Dr. F. Harms von der Not- 
gemeinschaft zur Verfügung gestellt war, entnommen waren, sage 
ich der Notgemeinschaft der Deutschen Wissenschaften meinen Dank. 


1) R. A. Millikan, Phys. Rev. 35. S. 1231. 1930. 


Würzburg, Physikalisches Institut. 
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Reibungsdispersion 
polarer Lésungen bei kurzen elektrischen Wellen u 


Von Wilhelm Müller 
(Mit 8 Figuren) 


1. Einleitung, Ziel der Untersuchung 
Die statisch gemessene Dielektrizitätskonstante (im folgenden 
abgekürzt: DK.) einer Lösung mit Dipolmolekülen setzt sich zu- 
sammen aus dem Orientierungseffekt der Dipole und dem Defor- 
_ mationseffekt der Moleküle. Letzterer ist allein maßgebend für den 
im sichtbaren Gebiet gemessenen Brechungsquotienten. Bei elek- 
_ trischen Schwingungen können von einer gewissen Frequenz ab die 
=. _ Dipolmolekiile infolge der Reibung an den Molekülen des Lösungs- 
mittels dem elektrischen Felde nicht mehr vollständig folgen, und 


lieh nur noch der Deformationseffekt übrig bleibt. Die DK. e muß also 
in einem bestimmten Frequenzgebiet von ihrem „statischen Wert“ & 
auf den „optischen Wert“ ¢,, der gleich dem (Quadrat des im sicht- 
baren Gebiet gemessenen Brechungsquotienten ist, herabsinken. 
i Die Dispersionsformel lautet nach Debye?): 


die Kreisfrequenz der elektrischen Schwingung bedeuten. 
Man kann Formel (1) auch schreiben: j 


4 ist die Wellenlänge der benutzten elektrischen Schwingung. 
2 


‘Fir 4=7' ist ¢= Demnach bedeutet 2’ diejenige Wellen- 


länge, bei welcher ¢ gerade das Mittel aus der statischen und der 
optischen DK, ist. Die Relaxationszeit r ist definiert als die Zeit, 


in welcher die Orientierung der Dipole nach Entfernung des elek- 
7° 
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trischen Feldes infolge der ungeordneten Wärmebewegung auf den | 
eten Teil zuriickgeht. 

In der vorliegenden Arbeit wird die Relaxationszeit für ein 
einfaches Dipolmolekül (Nitrobenzol) in einem dipolfreien Lésungs- = 
mittel (Shellöl) in ihrer Abhängigkeit von der Temperatur unter- 
sucht. Die benutzte Methode besteht darin, daß ¢ der Lösung bi } 
konstant gehaltener Wellenlänge (2 = 60cm) und bei verschiedenen i = 
(zwischen 0° und 40°C) gemessen Die Grenz 


den Formeln (2) und (3) den Wert von r ‘beveshnen. ; 
Das ist deswegen von Interesse, weil die Relaxationszeit sich — a4 
nach Debye aus dem Molekülradius und der Zähigkeit der Flüssig- = 


keit nach den folgenden Formeln berechnen lassen soll. 


& 


k die Boltzmannsche Konstante a, 37.1071 Erg/Grad), 
T die absolute Temperatur, 
n die Zähigkeit der Flüssigkeit in Poisen (Dyn. sec/cm?), 
a der Molekiilradius in Zentimeter. 
Die Formel (5) ist die Stokessche Formel für die Reibungs- 
konstante einer in der Flüssigkeit rotierenden makroskopischen Kugel __ a 
vom Radius a. 
Die Erfahrung hat gelehrt, und auch im folgenden werden wir 
es wieder bestätigt finden, daß die Theorie in dieser einfachen Form 
noch nicht richtig ist. Bei den rotierenden Bewegungen der Dipo- 
moleküle sind offenbar die Voraussetzungen der Stokesschn 
Formel (5) nicht erfüllt. Wir werden zunächst einmal so rechnen, 
als ob keine Assoziation der Dipolmoleküle vorhanden sei, d.h. als 
ob jedes Dipolmolekül sich frei in der Lösung drehe. Nach den 
Untersuchungen von van Arkel und Snoek?) soll ja Nitrobenzol 
keine oder nur sehr wenig Assoziation zeigen. Ferner werden wir 
annehmen, daß der Reibungswiderstand in der von dem Dipol- 
molekül mitgerissenen Flüssigkeit der Zähigkeit der Flüssigkeit 
entspricht, wie man sie bei der Bewegung eines makroskopischen 
Körpers experimentell beobachtet. Die Abweichung von der 
Stokesschen Formel erklären wir dadurch, daß bei der Be- 
wegung eines einzelnen Dipolmoleküls die umgebende Flüssigkeit 


ty nicht in derselben Weise mitgerissen wird, wie bei der Bewegung _ 

= eS einer aus vielen Molekülen gebildeten großen Kugel. Die Stokessche _ 
; - Formel setzt nämlich voraus, daß an der Oberfläche der rotierenden 
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Kugel die Bewegung der umgebenden Flüssigkeit ganz die gleiche ist 
wie die der Kugeloberfläche. Diese Voraussetzung lassen wir bei den 
rotierenden Molekülen also fallen. Wir nehmen an, daß die Flüssig- 
keit an der Oberfläche des Moleküls gleitet. Das Verhältnis der Ge- 
schwindigkeit der Flüssigkeit an der Oberfläche des Moleküls zu 
der Geschwindigkeit der Moleküloberfläche selbst werden wir den 
„Mitbewegungskoeffizienten“ & nennen. Er läßt sich aus den Werten 
von r und n berechnen. 

Um die von Debye gegebene Theorie der Reibungsdispersion 
zu prüfen, sind Messungen an verdünnten Lösungen von Dipol- 
substanzen in einem dipolfreien Lösungsmittel nötig. Es liegt 
bereits eine große Zahl von Untersuchungen an reinen unverdünnten 
Substanzen vor. Dagegen ist bisher nur eine brauchbare Arbeit 
erschienen, in der verdünnte Lösungen untersucht wurden. R. Luthi*) 
maß die DK. von Lösungen von Nitrobenzol in einem hochviskosen 
Mineralöl (Shell BL 3) in einem Wellenlängengebiet von 4 = 18 cm 
bis A= 318 m. Aus seinen Messungen geht hervor, daß in dem 
Bereich, in dem die DK. absinkt, die Dispersionsformel (1) giiltig ist. 
Wir legen sie deswegen den im folgenden ausgefiihrten Berech- 
nungen zugrunde. Luthi fand merkwiirdigerweise bei der Lösung 
des Nitrobenzols in Shellöl außer dem Gebiet starker Dispersion in 
der Gegend von 4=0,5m noch ein zweites Gebiet bei A un- 
gefähr 100 m, wo ein allerdings weit geringeres Absinken der DK. 
stattfindet. In der vorliegenden Arbeit beschäftigen wir uns nur 
mit dem Dispersionsgebiet bei den hohen Frequenzen. Unsere 
Messungen scheinen aber (vgl. S. 110) auch die Existenz des lang- 
welligen Dispersionsgebietes zu bestätigen. 


2. Methode der Messung der DK. 


Die Apparatur wird ausführlich beschrieben werden in einer 
demnächst erscheinenden Arbeit von E. Plötze. Wir bringen des- 
wegen hier nur das Wesentliche. Die Apparatur bestand aus einem 
Paralleldrahtsystem (Lecheranordnung) mit einer festen und einer 
verschiebbaren Brücke. Die bewegliche Brücke saß auf einem Prä- 
zisionsschlitten, der hydraulisch verschoben werden konnte. 

Der Sender bestand aus einer Osramspezialréhre in Barkhausen- 
kurzschaltung. Der Aufbau des Senders ist von W. Kroebel*) an- 
gegeben. Der Sender war äußerst lose mit dem Lechersystem gekoppelt. 
Als Indikator wurde ein von J. Brentano?) beschriebener Röhren- 
indikator verwendet. Er bestand aus vier zu einer Wheatstoneschen 
Brückenanordnung geschalteten Widerständen, von denen zwei durch 
Elektronenröhren gebildet wurden. Gitter und Kathode der einen 
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Fig. 1. Photographisch aufgenommene 
Resonanzkurve 


= 


instrument diente ein Deubnersches Schleifengalvanometer. Die 


Fig. 2. 


Flüssigkeitskondensator 


graphisch aufgenommen, die 
Mitte des Fadens wurde durch — 
einen engen Spalt abgegrenzt 
und auf die mit der Brücke fest 
verbundene Platte projiziert. 
Die Brücke bewegte sich zwischen 
zwei Anschlägen. Zu Anfang © 
und zu Ende der photogra- 
phischen Aufnahme lag sie kurze 
Zeit an je einem dieser An- ; 
schläge, so daß der Lichtschlitz — : 
hier als starker schwarzer 
Streifen erscheint. Fig. 1 zeigt 
das Positiv einer photographisch __ 
aufgenommenen Resonanzkurve. 
Die Entfernung von der Grenz- 
linie im Anfang bis zum Reso- ae 
nanzmaximum wurde mit einem 
Komparator gemessen. - 
Vor den DK.-Messungen wurde die 
Wellenlänge der benutzten Schwingung 
bestimmt. Dazu wurde die Entfernung 
zweier benachbarter Resonanzstellen ge- a 
messen. Diese Entfernung ergab sich 


zu = 30,0 cm. Also ist 2 = 60,0 cm. 


In einem Abstand von etwa °/,A 
von der festen Briicke befand sich me 20 
Flüssigkeitskondensator. Seine Form ist _ 
aus Fig. 2 ersichtlich. Das zylindrische Er 
Gefäß war von einem weiteren ange- __ 
schmolzenen Glasrohr umgeben. Der 
Raum dazwischen war evakuiert, damit 
der Flüssigkeitskondensator thermisch 
gut isolierte. Beiderseits waren diebeiden 
Platindrähte Pt, Pt eingeschmolzen. Ihre ox 


> 
inneren Enden hatten einen Abstand von © . 
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Flüssigkeit wurde in das untere kugelförmige Gefäß eingefüllt. Dieses 
befand sich in dem Olbad O. Das horizontale Rohr R konnte mit einem 
evakuierten Ballon verbunden werden, um die Flüssigkeit hochzusaugen. 
Um gleichmäßige Temperatur zu bekommen, wurde die Flüssigkeit da- 
durch gerührt, daß man sie abwechselnd hochsaugte und wieder in 
das Vorratsgefäß zurückfallen ließ. Die Temperatur wurde mit dem 
Thermoelement Th gemessen. Man konnte bis zu Temperaturen 
von + 2°C heruntergehen. Dann aber begann die Lösung sich zu 
entmischen; außerdem wurde sie so zähe, daß die eben beschriebene 
Methode des Rührens nicht mehr angewendet werden konnte. Auf 
Messungen bei tieferen Temperaturen als + 2°C mußte deshalb ver- 
zichtet werden. 

Um konstante Verhältnisse im Sender und im Indikator zu 
haben, wurden die Röhren ungefähr !/, Std. vor Beginn der Messung 
eingeschaltet. Während der Messung selber, die höchstens */, Std. 
dauerte, änderten sich dann die Betriebsdaten des Senders nur sehr 
wenig, wie Tab. 1 zeigt. Eine Kontrollmessung ergab, daß diese 
Änderung keinen merkbaren Einfluß auf die Resonanzlage, also auf 
die Frequenz des Senders hatte. 


Tabelle 1 


Vor der Messung Nach der Messung 


Gilterspannung .......... 138,0 Volt 137,6 Volt 

Die MeBapparatur wurde dadurch geeicht, daB man nacheinander 
zwei Flüssigkeiten mit bekannter DK. in das MeBgefiB brachte. 
Vorher wurde die Resonanzkurve in Luft aufgenommen. Gemessen 
wurde die Verschiebung der Resonanzstellung gegen Luft. 


Tabelle 2 


Verschiebung 
Substanz DK. der Resonanzstellung 
gegen Luft 
| 4,81%) 11,2 + 0,1 mm 


Die Eichkurve in Fig. 3 stellt die Abhängigkeit der Verschiebung 
der Resonanzstellung gegen Luft von der DK. dar. ee 
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10°/, Nitrobenzol in Shell K 20 


. 


Resonanzstelle | Verschiebung 
(gemessen von der der Resonanzstelle DK. 
linken Anse shlagsmarke) gegen Luft 
Luft . | 2 10 + 0,05 mm 0,0 — 
Lösung bei 17, 0°C 16,60 9. 2 + 0,1 mm 3,93 + 0,04 
Lösung bei 30,5° C 17.10 ; Se 4,14 
Lösung bei 39,00 C 17,05 > | 9,7 4,14 


In ähnlicher Weise wurde 
auch bei den Messungen jedes- 
mal zuerst die Resonanzstellung 
für Luft beobachtet. Darauf 
wurde die Flüssigkeit hoch- 
gesaugt. In Tab. 3 ist als Bei- 
spiel ein Messungsprotokoll vom 
15. 3. 34 angeführt. 

3. Messung der Zähigkeit 
Die Zihigkeiten der Lö- 
sungen und der reinen Lösungs- 
mittel wurden mit einem Ost- 
waldschen Viskosimeter er- 
| mittelt. Das Viskosimeter stand 
Wy erscmeoung gegen in einem Wasserbad, dessen 
a 50 100 mm Temperatur bis auf + 0,02° C 
Fig. 3 Eichkurve konstant gehalten wurde. B: 7 
4. Die verwendeten Substanzen 

Das Nitrobenzol verdankte ich der Freundlichkeit von Herrn 
Dr. R. Signer, der es im hiesigen Chemischen Laboratorium in 
reinster Form darstellte. 

Das verwendete Paraffinöl wurde von der Firma E. Merck, 
Darmstadt bezogen. Es ist im Merckschen Katalog mit: „Paraffin 
flüssig, nach Roth“, bezeichnet. 

Die Versuche ergaben, daß die Lösungen von Nitrobenzol in 
diesem Paraffinöl bei der verwendeten Frequenz noch keine Reibungs- 
dispersion zeigten. Deshalb ging ich zu einem zäheren Lösungs- 
mittel über, und zwar verwendete ich nach dem Vorgang von 
R. Luthi (a. a. O.) ein hochviskoses Mineralöl. Dieses wurde mir 
freundlicherweise von der Firma Rhenania-Ossag in Form des 
Shell-Kabelisolieröls K 20 zur Verfügung gestellt. Es ist ein klares 
Ol und besitzt bei 20° C eine Zähigkeit von etwa 60 Poisen. ae 
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er Lösungen wurden so hergestellt, daß man 1 Volumenteil 
der Dipolflüssigkeit mit 10 Volumenteilen des Lösungsmittels mischte. 


2 rk 5. Lösung von Nitrobenzol in Paraffinöl 


me Die DK. von reinem Paraffinöl wurde 12mal bei verschiedenen 
"Temperaturen gemessen, die DK. von einer Lösung des Nitrobenzols 
in Paraffinöl 9mal und von einer zweiten in genau gleicher Weise 
hergestellten Lösung 14mal. Die Resultate sind in Fig. 4 graphisch 


4 


€ 


{ 2-600cm 


1 l 1 


Fig. 4. DK. von Paraffinöl (IT) und einer Lösung von Nitrobenzol in Paraffinöl (I) 


dargestellt. Man sieht, daß in dem Temperaturintervall + 6° C bis 
+44° C die Änderung der DK. sowohl beim reinen Paraffinöl, 
als auch bei der Lösung noch innerhalb der Fehlergrenzen der 
Messungen liegt. Die DK. des reinen Paraffinöls ergibt sich als 
& = 2,18 + 0,04. Mit der Clausius-Mosottischen Formel ist 
dieses Resultat in Einklang. Nach ihr wäre in dem Temperatur- 
intervall infolge der Dichteänderung eine Änderung der DK. zu 
erwarten, die innerhalb der Fehlergrenzen liegt. Der Brechungs- 
exponent des Paraffinöls bei Natriumlicht wurde mit einem Abbeschen 
Refraktometer bestimmt zu 1,485. Daraus berechnet sich die optische 
DK. &, = 1,485? = 2,20. Das Paraffinöl ist eine dipolfreie Flüssig- 
keit. Deswegen erfüllt es die Bedingung &, = &. 

Die DK. der Lösung von 10 Volumenprozent Nitrobenzol in 
dem Paraffinöl ergibt sich zu &, = 4,16 + 0,04. Wenn bei der be- 
nutzten Frequenz schon eine Reibungsdispersion an der Lösung zu 
beobachten wäre, so müßte man bei verschiedenen Temperaturen 
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verschiedene Werte von & bekommen; denn die starke Änderung 
der Zähigkeit mit der Temperatur müßte zur Folge haben, daß man 
sich bei veränderter Temperatur auch an einer anderen Stelle der 
Dispersionskurve befiinde. Die Temperaturunabhängigkeit der DK. 
beweist also, daß wir bei der benutzten Frequenz noch außerhalb 
des Dispersionsgebietes sind. Der gefundene Wert der DK. ist 
demnach der statische Wert ¢,. Der Brechungsexponent der Lösung 
für Natriumlicht ergab 
sich zu 1,49. Daraus 
ergibt sich die optische 
DK. &, = 1,49? = 2,22. 
Für die Lösung besteht 
also zwischen «, und ¢, 
ein sehr erheblicher 
Unterschied. Es ist zu 
erwarten, daß bei höhe- 
ren Frequenzen als der 
von uns benutzten, Rei- 
bungsdispersion zu beob- 
achten wäre. 

In Fig. 5 ist die 
Zähigkeit sowohl des 
reinen Paraffinöls als 
auch der Lösung in ihrer 
Abhängigkeit von der 
Temperatur dargestellt. 

Die Zähigkeit der Lö- 
# % gung ist durchweg kleiner 
Fig. 5. Zähigkeit von Paraffinöl (7) und einer „]s die des reinen Paraf- 

Lösung von Nitrobenzol in Paraffinöl (IT) 


finöls. In beiden Flüssig- 
keiten ist die Änderung mit der Temperatur ganz außerordentlich 
groß. In der Lösung sinkt sie von einem Wert 1,6 bei +6°C 
auf 0,25 bei +36°C. Eine derartig starke Änderung hätte sich 
in der DK. bemerkbar machen müssen, wenn man im Gebiet der 
Reibungsdispersion wäre. 

Von der Lösung Nitrobenzol in Paraffinöl wurde die DK. noch 
bei einer Wellenlänge von 4 etwa 500 m von Herrn cand. phys. 
A. Riedinger nach der Schwebungsmethode gemessen. Hier ergab 
sich der Wert «, = 4,49, der wesentlich größer ist als der von mir 
gefundene Wert &, = 4,16. Ähnliches hat R. Luthi für die Lösungen 
von Nitrobenzol in Shellöl beobachtet. Es scheint demnach, dab 
für diese Lösungen schon bei verhältnismäßig langen Wellen ein 
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Dispersionsgebiet existiert, wo die DK. um einen geringen Betrag 
4,49—4,16 = 0,33 absinkt. Erst bei kürzeren Wellenlängen erfolgt 
dann der Abfall auf die optische DK. &,. Um diese beiden Dis- 
persionsgebiete zu erklären, werden weitere Untersuchungen 
nötig sein. 


us 6. Lösung von Nitrobenzol in Shellöl 


Mit reinem Shellöl wurden drei Messungen, mit einer Lösung 10 
und mit einer zweiten in genau gleicher Weise hergestellten Lösung 
8 Messungen bei verschiedenen Temperaturen zwischen + 2°C und 
+ 40° © ausgeführt. Die Resultate sind in Fig. 6 wiedergegeben. 


ag“ 


Do wr 
Temperatur 
Fig. 6. DK. von Shellöl (IT) und einer Lösung von Nitrobenzol in Shellöl (7) 


Während die DK. des reinen Shellöls sich mit der Temperatur 
nicht merklich ändert, zeigt die der Lösung einen starken Abfall 
nach den tiefen Temperaturen hin. Die DK. des reinen Shellöls 
beträgt 2,30. Dieser Wert stimmt sehr gut überein mit der 
optischen DK. Der Brechungsexponent bei Na-Licht wurde nämlich 
zu 1,510 bestimmt. Also ist ¢, = 1,51? = 2,28. Das benutzte Shellöl 
ist demnach wie das Paraffinöl dipolfrei. 

Um die Punkte der für die Lösung gefundenen Kurve aus- 
zuwerten, müssen wir die Grenzwerte ¢, und &, kennen. Wir dürfen 
jedenfalls, wie im Fall des Paraffinöls, beide als von der Temperatur 
nahezu unabhängig ansehen. Der optische Wert ist innerhalb der 
Fehlergrenzen gleich dem des reinen Shellöls, also e, = 2,30. Um 
den Wert e, abzuschätzen, nehmen wir an, daß die Differenz gegen 
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das Lösungsmittel Shellöl ebenso groß wie gegen das Lösungsmittel 
Paraffinöl ist. Da die DK. der beiden Lösungsmittel sehr wenig 
verschieden ist, und da wir beide Male dieselbe Konzentration ge- 
wählt haben, darf man diese Annahme mit ziemlicher Sicherheit 
machen. Beim Paraffinöl war die Differenz 4,16 — 2,18 = 1,98. 
Demnach berechnet sich für das Shellöl ¢, = 2,30 + 1,98 = 4,28. 
Die Kurve in Fig. 6 ist so gezeichnet, daß sie die beiden Geraden 
e=& und ¢ = «, asymptotisch berührt. Den folgenden Berechnungen 
sind als s-Werte die Punkte auf dieser Kurve zugrunde gelegt. 

Wir können nun für jede Temperatur den Wert von 4’ und 
damit den Wert von r berechnen. In Tab. 4 sind die so gewonnenen 
Werte zusammengestellt. 


Tabelle 4 


Temperatur in ° C € X cm T sec 


0 3,25 62,4 2,28 - 107" 
10 3,56 45,3 1,65 - 10-10 


3,84 32,1 117-107" 
30 4,08 21,2 0,77-10-% 


Wir berechnen weiter nach Formel (4) aus der Relaxationszeit r 
die scheinbare Zähigkeit 7’ der Lösung: 


-k-T 
Poisen . 


5 Der Radius a des Nitrobenzolmoleküls berechnet sich angenähert 
> aus dem spezifischen Gewicht des reinen Nitrobenzols s, dem Mole- 
u kulargewicht M und der Loschmidtschen Zahl N. Er ergibt 
sich zu 
3 — 
M 
= 2,8.10”® cm. 


Ferner wurde die wirkliche Zähigkeit 7 der Lösung bei den in Be- 
tracht kommenden Temperaturen direkt gemessen. Die gefundenen 
Werte sind in Fig. 7 dargestellt. In Tab. 5 sind die aus der Re- 


Er laxationszeit berechneten scheinbaren und die richtigen Werte der 
“ Zähigkeit der Lösung nebeneinander gestellt. 
Tabelle 5 
4 rm P n n 4 
= ratur in sec agg 7 
== = = Poisen = I = 
0 2,28 10-0 32 -10-% | 52 | 06-103 
10 1,65 - 1071 2,3 - 10? 18 
20 1,65 - 107? 7,4 2,2 - 10° 
30 0,77 - 10710 LE | 3,2 3,4 - 
= 
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Wie man aus dieser Tabelle sieht, unterscheiden sich die berech- 
neten Werte von den gemessenen Werten um rund 3 Zehnerpotenzen. 
Nach Formel (4) sollte außerdem die Relaxationszeit der Zähigkeit 
der Lösung proportional sein. Während aber die letztere in dem 
Temperaturintervall von + 30° C bis + 2°C auf das l4fache an- 


steigt, wächst die Re- 
laxationszeit nur bis 
auf das 21/,fache. Wir 
haben also anzunehmen, 
daß zwischen dem ge- 4, 
lösten Dipolmolekül und 

den umgebenden Mole- zy 
külen des Lösungsmittels 

eine Gleitung stattfindet. 4 
Wenn « das Verhältnis 
der Geschwindigkeit an 
der Oberfläche des Mole- 
küls zu der Geschwindig- 
keit der Molekülober- 
fläche selbst ist, so gilt 
für die Relaxationszeit 
die Formel 


90 


also nicht die Zähigkeit 
der Lösung selber, son- 
dern das Produkt =«:n, 
und der Quotient = ist nichts anderes als der Mitbewegungs- 


t-kT . 
gibt 


3 


koeffizient «. 
In der letzten Spalte der Tab. 5 ist der Mitbewegungskoeffi- 
zient & als Funktion der Temperatur angegeben und in Fig. 8 ge- 
zeichnet. Wie man sieht, steigt er mit zunehmender Temperatur 
stark an. Man kann das verstehen, wenn man bedenkt, daß bei 
höheren Temperaturen nicht nur das gelöste Dipolmolekül öfter und 
heftiger gegen die Moleküle des Lösungsmittels stößt, sondern 
außerdem der Zusammenhang der Moleküle des Lösungsmittels unter- 
einander auch gelockert wird, wie die starke Abnahme der Zähig- 
keit mit zunehmender Temperatur zeigt. 
x von der Lösung des Nitrobenzols in Shellöl sowie von 


/emperatur C 
Fig. 7. Zähigkeit von Shellöl (7) a 
und einer Lösung von Nitrobenzol in Shellöl (ID 


| hed 13 .. 
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ergaben sich die Werte: 


so erhalten wir fiir den Mitbewegungskoeffizienten bei 20° C den 
2 7 
Wert a = li 22.1073, Das ist etwa derselbe Wert, 
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dem reinen Shellöl wurde von A. Riedinger die DK. bei einer 
Wellenlänge von A etwa 500m (yn = 6 - 10°. sec”!) gemessen, Es 


é= 2,40, 
= 4,85. 
reines Shellöl stimmt mit dem bei kurzen Wellen 


— liegt der 


reines Shellöl 
Nitrobenzollösung & 


Der Wert für 
gefundenen Wert 2,30 einigermaßen überein. 


= a 


Fig. 8. Der Mitbewegungskoeffizient von Shellöl mit Nitrobenzolmolekeln 


Wert der Lösung auch hier wieder weit über dem (Grenzwert 
& = 4,28, den meine Messungen ergeben. Der Unterschied ist 0,57. 
Diese Differenz ist ungefähr von der Größe des Sprunges, den Luthi 
für das langwellige Dispersionsgebiet bei seinen Lösungen findet. 


7. Vergleich mit den Ergebnissen anderer Untersuchungen 


Wir vergleichen unsere Ergebnisse zuniichst mit denen von 
R. Luthi (a. a. O.). Luthi bestimmte im kurzwelligen Dispersions- 


gebiet bei der Lösung von 5 Gew.-°/, Nitrobenzol in Shell BL 3 
die Relaxationszeit des Nitrobenzolmolekiils zu r = 1,64 1071 sec 
bei 20°C. Aus diesem Wert berechnete er die Zähigkeit der 
Lösung und fand n= 0,025 Poisen. Die Zähigkeit des reinen 
Shellöls ist » = 23 Poisen bei 20°C. Die Zähigkeit der Lösung 
von 5 Gew.-°/, läßt sich zu etwa 1/, bis !/, des Wertes von dem 
reinen Lösungsmittel schätzen. Rechnen wir zB. mit „= 11,5, 
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den unsere Lösung von 10 Vol,-°/, in dem Shell-Kabelisolieröl 
K 20 bei 20° C zeigt. 

Nach einer ganz anderen Methode hat neuerdings P. Debye‘) 
die Relaxationszeit von Dipolmolekülen in Lösung bestimmt, nämlich 
aus der infolge der Absorption der Energie der elektrischen 


_ Sehwingung auftretenden Erwärmung der Lösung, die nach dem von 


J. Malsch®) angegebenen Verfahren ermittelt wurde. Unter den im 
Debyeschen Institut gemessenen Lösungen befinden sich auch Nitro- 
benzollösungen von der Konzentration 3 Vol.-°/, in einer Reihe von 
dipolfreien Lösungsmitteln von verhältnismäßig sehr geringer Zähigkeit. 
In der Tab. 6 sind die von Debye für Nitrobenzol gefundenen Werte 
von r angegeben. In der letzten Kolonne habe ich die daraus von 
mir berechneten Werte des Mitbewegungskoeffizienten « hinzugefügt. 


Tabelle 6 


Zähigkeit Relaxationszeit Mitbewegungs- 


Lösungsmittel 
8 sec koeffizient « 


326 1,9 10712 8,6 107° 
4,4102 9,8 - 107% 
 Cyelohexan . . . 2,8. 1071? 5,0 10? 
Dekelia. .- . 5,0 - 107? | 3,0 + 1072 


Die Mitbewegungskoeffizienten in diesen Lösungsmitteln von 
geringer Zähigkeit und von kleinerer Molekülgröße sind also durch- 
weg viel größer als der in der Shellöllösung gefundene Wert bei der 
gleichen Temperatur. Man kann das verstehen, weil die kleineren 
und leichter beweglichen Moleküle der von Debye gewählten Lö- 
sungsmittel natürlich leichter mitgerissen werden als die großen 
ziemlich fest miteinander zusammenhängenden Moleküle des Shellöls. 
Der für Benzol als Lösungsmittel gefundene Wert des Mitbewegungs- 
_koeffizienten ist besonders groß. Das deutet darauf hin, daß die 
Benzolmolekiile an den Dipolmolekiilen besonders fest haften. 
Dafür sprechen auch andere Erscheinungen nach J. Snoek"), 


Zusammenfassung 
1. Bei konstanter Wellenlänge 2 = 60,0 cm (v = 5 10%. sec') 
wurde nach der zweiten Drudeschen Methode die DK. von fol- 
genden Flüssigkeiten bei verschiedenen Temperaturen im Intervall 
+ 2°C bis etwa + 40° C bestimmt. 
a) Von einer Lösung von 10 Vol.-°/, Nitrobenzol in Paraffinöl. 
b) Von einer Lösung von 10 Vol.-°/, Nitrobenzol in Shell- 
Kabelisolieröl K 20. 
c) Von den beiden reinen Lösungsmitteln. 
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2. Während die DK. der reinen Lösungsmittel und auch der 
Nitrobenzollösung in Paraffinöl innerhalb der Fehlergrenzen sich als 
von der Temperatur unabhängig erweisen, zeigt die Nitrobenzol- 
lösung in Shell-Kabelisolieröl K 20 einen starken Gang mit der 
Temperatur, der darauf hinweist, daß die verwendete Frequenz für 
diese Lösung bei Temperaturen zwischen 2° und 40° C im Gebiet der 
Reibungsdispersion liegt. 

3. Aus den gemessenen Werten wird die Relaxationszeit des 
Nitrobenzolmoleküls in der Shellöllösung und der Mitbewegungs- 
koeffizient & des Lösungsmittels mit dem Nitrobenzolmolekül als 
Funktion der Temperatur berechnet. Es zeigt sich, daß der Mit- 
bewegungskoeffizient mit der Temperatur stark ansteigt. 

4. Messungen der DK. bei wesentlich längeren Wellen 4 ~ 500 m 
(v = 6-105 sec”') an denselben Flüssigkeiten scheinen den von 
R. Luthi gemachten Befund zu bestätigen, daß die Lösungen 
des Nitrobenzols in einem hochviskosen Lösungsmittel außer dem 
Dispersionsgebiet bei kurzen Wellen noch ein zweites bei wesent- 
lich längeren Wellen haben. Die Änderung der DK. in diesem 
langwelligen Dispersionsgebiet ist aber jedenfalls wesentlich kleiner 
als in dem Hauptgebiet der Dispersion bei den kurzen Wellen. 


Es ist mir eine angenehme Pflicht, meinem hochverehrten 
Lehrer Herrn Geh. Rat Prof. Dr. G. Mie für die Anregung zu dieser 
Arbeit, sowie für viele wertvolle Ratschläge bei ihrer Ausführung 
meinen herzlichsten Dank auszusprechen. Herrn Dr. E. Plötze 
danke ich herzlich für seine Hilfe bei der Ausführung der Messungen. 
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